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Résumé. — L’agitation thermique développe des forces de rappel qui n’obéissent pas à la loi 
de Hooke. Elle n’est pas analysable en oscillations harmoniques indépendantes, au nombre de 39€ 
dans un cristal formé de 2T atomes. Néanmoins, elle peut être résolue en composantes harmo- 
niques qui comprennent des oscillations fondamentales, au nombre de 391, et des oscillations 
secondaires qui sont au sens large les harmoniques des oscillations fondamentales. Les oscillations 

ÿ fondamentales et secondaires développent des forces sensiblement constantes au cours du temps 
qui constituent les tensions thermiques. 


Abstract. — Thermal agitation produces interatomic forces which are notin concordance with 
Hooke’s law. It cannot be analyzed into independant harmonic oscillations, which are 39 in 
number in a crystal composed of #t atoms. Nevertheless, it can be resolved into harmonic com- 
ponents which include fundamental oscillations, to the number of 3 AT, and secondary oscillations 
which in the broad sense are the harmonics of fundamental oscillations. Fundamental and 
secondary oscillations develop forces sensibly constant in course of time which constitule thermal 
strains. 


Je fixerai l’origine des coordonnées sur la posi- positions des atomes dans le cristal en cet état 
tion moyenne d’un atome et prendrai des axes de seront définies par les translations 
coordonnées rectangulaires qui forment des angles FES 
: 2 At: E Mo + Jo; Po + k,; 
constants avec trois rangées réticulaires (non copla- 
naires) du milieu cristallin libre de toute con- les déplacements entraînés par la dilatation qui a 
trainte, à une température constante et uniforme. lieu entre l’état inerte et l’état atteint à la tempé- 
Les positions des atomes seront définies par trois  rature absolue @ étant : 
translations : 5 ; à à 
#j,0 = M +] —Mo—Jo; 
m+j+uw; p+rkh+uw; q+il+uw;..…. 


M, P, 4, ..., étant des translations du réseau cris- 
talhn, j, k&, 1, ..., (1) des vecteurs inscrits dans la 
maille, u?, u?, uf, ..., les élongations des oscilla- 
tions faites par les atomes situés sur les positions 
moyennes M + j, p + k, q + 1, ...; la masse 
d’un atome en position j dans la maille sera désignée 
par. u;, la masse d’un atome en position & par 
W, -.., et la masse d’un motif cristallin par pu: 


to =p+k—p,—k: ... 


Il sera entendu que le milieu cristallin est en 
équilibre thermique, c’est-à-dire : à température 
constante et uniforme, et que toutes les grandeurs : 
les énergies et leurs dérivées, les translations m, 
les élongations u*, les dilatations »$, etc., se 
rapportent à la température donnée. 

L'énergie potentielle, %, d’un cristal libre de 
toute contrainte, à la température ©, peut être 
u = 2 Uj. exprimée par une série [1-2] 


m 
J’appellerai atome mj celui qui se trouve sur la  V = We, + > 2 D Bin 
position m + j, atome pk celui qui se trouve sur la 


| position p + k, ... Et j'envisagerai le milieu re ÿ SAS CT £ a 20 
cristallin en état inerte, soit à la température du 2 mp 5e «6 
"ZE ’agitati éteinte. Les " 4 
zéro absolu et l’agitation thermique étei os, D’ ” ue nf 


e , re à 3 | mpa jkt opy 
(1) L’origine étant sur la position moyenne d’un atome 


l’un des vecteurs j est nul. JE 
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2 
a, , y IS (1) 
où 
r lo , 1 ) 
We, ZDATE + D'DUPET + (2) 
2:mD jk 9 : mpa ik 
(m+jÆp+k) (m+jÆp+k#q+lm+)) 
à — D (6) 
BED C 
DOTE 
1 mpa DB (Y 
D DE DTA CT ee € 
DL gai Br La TEEN EAN NE S) 
D Vogue HARTE 
RE ENOR f+222D JE ri. (4) 
mpa mpq mpar ,Ù 
VUE jkl j 7 
DID A 2 2 > Es ke + (5) 


et W, étant l’énergie potentielle du cristal en état 
inerte : 


C Vue: Se DAC 
10/@ 06 Der mx 2b m . D \ 
dti dTk (Æÿ = M6 + jo; 2x = P + k) 


(6) 


M 
CE) 


De 
De | 
% el 


2 = Po + ko; Lie do + Lo) (7) 


= Mo + Jo ; 


Han FA [ > Wo | 
4 da%" . (&ÿ =m +j; x = p+k) (6 bis) 


mra = 
jk 
Dir} x: | 


d We | 
da?” .. da (a? =m+j:aË =p+k: 


a —q +) (7 bis) 

He PE RATER LS el ee). ù 
SCI f-ECTÉE-0 (8) 
SD-SD-SDE 0 ( 


I. Les forces et les tensions développées par la 
dilatation thermique. — Au cours de la dilatation 
thermique qui a lieu entre l’état inerte et l’état 
pris à la température © par le milieu eristallin 
libre de toute contrainte, deux atomes subissent un 
déplacement relatif 


8 ! 
m 7. < : À 
5,0 KO —= = xxP k > ADN ES Les Po — kb) 


(10) 


La dilatation thermique conserve la structure 
périodique du milieu cristallin, Chacune de ses 
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composantes 4$ prend la même grandeur-en toute 
parie du cristal. Mais si le motif cristallin est 
polyatomique, la dilatation rest pas, en général, 
strictement linéaire : elle ne se réduit pas aux 
déplacements relatifs 


> 18 (60 + jo — po — ko) ; 


les réseaux 7, k, ..., formés par les atomes en 
mêmes positions j, k, dans la maille, se déplacent 
les uns par rapport aux autres : À — A" est le 
déplacement relatif des deux réseaux 7 et X. 

Comme toute déformation, la dilatation déve- 
loppe des forces de rappel entre les atomes. Celles 
qui sont appliquées sur un même atome, M], exer- 


, D, 
cées les unes f45 # par un seul atome pk, les 
œ 


autres far M par deux atomes distincts pk et 
ql (@), ..., ont une résultante : 


M—D 


pen ie 


M—DI 
DEPEBLLE" DMC 
? pa kl 


LR É ms 
Dai Ja = m + ÿ 


et compte tenu des relations (3), (8), (9) et (10) 


(11) 


I 


M—D 


Nr AR A PEN B 
hi ZE CE ÉLUS — 0) 


+ > 188 + j5 — po8 — k8)] 


Le 2 D sd LO — x) 
a 2 1$ (mod 7,8 Po kGÿ)] 
x LOT — 87) + D mot + Jot — pot — KoË)] 
+) 
HDi — 2 DHL __ x) 
SE 2 18 (m8 Æ 100 — po — K)] 


x LOT — 2) EE mat 2 jo — quf — do) 


m PR ; 
La force — B; se répète sur tous les atomes en 
même position dans la maille : 


My ns: 


Supposons chaque atome sur un centre de symé- 
trie (cristaux cubiques formés par les métaux, sel 
gemme, etc.….). À tout atome, pÆ correspond un 
atome p'k, à tout atome gl un atome g'l, ..., sur 
des positions telles que 


m +j—?p E= (m + j 


(?) Les forces jai #1 dérivent des énergies potentielles 


triatomiques WR1 (2). 


No1 
eb 
mp M—p mp’. mp 4 
CEE ÿ k. NOTES ik. 
0& 8 Coë ( De De 


. En outre, quand le milieu cristallin se dilate, 
libre de toute contrainte, chaque atome reste sur 
un centre de symétrie. La dilatation thermique est 
strictement linéaire (A — A# — O). 

En conséquence, les forces far +; ai ti (42-15), 
exercées sur un même atome, M", sont deux à deux 
égales et opposées : 


M—p M—D m—p"q" M—PA 
MR == k: LRÉRUNEREEN FRERE 


leur résultante (11) est nulle : 
sp: = 0. 


L'énergie potentielle W,4, (1), rapportée aux posi- 
tions moyennes des atomes, a ses dérivées pre- 
mières nulles (11-14). Elle est donc minimum ou 
maximum selon que l’édifice formé par les atomes 
en positions moyennes est en équilibre stable ou 
instable sous les tensions thermiques qu’il supporte. 
Si l’équilibre était instable, il serait rompu par 
l’agitation thermique. Ainsi, les égalités (14) 
expriment la condition de stabilité du milieu cris- 
tallin soumis aux tensions thermiques. 

Chez la plupart des espèces cristallines les atomes 
ne sont pas situés ou ne sont situés qu’en partie sur 
des centres de symétrie. De ce fait, la symétrie 
du milieu cristallin n'implique plus les égalités (14). 
Mais, si les forces —- B° (11) n'étaient pas nulles, le 
milieu cristallin serait instable. Les égalités (14) 
restent donc satisfaites quelles que soient la struc- 
ture atomique et la symétrie du milieu cristallin. 

Mais, de même que toute déformation, la dila- 
tation thermique développe des tensions 1in- 
ternes T,Y («, y = 1, 2,3). Si West la densité, par 
volume unitaire, de l’énergie prise par la dilatation 

thermique 


(14) 


Les tensions thermiques peuvent encore être 
définies de la sorte. Envisageons le cristal divisé par 
un plan P en deux parties À et B. La tension, Tiï 4m, 

—> 


dirigée suivant l’axe de coordonnées O1*, exercée 
re 

dans le plan P normal à l’axe Oa* par la partie À 

du cristal sur l’autre partie B, est égale à la résul- 

# A 1 el . . M—D 

tante, rapportée à l’aire unitaire, des forces fi k, 


1 , ..., exercées par les atomes, ph, ql, ..., 


situés dans la partie À, sur les atomes, My, situés 
dans la partie B [3]. Cette seconde définition ne fait 
pas état des forces far tt qui sont exercées sur un 


atome, mj, contenu dans la partie B, par deux 
atomes pk et gl situés l’un dans la partie À l’autre 
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dans la partie 2. Mais la contribution que de telles 
forces apportent aux tensions thermiques reste 
toujours minime. Nous avons : 


Don Ta Dee 
les tensions partielles T,%, TX ., étant respec- 


GE QC 
tivement produites par les forces fy7 k, fy %1,... 
(04 (02 


IT. Les forces développées par l’agitation ther- 
mique. — Ainsi, les dérivées premières B° (11) de 

rs . . ee 
l'énergie potentielle restent nulles à toute tempé- 
rature (14). Sur chaque atome, mj, est exercée une 
force f}, développée par l'agitation thermique : 


7 D 
Era à 
m 
du; 
mp 6 pa 6 dt 
=—222Ci fu EEZDiiuut—. 
DIBTLNEL f Ÿ 2 loa ir BY CAIERE 


(17) 
Le principe de d’Alembert donne les équations : 


us + X CT Eu 07 (18) 
où 
œ, se pa kl 6 Dik me “ii 
a kl RY 
DD D Prune RS) 


3 l'or kis By 


Le milieu cristallin et son champ de force ont une 
structure périodique, et les atomes exercent les 
uns sur les autres des forces de rappel qui décrois- 
sent rapidement quand leur distance augmente ; le 
rayon d’action sensible de ces forces est borné. En 
conséquence les oscillations atomiques du milieu 
cristallin forment des trains d’ondes planes [4]. 
Il convient donc de rechercher les solutions des 
équations (18) qui représentent des oscillations pro- 
gressives transportées par de telles ondes. 

Pour fixer les idées, envisageons un cristal, de 
masse M, comptant 9T atomes, répartis dans N 
motifs cristallins formés de g atomes 


INT A = Ne. 
III. Première analyse. — Les théories classiques 


de l’agitation thermique font seulement état des 
forces de rappel conformes à la loi de Hooke : 


m nt 
es AU 20 
fa ei Gi CRC (20) 


C'est négliger le second nombre, D (19) des 


équations de d’Alembert (18). 
Les équations ainsi réduites, 


.m MD 8 
à A p 9/ 
He A CES LE, (21) 
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admettent la solution : 
= jee \ 
ui} = DU pe. (22) 
ST 
où (°) 


1 fie à à 
WT 87 — 3 Asr VE LU se exp.i 27 (ver t + 057 — Sm) 
1 


* 
LU sr exp.— i 2r( vs t + 057 — Sm)) 


= A$r ci sr cos 27T(Vsr t + Ûs- + Dia, Sr — Sm) 


(23) 

1 2W$s- 
A$r de. mire OST T VS ( ) 
ge = PJ E : (25) 


Us” 


d* est le module et 2x9. l'argument de C5, 
c’est-à-dire : 


UF gr = Pla.sr EXP. 2 TPja,Sr, 
et 


* , : j 
2 sea Cas = D pese çasr = 1, (26) 
ja jx 


donc : 
Z ule%5+ cdsr) = u. 
F 104 


Les oscillations u,%s (23), de même fré- 
quence vs, pilotées par le même vecteur d’onde $, 
faites par tous les atomes du cristal, constituent 
une oscillation harmonique fondamentale, S+, de 
l'agitation thermique. 4, est l’énergie quantifiée 
de cette oscillation, 44, son amplitude moyenne 
prise sur les élongations de tous les atomes et sur 


. . Ver . M4 
toutes les directions de ces élongations, et cs Asr . 


mesure l'amplitude suivant l’axe de coor- 


— 
données O2 des oscillations u,%,. faites par les 
atomes en position j dans la maille. 

Je fais état des vecteurs d’onde fondamentaux, 
ceux qui ont les plus petits modules. En même 
nombre, NV, que les motifs du cristal, deux à deux 
égaux et opposés, ils sont inscriptibles dans une 
première zone de L. Brillouin, l’origine au centre 
de la zone, les extrémités formant un réseau de den- 
sité uniforme. Chaque vecteur d’onde fonda- 
mental S pilote 3g oscillations fondamentales. Donc 
dans la somme (22), S pris comme indice court 
de 1 à V, et l'indice +, désignant l’une des 3g oscil- 
lations pilotées par le vecteur $, court de 1 à 3g. 
Cela fait en tout 3 9T oscillations fondamentales ({). 


* ‘ d 
(5) Us est le complexe conjugué de K*5z. s 
(*) Abstraction faite de la structure périodique du milieu 
cristallin, les pulsations, w, et les coordonnées vectorielles, 


5”, des oscillations fondamentales sont déterminées par 
les équations séculaires 


. MD | ——— np 
Gi = Z Ci k]Vh; ue @ 
œ KB a & 


donc par les valeurs caractéristiques et les vecteurs propres 


No 


Les 3g composantes ©”. des g vecteurs C/s- qui 
se rapportent à la même oscillation harmonique, 
peuvent être prises pour celles d’un vecteur. 
unique, ©, représenté par une matrice C 
d'ordre 3g x 1. Le vecteur d’onde $ donné, les 
carrés w? des pulsations et les vecteurs & sont 
égaux respectivement aux valeurs caractéristiques, 
et aux vecteurs propres, normés à l’unité (26) 
d’une matrice T(S) (dite de Fourier), hermitienne, 
définie positive (°), d'ordre 3g X 3g: | 


wo? 6 = T(S)C. (27) 
La matrice l'(S) a pour éléments : 
Ci *€ 
DÉS Coxp.i2x S(m — p) (28) 
! ? Vu Ux 


j et k pris comme indices courent de { à g; 
Mali 120. 

T, w? et € sont des fonctions périodiques du 
vecteur d'onde S, et leurs périodes sont celles du 
réseau polaire. : 

En outre 


T(— S) = T*(S) 


de la matrice dynamique, définie par Max Born [5], 
d’ordre 39T X 39L, d'éléments : 


MD | ,—— 
Ci k [Vu Ux. 


Et le nombre des oscillations fondamentales est donné 
par le rang de la matrice dynamique. Il est donc déterminé 
par le nombre des relations linéaires entre les coefficients de 


rappel (À £ (6 bis). Car toute relation linéaire entre ces 


coeflicients entraîne une relation linéaire entre les colonnes 
de la matrice dynamique. 


D’après leur définition, les coefficients de rappel C5 £ 


sont liés par trois relations linéaires (8), imposées par 
l'égalité de l’action et de la réaction. Les oscillations fonda- 
mentales sont donc, en général, au nombre de 3297 — 5, 
Toutefois, si le milieu cristallin à une structure atomique 
simple et une symétrie élevée, les coefficients de rappel 
peuvent satisfaire trois autres relations linéaires, décelées 
par Max Born [6]. En ce cas le nombre des oscillations 
fondamentales tombe à 39T — 6. Mais, même si le cristal 
est minuscule, 3 et 6 sont négligeables par rapport à 321. 

Selon une conception courante, les oscillations fonda- 
mentales seraient toujours au nombre de 39T — 6 : un 
cristal formé deT atomes est un système à 3 degrés de 
liberté, et six de ces degrés seraient pris par les 3 trans- 
lations du milieu cristallin suivant les axes de coordonnées 
et par les 3 rotations autour des mêmes axes. Il ne resterait 
donc que 3T — 6 « degrés d’oscillation ». Mais la relation 
entre le nombre de degrés de liberté et celui des oscillations 
fondamentales prête à confusion. Une translation, une 
rotation, faite par le cristal en bloc, libre de toute con- 
trainte, ne développe nulle force de rappel entre les 
atomes [2]. On ne peut donc définir des forces de rappel 
qui seraient produites par les translations et les rotations 
du cristal, puis, en posant ces forces nulles, obtenir six 
relations linéaires satisfaites par les coefficients de 


rappel CruËe 


(5) Le milieu cristallin supposé stable. 
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donc : fréquences sont lv + vel, leurs vecteurs 
&(— S) = w!(S) (29) ou + (S + si oO Soit u,%" la résultante de 
ne. outes ces oscillations mê 
eb si S' = — Set r! — +, nous avons : (& — 1,23) M 
7 * 
+ La troisième opération consiste en la résolution 
Vs = Ver;  Cem = Ur. (30) des équations : 


Les oscillations harmoniques pilotées par deux 
vecteurs d’onde égaux et opposés ont deux à deux la 
même fréquence et la même forme. 

Les vecteurs d'onde, S, conformes aux conditions 
cycliques (Max Born) satisfont les relations : 


D exp.i2r Sm — 10 
m 


D exp.&27(S — S’) m — Nos. (31) 

Ces relations supposent la structure du milieu 
cristallin, rigoureusement périodique, partout, 
jusqu’à la surface du milieu. Cette structure n’est 
plus périodique dans l’assise périphérique du cris- 
tal. Néanmoins, négligeant la perturbation des 
oscillations fondamentales dans cette assise je 
ferai état des relations (31). 


IV. Développement de l’analyse. — Pour déve- 
lopper l’analyse harmonique de l'agitation ther- 
mique, je résoudrai les oscillations faites par les 
atomes, et les forces qu’elles développent, en com- 
posantes dont les expressions sont des poly- 
nômes homogènes par rapport aux amplitudes 
moyennes A4. (24) des oscillations fondamen- 
tales (23) ; et je classerai ces composantes selon 
l’ordre de leur degré croissant par rapport aux 
mêmes amplitudes. Afin d'éviter de le répéter à 
chaque fois, il sera entendu une bonne fois pour 
toutes, que le degré de toute grandeur prise en 
compte : élongation, force de rappel, etc. est, sauf 
indication contraire celui qui concerne les ampli- 
tudes A.. 

La force ® (19) est précisément analysable en 


forces On À » de degrés n distincts par rapport aux 


amplitudes A4. : ñ — 2, 3, .... Et l'analyse har- 
monique de l'agitation thermique peut être déve- 
_loppée par des opérations successives. 
La première opération est la résolution des équa- 
tions de d’Alembert sans second membre (21). 
En seconde opération, on porte dans le deuxième 
membre la force 


mm al 
Dors > 2 D 


: pa 8 


q 6 Y 
lux WU. (32) 


î 
œ 


2 


Les équations 
… m M—D  pG m 
ui +2 Ci kUy —=Pri 


harmo- 
leurs 


déterminent de nouvelles oscillations 
niques que je qualifierai de secondaires ; 


.. .m M—D 
bi Us à + à Ci, & up = Ps: (33) 
où 
7m MPa DB, dv 
DR di 
1 mpoar (e Y 5 
te re, kl ? q v 


La solution des équations (33) donne de nou- 
velles oscillations secondaires : elles ont les fré- 
quences [vs + Vs + vw-| et les vecteurs d’onde 
+ (S + S" + S”) ; toutes celles qui sont effectuées 
par un même atome, mnJ, forment une oscillation 
rlobale uw a(ee="1,2,3); 

Et ainsi de suite : on porte successivement dans 
le second membre des équations de d’Alembert (18) 
les forces Pni par ordre de degré n croissant, et l’on 


résout successivement les équations 


. M m—D LB m 
p 
HAT ER os 


m . » 
La force ®,; est exprimée par la somme de tous 
œ 
les termes de degré n qui sont contenus dans la 
série D? (19) où 
œ 
8 œ 
Uÿ = Uk +... + Uni, 
et u,_1%" est la résultante de toutes les oscillations 
secondaires, de degré r7 — 1, déterminées par 
l'équation 
m 


_ m M—D pê k 
bi Una t + 2 CE Elnnak = Quai 


Les équations de d’Alembert (18), 
M MD m m 
ju À + 2 Ci Eut =péitosi+..: = 
_ 


sont ainsi satisfaites par la série : 
m* mn mn m 
HIS Une 1 Me ni Use 


Toute force w,7 peut encore être résolue en 
œ 
forces élémentaires, variables dans le temps selon la 


(5) Si la fréquence est vs. + vs, le vecteur d’onde 
est S + S’;si la différence vs, — vs est positive, le 
vecteur d’onde est :S$ —— S’ ; si la différence vs: — Vs” 
est négative (fréquence vs — vs.), le vecteur d’onde 
est — (S — S’) == S” 


ee 13, 
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loi sinusoidale, de fréquences [vs + vs + ...|. 
Et ÿ’appellerai ordre d’une force élémentaire le 
nombre { des termes contenus dans la somme qui 
définit sa fréquence, vs: + vs + ...,où le même 
terme vs, Vs, -.-., peut se répéter plusieurs fois. 
Les forces élémentaires qui appartiennent à Ja 
force On 5 sont d'ordres 


m 


r Q Û m { 
loutefois, les forces Pris Pais 


pairs, comprennent des forces élémentaires d’une 
autre sorte. Ce sont des forces qui restent cons- 
tantes pendant des temps de durées limitées et 
passent instantanément d’une intensité à une 
autre ; cette variation discontinue de leur intensité 
résulte des échanges d’énergie, par quantum Ay., 
entre les oscillations fondamentales. Toutes les 
forces élémentaires de cette sorte (1 — 0), qui sont 
appliquées sur un même atome se compensent 
exactement a chaque instant (7) [9]. De telles 
forces n’entretiennent pas des oscillations ; elles 
constituent les tensions thermiques. Il convient 
donc de les éliminer pour faire l’analyse harmonique 
de lPagitation thermique. 

Ce sont les forces élémentaires périodiques, 
d'ordre 2, 5, , qui produisent les oscillations 
secondaires : toutes les forces élémentaires, de fré- 
quences [vs + vs = ...|, égales, entretiennent 
une oscillation secondaire de la même fréquence. 

Mais toutes les forces élémentaires du même 


, de degrés 


ordre ! ne sont pas du même degré ; elles font 
partie de forces En? différentes : 
DRE PIRE EUR 
m 
Par exemple, toutes les forces 927, çs s de 


degrés pairs, comprennent des forces élémentaires 
d'ordre 2. I] en résulte que toute oscillation secon- 
daire peut être résolue en composantes de différents 
degrés. Si l’on définit aussi l’ordre d’une oscillation 
secondaire par le nombre / des termes inclus dans 
la somme, ve + Vs + qui donne sa fré- 
quence, les composantes des oscillations d’ordre / 
sont de degré 


EU AE AE PR OS OU 


ê5 , de 

degrés impairs, comprennent des forces élémen- 
m Q 2 

taires fi, , À ., qui ont les fréquences vs. 


Se 


des ue fondamentales. 


Fe 3 m 
En particulier, les forces CRE Ps) 


cu } Comme se compensent les forces élémentaires déve- 
loppées par la dilatation (14). 


N°1 


Une oscillation fondamentale S+ est, en fait, 


m 
entretenue par des forces de rappel, Ji ç. COM- 
plexes : à 

/ on 
Î j 9Fr 
ee 2 TO “has ee FRS (35) 
"mn f 
j 5 in r 36 
" &, ST Ro œ AIS STE { i) 
. FAUOUTS nr 
Par l'intermédiaire de ses composantes, f35 
mn , 
re ds la force f5 … est fonction de toutes les 
œ DT % ,DT 
; m—p 
dérivées de l’énergie potentielle (1-2): Ci 
DE, ..., (hormis les dérivées premières qui sont 


nulles), et des coordonnées vectorielles &”s. et des 
amplitudes 44, de toutes les oscillations fonda- 
mentales (23). 

De ce fait, les oscillations fondamentales ne 
peuvent être déterminées qu’au terme de l’analyse 
et par des approximations successives. Pour définir 


les forces Î35 er Î5 tete il est nécessaire d’expri- 


mer les ee acoue secondaires en fonction des 
oscillations fondamentales. Et, en éliminant le 
temps des équations 


. m mn 
MU sr — Tose 
on aboutit à des équations séculaires, semblables 
aux équations (27), mais où les coeflicients des 
inconnues, &Ÿ, renferment comme paramètres les 
coordonnées vectorielles ©”. et les amplitudes 4%, 
de toutes les oscillations fondamentales. Enfin, on 
détermine les oscillations secondaires en portant 
dans leurs expressions les fréquences w., les coor- 


données vectorielles €”. et les amplitudes 4,4. des 
oscillations fondamentales. 

Afin d’alléger les calculs, il convient de procéder 
par analyses successives, de précision croissante. La 
première analyse est celle des théories classiques 
(équations (21)). La seconde a pour objet la réso- 
lution des équations 


ui + DC EUË = qi tou. (7) 


La troisième consiste en la résolution des équa- 
tions 


..m M—D @ m M m m 
x À : 1] 12 + . Q * 
mur + DOI EU qi toi tqar test (88) 


et ainsi de suite, jusqu’à ce que soit atteinte la 
précision exigée. 


Manuscrit reçu le 14 octobre 1958. 
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L'ABSORPTION ET LA DISPERSION ULTRASONORES 
DANS LES GAZ MONOATOMIQUES : 
LES TROIS SOURCES D'IRRÉVERSIBILITÉ CLASSIQUES 


Par Mme M. HUETZ-AUBERT et M. J. HUETZ 


Résumé. — L'étude des effets de dispersion et d’absorption dits classiques ou encore de trans- 
lation, puisqu'ils affectent ces degrés de liberté moléculaire, est utile pour séparer, des effets 
globaux seuls accessibles à l'expérience, ceux qui sont dus à la relaxation intramoléculaire. Comme 
cette dernière ne peut affecter le comportement des monoatomiques, le contrôle expérimental est 
facilité si on se restreint à ces gaz. La viscosité et la conduction sont les effets classiques les plus 
importants, mais le rayonnement conduit à des équations de dispersion et d’absorption identiques 
à celles que l’on trouve pour la relaxation. L’irréversibilité ne diffère donc pas essentiellement 
suivant son origine inter ou intramoléculaire. 

On peut äonc conclure que la théorie des effets d'absorption et de dispersion de translation est 


actuellement vérifiée : nous traiterons de la relaxation intramoléculaire dite aussi structurale 
dans un prochain article. 


Abstract. — The study of dispersion and absorption effects also called ‘‘ classical effects ” 
or effects of translation since they affect these molecular degrees of freedom is most useful toseparate 
the global effects which are the only ones accessibles to experiment, from those which are due to 
intramolecular relaxation. As this latter cannot affect the behaviour of monoatomic gases, the 
experimental control is easier if it is confined to these gases. Viscosity and conduction are 
the most important classical effects but radiation leads to equations of dispersion and absorption 
identical to those which are found by relaxation. Thus irreversibility does not essentially differ 
according to its inter or intramolecular origin. Thus it can be concluded that the theory of effects 
of translational absorption and dispersion is now verified : in a forthcoming article intramolecular 


relaxation, also called structural relaxation will be treated. 


Trois équations définissent les mouvements d’un 
gaz si on se limite à la propagation unidirec- 
tionnelle, ce que nous ferons dans tout cet article. 

Si l’on désigne par l'indice e les petites variations 
des variables, pression p, masse spécifique & et 
température 0, autour de leur position d'équilibre 
Po Po do, l'équation du mouvement et l’équation 
de continuité sont respectivement, u — dx/dt étant 
la vitesse selon l’axe Ox : 


po(dufdt) = — (dpe[dx) en négligeant : u(du/dx) 
Pol(dufdx) = — (d2e di) u(de dx). 


À ces deux équations il faut ajouter la classique 
relation d’état. 

Les équations (1,1) et (1,2) sont valables, en 
première approximation, aussi bien pour les solides 
que pour les gaz, le fait que, dans les fluides 
gazeux, les ondes soient longitudinales, alors 
qu’elles sont à la fois transversales et longitu- 
dinales dans les solides, ne changeant rien à leur 
établissement ni à leur aspect. La troisième en 
revanche introduit la thermodynamique par une 
nouvelle variable, la température. C’est par elle 
que les distinctions entre gaz, d’une part, solides 
et liquides, de l’autre, commencent à s’accuser. 

_ Éliminons, par exemple p. et & entre (1,1) et 
(1,2) et l’équation de définition : 


(1,1) 


en négligeant : (12) 


2e 


_ À (1,3) 
4 Po Pe 


Il vient : 
due? Œ (I 4) 
Po dr? + NE > + 


équation de propagation typique des ondes sonores 
planes, de célérité : 


US = 1/00 % 


x étant un coefficient de compressibilité. 

pe est fonction de l’évolution envisagée. 
Comment définir cette évolution ? Pour les gaz, 
dans le domaine des évolutions possibles, les limites 
que constituent l’évolution isotherme et l’évo- 
lution isentropique sont particulièrement com- 
modes à manier : aussi a-t-on tout de suite essayé 
de supposer que la propagation des ondes sonores 
était, soit isotherme, soit isentropique. 

La première hypothèse est due à Laplace, qui 


pose : 
1 
D (£:) : 
Po Pe 0 


US — 1/p0 Lo: 


D'où : 
(1,5) 


La discordance avec les résultats expérimentaux, 
l'absence de dispersion de la vitesse du son dans 
tout le domaine sonore malgré ce que lon savait 
sur la lenteur des échanges thermiques ont fait très 
rapidement abandonner la propagation isotherme 
pour l’autre hypothèse extrême : la propagation 
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isentropique, à laquelle on s’est tenu. En ce cas : 


US = 1/00 X8 = (pe lee (1,6) 


lPindice zéro pour U, rappelant que cette for- 
mule (1,6) n’est valable que si les conditions d’isen- 
tropie introduites par la définition même de y, sont 
rigoureusement respectées expérimentalement. 
Nous utiliserons par la suite sans indice 0 les cha- 
leurs spécifiques c et C'et leur rapport y = CJc. 
Ces paramètres seront toujours relatifs à un gaz en 
équilibre intramoléculaire et devront avoir été 
déterminés, conformément à leur définition, lors 


d'une évolution réversible. Les quelques discor-. 


dances résiduelles, inévitables puisqu'il s’agit d’une 
limite idéale, restent calculables en tant que 
simples corrections. 

Examinons successivement les causes de per- 
turbation du schéma exprimé par (1,6) qui uti- 
lise (1,1) et (1,2) et l’isentropique. 

Les unes viennent d’une toute autre origine que 
des hypothèses thermodynamiques sur la nature 
des compressions et détentes ou sur la valeur du 
coefficient de compressibilité. La considération de 
l'équilibre de l’élément de gaz dr exprimé par (1,1) 
néglige par exemple les forces sur la surface laté- 
rale ; s’il y en à, dues à la viscosité du milieu, 
c’est au niveau de la mécanique rationnelle qu’il 
faut calculer les corrections, et on peut, comme 
Stokes, rechercher directement les termes correctifs 
sur (1,1). 

Le manque d’isentropisme peut aussi être d’ordre 
thermodynamique. Une particularité affecte cepen- 
dant la propagation des ondes ultrasonores ou la pro- 
duction d’ondes stationnaires. Le gaz n’est plus en 
équilibre de température en toutes ses parties : il 
n’est plus question de confiner l’isentropisme à sa 
définition habituelle d'absence d’échange calori- 
fique entre le gaz et l’extérieur ; c’est pour chaque 
tranche élémentaire qu’il doit être réalisé. Les 
perturbations tiennent à deux causes principales : 
la conduction de la chaleur et le rayonnement ; 
nous conviendrons de négliger les phénomènes de 
convection. 

Les échanges de chaleur par conduction ou 
par rayonnement, les frottements visqueux ont 
cependant tous un point commun :ils soustraient de 
l’énergie mécanique à l’onde sonore et la dissipent, 
introduisant lirréversibilité dans la propagation, 
L'hypothèse des compressions et détentes isen- 
tropiques doit donc subir une double retouche. Les 
évolutions ne sont plus rigoureusement isentro- 
piques, non seulement par manque d’adiabaticité, 
mais aussi par manque de réversibilité. Nous cal- 
culerons successivement l'effet de ces trois causes 
sur la propagation sonore ou ultrasonore. 


[,1. La viscosité : Théorie de Stokes. — Stokes 
modifie l’équation (1,1) du mouvement en intro- 
duisant les forces exercées par la viscosité sur la 


No1 


surface latérale de l’élément, ce qui conduit à : 
olduJdt) + (dpeldx) — Cd? u[dx?) = 0. 


Admettons toujours le milieu justiciable des 
équations (1,2) et (1,3) pour un phénomène 1sen- 
tropique, donc : 


(1,7) 


RAS 


et soumettons-le à des ébranlements longitudinaux 
harmoniques de pulsation © ; les solutions pour 7e, 
ee, u seront de la forme : 

pe = Pe exp [i(ot — kx)] 


ee — Re exp [j(ot — kx)] (1,8) 
uw — UÙ exp [j(ot — kx)] 
k étant lui-même imaginaire : 
Let ee (1,9) 
L'équation de Stokes devient alors : 
60 OU — Pe k + j € Uk? = . (1,10) 
Les équations (1,2) et ([,3) conduisent à : 
Po Uk = ORe 
PelRe = 1/[(po Xs). (1,11) 
D'où l’on déduit d’après (1,10) : 
Pie: "20 
Re HE J ; Po (1,12) 
et : 
o 1 pAl ; 
ke = 5, + it] 
En posant : 
1 
qe 
il vient 
ARR (1,14) 


©? + w? 


que nous pouvons écrire sous une forme plus 
générale : 


RER US. (1,14 bis) 


Toutes les perturbations apportées à la propa- 
gation du son que nous étudierons ultérieurement 
conduisent à une équation du type (1,14 bis). 

D’après (1,8) et (1,9) la partie réelle k, est liée à la 
vitesse de propagation U par la relation : 

7? 

mi 
« représente le coefficient d'absorption. On obtient 
ces deux quantités en identifiant (1,14 bis) à : 


Me (1,15) 


KR = kB — 42 — 2jk à. 


D'où 
ie — FAI 
[1 + (1 + Si) 
L à (1,16) 


No1 
Soit, dans l'hypothèse où SR € 1 : 
Ur = of 
DIE SET UE (1,16 bis) 
en à 
R2 


L'application des formules (1,16) donne dans 
le cas particulier de l'équation (1,14) : 

2\1 = li 
Po 12 (1 € SE 1 Je 
= EE ra) —— 


2 
2C0y 1 ea ee 
CO 
o\L 1 
: Les 
7 1 o ( 2 
"vit En PNEU 
OY 
que — 2e @? + wo? 
Po @v 1 + [1 + (w?/o,)]1/? 
à 
= aug + er : (1,17) 


@y  1+ [1 + (w?/o$)]#/? 


en reprenant l’expression (1,6) de la vitesse U,. 

On aboutit donc à une véritable dispersion et à 
une absorption que nous envisagerons pour deux 
cas particuliers. 


a) © € y. | 
C’est le cas le plus fréquent pour les gaz. On 
trouve en première approximation : 


wo? € | | = 1 a (152) 
Ce RTE Bail 24, 8 © 
OS 
AU = U — — |. 
2. 8 < 


En particulier, si © —>0, « — 0 et U devient 
égale à U,, valeur de la vitesse pour une propa- 
gation rigoureusement isentropique. 


b)'o> 0, | 
Cette autre hypothèse simplificatrice conduit à : 


1e SE 
leule 
La vérification de ces formules implique la con- 


naissance du paramètre € non directement lié aux 
données expérimentales. 


(1,19) 


1 
De A =, FT 
Ov 


LE PARAMÈTRE €. — D’après les équations de 
Lamé, appliquées à un mouvement unidirectionnel, 
l'effort visqueux n, est donné par : 

il 
Nx — UE + 2U (ex ma c.) (1,20) 
u et x’ étant les coefficients correspondants à ceux 
de Lamé, €, la vitesse de dilatation cubique, #, la 
vitesse de dilatation linéaire selon l’axe Ox. Et 
l’on a : 


ÉChrMEx elu Ce 
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L’équation de Stokes devient alors : 


_ 1 pe le 
Pe = + (u Pate (1,21) 


X8 Po 
uw" + 2 est le coefficient de viscosité d’ensemble, 
u le coefficient de viscosité classique habituellement 
mesuré. Stokes pose en outre, par analogie avec 
l’état d'équilibre statique, que la vitesse de dila- 
tation cubique ne donne lieu à aucun effort vis- 
queux : 


TETE (1,22) 


En tenant compte de (1,7), de (1,21) et de (1,22), 
il vient : 


ET (1,23) 
HR Rare 
EE Xs Po 3 Po 


C se déduit alors immédiatement de la mesure de p. 


La vérification des expressions (1,17) ne peut se 
faire puisqu'il est impossible de séparer dans les 
résultats expérimentaux la part due à la viscosité 
de celle provenant d’autres effets dont le manque 
d’adiabatisme que nous examinerons maintenant. 


1,2. Les pertes de chaleur par conduction : équa- 
tions de Kirchhoff. — Chaque élément de volume 
soumis à des compressions et détentes perd en fait 
une quantité de chaleur dQ telle que, par unité de 
masse : 


dO:= cdô + Tdn : = 0, «1x0 (1,24) 


«, coefficient de dilatation : 


1 [dv 
= 5), 
en supposant applicable ïei cette relation valable 


uniquement pour des phénomènes réversibles. 
D’après l'équation classique de la conduction : 


do _ E2°0 
JMS 


(1,25) 


Ë pouvoir conducteur spécifique du gaz considéré. 
L’élimination de Q et de 0 entre ces deux équa- 
tions (1,24) et (1,25) et la relation : 


d6 — (20/dp) dp + (d0/dv), dv 


sOIt : 
X0 al Xe 1 
d8 = — d — do — dp — d 1,26 
er AO os (1,26) 
conduit à : 
EpUe 1 de] D 7 (or) 
00 € dx? 20 Xe dx? dt Po Xs dt + 


Supposons toujours ce milieu soumis à des vibra- 
tions planes et harmoniques d'équations (LB). 
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L'’équation (1,27) devient : 
‘RE 16 E k? 
Ps] He | Re | nes | à 
HS PS SE PoXs  PêXe © 
D'après (1,1) et (1,2) : PB JRe = «w?Jk?. 
Done : 


k° [: | 
— (©) 
Po 
k? est donc donné par une équation du 4 degré. 


Dans les parenthèses de (1,28), remplaçons Æ? par 
w?JU? et posons : 


At. a k°? 
TOME | Co) 1,28 
io tee (1,28) 


o € 


(1,29) 
Il vient : 
Y ©? + où — j wuwc(y — 1) 
©? + owë 
Les perturbations restant faibles, U peut le plus 
souvent être remplacé par U, dans l’expres- 
sion (1,29), « est alors indépendant de © et, en 
utilisant l’approximation (1, 16 bis) : 
4 w? 


(1,30) 


k2 — ©? po Xs 


es 9 U?2 ©? + où oelY ‘ 
4 U ; ©? Oe X9 — Xs 
EN PIC REE 
, O2? + wo 
EN I TU RUES |; 
0 2 a GO (1,31) 
Les équations (1,31) diffèrent des  précé- 


dentes (1,17) par le fait que «x (1,29) est en toute 
rigueur fonction de U donc de la fréquence 
imposée ©. 

D’après (1,31), lorsque © — 0, soit © < «: 


= 10 À Y + 41\ w? 
eg UE oe W 1) l ; Fe ] | 


œ —> OÙ 
; Y — 1\ w° 
et ua f1 —( 5 ) ©] 
U =, 
Lorsque © — co soit © > ox : 
Le Le Fe YŸ 4 (1 Y +1 es) 
RAADRE Vy 2 UE 


ax tend vers une limite finie et 


TEE US [1 i Y —1 ee] 
Vy Pro 


No1 


effets visqueux et conductifs, il suffit d'admettre 
comme valables les équations (1,2) et ([,8) et 
d'adopter simultanément les équations (1,7) et 
(1,27). k est alors donné par une équation du 4 degré 
que l’on doit résoudre par calcul numérique (°) : 


RUE [! Co | 
Po € LY Po Uÿ 


— k? [1 


Li (e- =] D EP 1 _ 
Î Jar Po E 20 c : Uz ( (1,32) 
et que l’on peut mettre, dans le cas de gaz mono- 
atomique où : 


Y=5/38 EE (5/2)ue 
sous la forme : 


k4 US$ Ê Il sue ci cou }] 
Oo? 2 US Po 3 VUS po 


k2U 23.{ w 
ENT Ê Ne AE] | 

L'ensemble des effets de la viscosité et de la 
conduction sur l’absorption et la vitesse est alors 
fixé par les équations (1,32) ou (1,32 bis) de Stokes- 
Kirchhoff. 

Les équations du mouvement utilisées jusqu'ici 
s'appliquent au gaz considéré comme un milieu 
continu. Certains auteurs [1] ont montré qu’en 
partant de la théorie cinétique, dans le cas d’un 
gaz monoatomique visqueux et conducteur, les 
équations (1,7) et (1,25) auxquelles reste adjointe 
l'équation (1,2) de continuité devaient être rem- 
placées par : 


1 — 0. (1,32 bis) 


du dPe - DUT? : Din 
Po dt dx D 4 3 3e NT: 
2) 2 
+2, mo à 
Fa Re 0 (1,33) 
et 
dQ LE 0 DO EU eo te 
po 3 1 08 2 a ©2 > ne A SA 
: FO 3 po dx 3 me 


où €, €2, Egets — p /u dépendent de la loi de force 
intermoléculaire adoptée. Le tableau I donne les 


{:) Lorsqu'on se limite à des fréquences relativement 
basses, on peut remplacer dans (1,32) k4 par otJU4 et 
négliger le terme en EC. Adoptant la simplification (1,16 bis) 
on parvient à la relation très couramment employée: 


U — (1 /00%0)!/?, valeur dela vitesse des ondes pour U = Uy 
une propagation isotherme. DUT £ 
Si l’on veut tenir compte simultanément des "202 Ë Mdr Aie |. (1,32 ter) 
TABLEAU I 
MAxWELL LENNARD-JONES , SPHÈRES 
ELASTIQUES 
e, — 45 J8 (25/8) X 1,008 — (25/8) X 1,017 
Co 3 3 X 0,855 3 X 0,800 
c; 2 2 x 1,007 2 X 1,044 
SU] 1,000 1,004 1,009 
2,81 3,15 3,24 


(voir formule (1,43)) 


No 


valeurs de ces paramètres dans le cas de forces 
r, distance séparant deux molé- 
cules) et en 1 /r° — 1 fr? proposées respectivement 
par Maxwell et par Lennard-Jones et dans le cas 


répulsives en 1 /r5 ( 
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TABLEAU Il 


du schéma assimilant les molécul 
élastiques. 


équations 
pour €1, £,, €; et $s les valeurs de Maxwell et 


11 
es à des sphères 


(1,33) et(1,34) en adoptant 


Ne 2, 
res y c) 


2 A? 
Viscosité CA 1 
1 1 fa ES VESTE 1 2 
’ EC) (oi au? (4 rs (Gi) 4l 
æ#' Va À «2 Ok &w?\1 
re | 142 "i+fi+) 
LS cx oÿ 
HO INT 5e 3 w? 
LE L HA Ne TRE AL2 -  — 
Res DATE on Ê 3 | où [1 uv | 
. 1 Fooyli ol À. 
oO > Ov ‘qu, | 9 IE de 
Conduction HU «? RO de 
@c(Y 1) (2 as = 
ape ec Hire + a O2 Y + où 
es es Le Dee &| “UMA {y = 1) 
[Q) << Oc DU ra Ï ) | a 9 «8 0 | ot 2 
1 rl | A2 Y — 1 où 
2 Wen = 1 [1 ER EU 2< [1 ES | 
© > @e JU, VS e(Y ) y © _ Me 
Rayonnement | 1 œo 1) que 2° ce 
DUO EN OE ; © + of Y 
Or — . 1 1 © ! 4 Van | US [1 | Mines l | 
F5 a ë V+ Or Du | mn 2rro: 1 Ro: 
l Y + 1.0} uw [1 y—1) À] 
OZ Or ‘AL, œrlY CE 1) [1 SE Ro = ui 1 (Y ) ©? 
TABLEAU IT (suite) 
VALEURS be æ@, 2aù ET Ù POUR © —> 0 ET © > co 
TER = w OO —> CO 
P£RTURBATIONS UE s 
ENVISAGÉES u | 2010) al * 2ax (*) al 
SR AMC ÉRRR = limite 00 
Viscosité 0 0 MALE finie) 
1e) 
1 \1 TE (2): ; (Ce 
Conduction 0 0 He ee AE DE (BE Xo/” | 20 Xe 
1 \1 1.q (Fe) 0 ee à Je 
Rayonnement | 0 0 (re AE LAS Po Xs 
(*) 24X représente la variation relative d’énergie entre deux tranches d’abscisses x et x + 2 


A 


avec y — 9/3 conduisent à : 


k6 US EI ou | L 10 ( oz }] 
©$ 127 VUS po RATER 

a ma ( ou | - 3 

APRES 2 


TETE 


(ee ] 
Dates 


242 23 > | 
UT 
0ro 


» 
(0) 


(1,35) 


et à un plus haut degré d’approximation à : 


kS US ÊE X 2 ( OL | 15 { ou }] 
æ | o72. us QUES 
È 8 L5e (Wobs 


67 { ou }] 
o 24 ] U$ Po 
L RANUA F 275 ( oz | | ot f Le ( ot ] 
DCS ME) tee 
k2U2 1 = 29 { ou ] , 
ARE | TG l U2 20 Ï 


Les équations (1,35) et (1,36) doivent être uti- 
lisées numériquement, les formules approchées que 
l’on peut en déduire pour Æ ne pouvant convenir 
quelles que soient les valeurs de l’expres- 
sion ouf UE p,. 

Les calculs précédents, tenant compte des per- 
turbations dues à la viscosité et à la conduction, 
montrent que le fait d'augmenter la fréquence 
aggrave l’irréversibilité, (voir tableau IT). Il reste 
à examiner maintenant le défaut d’isentropisme 
dû au rayonnement. 


. xs US$ FE ( ot je 


4 =:0. 


(1,36) 


1,3. Le rayonnement. —- Un élément de 
volume dr à la température 0 perd vers les milieux 
à la température 0, la quantité de chaleur rapportée 
à l’unité de masse : 


dQIdt = — qe(9 — 0) (1,37) 
en utilisant, issue de relations classiques, l’expres- 


sion de dQ applicable aux transformations réver- 
sibles : 


C 
Po Xp 


et l'expression (1,26) de d0, et, en dérivant (1,37) 
par rapport au temps, il vient : 


Je 
dO = — cd? de (1,38) 
Xp 


X6 d?p CARO EEE AO 
Gp D? Poup dE? MRC 
soit encore : 
tbe TO MOMIE TARDE 
dt? dt Xs po D?  XBPodt 


Intégrons en fonction du temps : 


. 2 ve 
Le — % Pe + Le — 9% pe = 0. (1,39) 


0 PoX 


Si le milieu est, là encore, soumis à des ondes 
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No1 


. 
. 


planes et harmoniques où p, et & ont la forme (1,8) 


Re (a — 4) — pit — 9%) = 0. (140 
Po nt 
En tenant compte de (1,1) et (1,2) 
Pe]Re Eu w?/k?. 
D'où la relation entre £ et © : 
, . . 4 
©? po (us = 1e x) = (o— 4). 
Et, en posant &, — q{Y, : 
2. 0% + © 6 — j @or(Xg — s) ; 
RERO" £ È SE of J 2 (1,41) 
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2 
uw = 4, [1 = L #| Re 
. Le fait d'augmenter la fréquence permet dans ce 
cas de retrouver, pour la vitesse des ondes, la 
valeur U,. 

Nous sommes ainsi devant trois formes d’irréver- 
sibilité dans la propagation ultrasonore plane ; 
toutes, malgré un mécanisme différent, aboutissent 
à une dispersion de la vitesse et naturellement à une 
absorption. w/p, apparaît comme un paramètre de 
similitude pour la viscosité et la conduction (2). 


(?) Les expressions telles que (1,32 bis) et (1,35), (1,36) 
montrent que k? U5]w? depend uniquement de l’expres- 
sion œufU$ 65. Si l’on admet U? et indépendants de la 
pression, ce qui est parfaitement justifié, même pour un 
gaz loin de l’état parfait, dans un large domaine de varia- 
tion de po, &/ps — wfps apparaît comme un paramètre 
de similitude. 

Par contre, de l’examen de la formule (1, 41) relative 
aux perturbations dues au rayonnement, on déduit que 
k?U5 Jw? est fonction uniquement de wo — wy{q. Or g 
n’est pas proportionnel à p4, mais au contraire sensiblement 
inversement proportionnel. 


0 duo. sfr d'a ftè sai 
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confirmant ainsi l'intuition suivant laquelle il 
importe seulement que la fréquence varie par 
rapport à la durée du libre parcours moyen, tandis 


Courbes théoriques : 


Équation (1,82 bis) 
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F1G. 4. — Argon 0 — 300 C. 


dres de grandeur relatifs aux effets que nous 
venons d'examiner ; nous prendrons l'air à titre 
d’exer nle pour donner un sens à la diffusion. 


TABLEAU III 


Fr NE 
o? 

POUR LES DIFFÉRENTS TYPES D'ABSORPTION CLASSIQUE 
DANS L'AIR A 20 0C, PRESSION 1 ATMOSPHÈRE 


œ 
VALEURS DE — — 
V 


œ 
CAUSES D’ABSORPTION SE MUPICMmEs sec 


Viscosité 99 
Conduction 38 
Diffusion 7,4 
Rayonnement 0,42 


(à 6 000 cycles/sec) 


Les ordres de grandeur relatifs sont à peu près 
les mêmes pour tous les gaz, sauf l'hydrogène et 
l’'hélium, où les irréversibilités par conduction sont 
nettement plus élevées. 
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que sa nature électromagnétique place le rayon- 

nement à part. 

Dans les mélanges gazeux, on devrait aussi tenir 
compte de la diffusion [2] dont l’effet quanti- 
tatif, bien qu’encore négligeable, est très sou- 
vent supérieur à celui du rayonnement. C’est 

y le type d’équation auquel on aboutit pour ce 

dernier qui nous l’a fait traiter, pour des com- 
paraisons ultérieures. 

Avant de passer à la comparaison avec les 
résultats expérimentaux, fixons quelques or- 


083 Courbes théoriques : 
0,30 Équation (1,82 bis) 
= Sn Se es (LES) à 
ee ” (1,36) 
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0,6 GE 
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Fic. 2. — Argon 6 — 300 C. 


1,4. Comportement expérimental des mono- 
atomiques. — Les figures 1 et 2 donnent les points 
expérimentaux de Greenspann [3] relatifs au 
rapport (U—U,)/U, et à l'expression «U,/o en 
fonction de w/p, pour l’argon et les courbes théo- 
riques calculées à partir des équations (1,32 bis), 
(1,35) et (1,36). Ont été représentés également sur la 
figure 1 les points expérimentaux indiqués par 
Boyer [4] et les valeurs théoriques calculées par cet 
auteur à partir d’une formule approchée tirée des 
équations (1,33) et (1,34) et limitée aux termes du 
1eT ordre : 


2 2 
at a, [1 + B.105 =] (1,43) 
0 


Pi #T 
où p, est exprimé en cm de mercure, et où B est 
une constante dépendant des lois de force inter- 
moléculaires choisies (voir tableau T). 

Les figures 3 et 4 sont relatives à l’hélium. 

Les courbes théoriques traduisant la for- 
mule (1,43) s’écartent de l’expérience comme des 
autres théories lorsque w/p, croît, en particulier 
lorsque p, décroît pour « constant, processus 


1% JOURNAL DE PHYSIQUE 


habituellement utilisé par les expérimentateurs. Or 
on devrait obtenir la convergence de toutes théories 
basées sur (1,33) et (1,34) lorsque p décroît, c’est-à- 
dire lorsque le gaz tend vers l’état parfait. I faut en 
effet prendre garde, quand on modifie la conception 
du milieu continu pour introduire la théorie ciné- 
tique, aux approximations faites ensuite dans le 
calcul. A ce titre, la tentative de Boyer semble assez 
décevante, en ce sens qué des développements 
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limités, type (1,43), ne sont valables que pour les 
valeurs du paramètre w/p, où ni U, ni « ne pré- 
sente de particularité intéressante. Les contra- 
dictions auxquelles on aboutit, lorsqu'on adopte 
de tels développements pour les domaines où lab- 
sorption et la dispersion deviennent notables, sont 


NO 


illustrées par Panomalie que nous signalons Gi- 
dessus. PT 

Une remarque semble nécessaire sur le sens des 
évolutions qu'impliquent ces effets visqueux et 
conductifs : ils augmentent la célérité U et l’absorp- 
tion «, infirmant ainsi le préjugé habituel suivant 
lequel lisentropisme est favorisé par une augmen- 
tation de la fréquence. 

Lorsqu'on utilise une masse de gaz constante et 
homogène soumise à des évolutions, harmoniques 
pour fixer les idées, l'augmentation de fréquence 
favorise bien l’isentropisme. Mais, si on transpose 
ce résultat aux évolutions créées par des ondes 
ultrasonores stationnaires en omettant de consi- 
dérer une masse de gaz proportionnelle à À pour y 
calculer les défauts d’adiabatisme, on assimile à 
tort l’ensemble du gaz soumis à la propagation à 
celui contenu dans un cylindre fermé par un piston 
mobile. Augmenter la fréquence pourrait sembler 
cependant favorable à l’adiabatisme d’une manière 
indirecte, l'amplitude imposée au gaz par les oscil- 
lateurs à quartz ou à magnéto-striction décroissant 
en général quand la fréquence croît. Mais le coef- 
ficient d'absorption «, donc le rapport des pertes 
par conduction à l’énergie mise en Jeu n’est pas 
fonction de l'amplitude, et l’adiabatisme relatif est 
indépendant de cette dernière. Pratiquement, dans 
tout le domaine où on forme des ondes sonores 
planes, il n’y a aucun intérêt à augmenter les fré- 
quences, au contraire. 

La considération des résultats précédents montre 
que, pour les monoatomiques, les équations de 
Stokes-Kirchhoff groupant le terme principal U, 
avec les corrections dues à la viscosité et à la con- 
duction sont en accord avec l’expérience dans un 
v/ps étendu jusque vers 200 mégahertz/atmo- 
sphère pour l’argon, soit 1 200 C. G. S. pour w/p,, 
Bien que de manière plus imprécise, on arrive à la 
même conclusion pour l’héllum et on peut done 
conclure, pour les monoatomiques, à une concor- 
dance honorable entre les théories de Stokes- 
Kirchhoff et l'expérience. C’est seulement quand la 
longueur d’onde imposée se rapproche du libre 
parcours moyen qu’un écart commence à se mani- 
fester, ce qui est tout naturel puisque le concept de 
milieu continu ne convient alors plus du tout à la 
nature du gaz [5]. 

Dans ce cas, la considération des équations de 
Burnett permet d’aller plus loin ; le calcul numé- 
rique complet développé par Greenspann [3]illustre 
une concordance convenable pour «U, / w et surtout 
pour (U—U,;) JU, jusque vers 3 000 mégahertz/- 
atmosphère pour v/ps (soit 20 000 C. G. S. 
pour ©/p5). L'introduction des termes d'ordre supé- 
rieur à AS dans ces équations est cependant déce- 
vante en ce sens que, si elle améliore bien la concor- 
dance avec l’expérience pour (U — 4) IU,, elle 
la perturbe au contraire pour «U,/o. Malgré cet 
obstacle on peut conclure des travaux entrepris 
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sur les monoatomiques que les effets de disper- 
sion et d'absorption sont remarquablement confir- 
més par les théories développées dans le cadre des 
effets visqueux et conductifs. 
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Dans un prochain article, nous traiterons du cas 
des polyatomiques, qui, après une perturbation 
rapide, présentent des déséquilibres entre leurs 
différents degrés de liberté. 


Manuscrit reçu le 30 juin 1958. 
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INFORMATION 


Réunion de Spectroscopie moléculaire. La prochaine 
réunion du Groupe européen de Spectroscopie moléculaire 
se tiendra à Bologne, Italie, du 7 au 12 septembre 1959. 
Pour tous renseignements, s'adresser au PT A. Mangini, 
ou Chimica Industriale, Viale Risorgimento 4, Bologna, 
Italie. 

Les principaux sujets traités seront : 1. Niveaux d’énergie 


relatifs à la dispersion rotatoire moléculaire serontTéga- 
lement étudiés. 

Les adhésions sont reçues jusqu’au 15 février 1959. Les 
résumés des communications (250 mots au maximum) 
doivent être présentés avant le 15 mars, et les communi- 
cations elles-mêmes avant le 30 avril. : 

Le comité scientifique du congrès est composé de 


et spectres ultraviolets. — 2. Spectroscopie infrarouge. — membres du Comité de Spectroscopie moléculaire de 
3. Spectroscopie Raman. — 4. Spectroscopie des micro- VI. U. P. A. C. 
ondes et résonance magnétique nucléaire. — Les problèmes 

ERRATUM 


Article de MM. Benoît, Grivet, Guibé et Sauzade, J. Phy- 


sique Rad., 1958, 19, 905. 


La photographie n° 3 a été placée par erreur figure 6, 


D'RJ07E 


La photographie n° 4 a été de même, placée figure 3, 
p. 906, et la photographie n° 6, placée figure 4, p. 906. 
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MESURE DE LA DIFFUSION MULTIPLE DES MÉSONS x 


S’'ARRÊTANT DANS L'ÉMULSION NUCLÉAIRE PAR LA MÉTHODE DE LA TANGENTE 


Par Bricirre DÉPAUX, 


Faculté des Sciences de Paris. 
Laboratoire de Physique. Enseignement. 


Résumé. — Nous mesurons la diffusion multiple de 40 mésons + qui s’arrêtent dans une émul- 
sion nucléaire G5 au moyen de la méthode de la tangente. La direction d’une cellule double est 
déterminée par la direction moyenne des deux cellules simples consécutives qui la constituent. 
Nous pouvons ainsi éliminer le bruit de fond entre cellule double et cellule simple. Les formules 
proposées pour l'élimination sont vérifiées par cette étude. La constante de la diffusion multiple 
déduite de l’étude des mésons II dont l’énergie est comprise entre 6 et 12 MeV est en bon accord 
avec la valeur théorique, dans la limite des erreurs expérimentales. La distribution des angles de 
diffusion multiple est gaussienne pour les angles inférieurs à 4 fois l’angle moyen, mais pour les 
angles supérieurs, il existe un écart significatif. Le rapport entre la différence seconde et l’angle 
de diffusion multiple est du même ordre de grandeur que celui trouvé par la méthode de la flèche. 


Abstract. — We measure the multiple scattering of 40 x-mesons which come to rest in a 
nuclear emulsion G5 by means of the tangent method. The direction of a double cell can be 
determined from the mean direction of two consecutive single cells forming the double cell. We 
can thus eliminate the noise between double and single cells. Formulas proposed for the elimination 
are verified by this study. The multiple scattering constant deduced from the study of IT-mesons 
of energy 6 to 12 MeV is within experimental error, in good agreement with the theoretical value. 
The distribution of scattering angles is gaussian for angles 4 times smaller than the mean angle, 
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but for greater angles, there is a significant excess. The ratio between the second difference and 
the angle of multiple scattering is of the same order as that found by the sagitta method. 


Introduction. — La mesure de l’ionisation et de 
la diffusion multiple des traces dans l’émulsion 
nucléaire permet de déterminer la masse et la 
vitesse d’une particule. La diffusion multiple a été 
mesurée par de nombreux auteurs [1-3] au moyen 
de la méthode de la tangente préconisée d’abord 
par Goldschmidt-Clermont [4]. La méthode de la 
flèche proposée par Fowler [5] est maintenant très 
répandue dans les mesures de la diffusion multiple 
des particules de grande et moyenne énergie. Les 
mesures par cette dernière méthode peuvent être 
effectuées très rapidement, une fois la trace alignée 
parallèlement à la direction du déplacement de la 
platine. Le bruit de fond provenant de la platine est 
pratiquement négligeable en utilisant un micro- 
scope bien conçu. L’étalonnage de la constante de 
la diffusion multiple dans une région d’énergie 
moyenne à donné une valeur [6] ne différant pas 
de plus de 2 % de la valeur [7] théorique, mais les 
écarts des points expérimentaux sont assez grands. 

Pour les traces présentant beaucoup de grands 
angles de diffusion et pour les traces ayant une 

faible densité de grains, la méthode de la flèche 


n’est plus pratique. Car, dans le premier cas, il est 


nécessaire de changer souvent d’axe de référence, 
et dans le deuxième, il est difficile d’aligner la 
trace le long de l’axe et d’estimer avec précision la 
flèche des grains. Il est donc préférable de recourir 
à la méthode de la tangente, dans les mesures de la 
diffusion multiple des mésons x, des mésons p et 
des électrons de faible énergie. 


Tsai-Chü [8, 9] et l’auteur [9] ont proposé des 
formules d'élimination pour la méthode de la tan- 
gente. La direction d’une cellule double peut être 
calculée à partir des directions des deux cellules 
simples qui la constituent. On peut donc calculer 
l’angle de diffusion multiple d’une cellule double et 
éliminer le bruit de fond entre cellules doubles et cel- 
lules simples. Le bruit de fond pour la cellule double 
doit être plus petit que pour la cellule simple. Nous 
avons appliqué cette technique à la mesure de la 
diffusion multiple des mésons x en fin de parcours. 
Pour avoir un bruit de fond assez grand pour être 
mesurable afin de vérifier les formules proposées, 
nous nous sommes servis, dans cette expérience, 
d’un microscope ordinaire. La constante de la dif- 
fusion multiple obtenue pour des mésons x dont 
l’énergie varie de 6 à 12 MeV est en bon accord 
avec la valeur théorique ; elle est dans les limite 
des fluctuations statistiques. 


Conditions expérimentales. — Nous avons repéré 
dans une plaque Ilford G5 de 7,5 x 15 em? et 
d'épaisseur 625 microns, exposée en un paquet 
d’émulsions pelées au faisceau de protons de 6 BeV 
à Berkeley, 40 mésons x positifs ou négatifs qui 
s’y arrêtent. Le faisceau est très intense, et nous 
avons cherché des mésons + dans la région de la 
plaque où il doit y en avoir le plus. Nous les avons 
identifiés par la désintégration x+ —> ut — e+ ou 
par la formation d'étoile sigma. Nous avons con- 
servé pour notre étude de la diffusion multiple les 


cran vie On AE 


No 


traces de mésons dont le parcours résiduel dans la 
plaque considérée est au minimum de 3 000 mi- 
crons, et dont l’inclinaison par rapport au plan de 
émulsion est suffisamment faible pour ne pas 
avoir à introduire de correction de profondeur 
dans la détermination des parcours. De plus, pour 
ces traces, les corrections dues à la déformation de 
l’émulsion peuvent aussi être négligées dans les 
mesures de diffusion multiple. 

Pour effectuer les mesures de diffusion multiple, 
nous nous sommes servis d’un microscope Stiassnie 
muni d’un objectif à immersion 100 Koristka et 
d’une paire d’oculaires X 4. Nous avons mesuré la 
direction de chaque portion de la trace en l’alignant 
avec un fil de réticule fin placé dans l’un des 

oculaires dont le mouvement de rotation était 
repéré avec précision sur un goniomètre gradué 
au 1/10 de degré. 

Nous avons partagé chaque trace en cellules 
consécutives de 44,4 microns, et nous avons déter- 
miné la direction de chacune de ces cellules à 
partir de 888 microns de la fin du parcours. Nous 
avons mesuré 50 cellules par trace. 


Distribution angulaire de la diffusion multiple. 
— L’angle de la diffusion multiple est, par défi- 
nition, l’angle entre les directions de deux cellules 


300 


200 


100 
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FIG. 1. 
consécutives, c’est-à-dire la différence A6 entre les 
‘angles 0 des directions de ces cellules. Nous avons 
‘étudié la distribution des valeurs absolues des 


angles A6; pour cela, nous avons réparti les 
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mesures en 9 sections : dans chaque section corres- 
pondant à une bande d’énergie bien déterminée, 
nous avons calculé la moyenne A6 des A0 des 
40 mésons après élimination des grands angles, et 
nous avons normalisé chaque A0 par rapport à 
cette moyenne. Nous avons porté alors sur l’histo- 


gramme de la figure 1 le nombre de A6/AÛ compris 
dans un intervalle donné en fonction de la valeur 
de cet intervalle. Pour les petits angles cet histo- 
gramme s’accorde bien avec une courbe de Gauss 
rapportée à la même surface, mais pour ceux qui 
sont supérieurs à 4 fois la moyenne, 1l y a un écart 
très appréciable. 

À chaque valeur de A0 obtenue, nous avons fait 
correspondre A0?, et nous avons calculé le rap- 


port VAU2/A6. Si nous éliminons les grands angles 
supérieurs à 4 fois l’angle moyen, ce rapport est 
égal à 1,29 + 0,01 ; valeur en accord avec le rap- 
port des moyennes gaussiennes correspondantes 
VxJ2 = 1,25. Au contraire, si nous conservons les 
grands angles, ce rapport est très supérieur à la 
valeur précédente et égal à 1,37 Æ 0,03. Ceci 
montre l’écart significatif d’une distribution gaus- 
sienne pour les grands angles. 


Détermination de la constante de diffusion mul- 
tiple. — Dans la région d’énergie sur laquelle ont 
porté nos mesures, des calculs utilisant les courbes 
de Voyvodic et Pickup [7] nous ont montré que la 
constante de diffusion multiple ne variait pas de 
plus de 1 %. L'erreur que nous introduisons en 
supposant qu’elle conserve une valeur constante est 
par conséquent négligeable par rapport aux autres 
sources d'erreur. Nous pouvons admettre la rela- 
tion: 

| A0, = At X e(1fpu)? + «? 
dans laquelle A6, désigne l’angle de diffusion mul- 
tiple, À, la constante de diffusion multiple, #& la 
longueur de cellule exprimée en unités de 
100 microns, pv le produit de la quantité de mou- 
vement du méson par sa vitesse, et e le bruit de 
fond des mesures exprimé en degrés comme A6. 

Nous avons groupé nos mesures en 9 sections 
comprenant sensiblement le même nombre de cel- 
lules, et nous avons déterminé pour chaque section 
le (1/pv)? moyen et A0 pour les 40 mésons ; nous 
avons procédé à l’élimination des angles supérieurs 
à 4 fois la moyenne. Les mésons p provenant de la 
désintégration des mésons II ont une longueur 
moyenne de 605 Æ 9 microns, valeur pratiquement 
égale à celle de Barkas [10] et al. ; nous avons 
utilisé leurs résultats pour calculer les valeurs de pv 
des mésons IT. X, reste constant, mais pv variant 
pour ces mésons IT très rapidement, il en est de 
même de A6,, de sorte que la variation de AG est 
grande vis-à-vis de sa fluctuation statistique et 
qu’une détermination de Æ, aussi bien que de € 
est rendue possible par la méthode graphique. 
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La figure 2 donne la courbe obtenue en portant 
en abscisse les valeurs de {4(1/pv)° et en ordonnée 
celles de A6. Cette courbe est une droite dont la 


1,0 


0,444 
(pu? 
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rc? 


pente est ÆX? et l’ordonnée à l’origine e?. Nous 
avons obtenu de cette façon K, = 24,3 Æ 0,3 et 
e — 0935. La valeur de X, est en bon accord avec 
la valeur [9] prévue théoriquement qui est de 24,5. 


Méthode d’élimination. —- Nous avons cal- 
culé [8, 9] à partir des cellules simples la diffusiôn 
multiple des cellules doubles. Une cellule double 
formée de deux cellules simples consécutives dont 
les directions 6, et 6,,, a pour direction 


0, — (1/2) (4 + 6;;;) 
(fig. 3). 


FIG. 


Les angles A06,, différences entre 2 angles 6, 
consécutifs, définissent l’angle de diffusion multiple 
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pour les cellules doublés. Nous obtenons deux séries 
d’angles qui se chevauchent, en associant les cel- 
lules simples i et i + 1 oui —1eti. ère 

Nous pouvons caleuler la diffusion multiple par 
les formules suivantes : | 


A6 — KÈ? &(1fpv)? + €? pour les cellules simples. 
A0 — K3 2(1/pv)? + <°]4 pour les cellules doubles. 


La constante X, doit être légèrement plus 
grande que Æ,, mais comme X;, elle reste prati- 
quement constante dans le domaine d’énergie 
considéré. Nous avons déterminé le rapport K,/K; 
en utilisant la relation théorique : 


0,94 t Fa 
p? + 0,30 


donnée [71 par Voyvodic et Pickup, formule dans 
laquelle t est exprimé en microns. Le rap- 
port K,/X, calculé pour les 5 sections considérées 
conserve une valeur constante de 1,039. 

Il a fallu également déterminer le bruit de fond 
pour les cellules doubles en fonction du bruit de 
fond des cellules simples. 

Par définition, 


HE che} 1,25 + 0,80 (ioge 
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D'après les formules [11] rappelées par D’Espa- 
gnat, nous avons : 


D? — (diffusion pure)? + €? 


Dj Dis, — 1/4 (diffusion pure)? — 


D'où 
A6 = 2 (diffusion pure)? + e?]4. 

Nous avons groupé les mesures par sections 
telles que le parcours moyen de chaque section soit 
le même pour les mesures des cellules simples et des 
cellules doubles, et corresponde ainsi à la même. 
énergie. Pour accorder les cellules simples aux cel- 


N°1 


lules doubles, il a donc fallu faire légèrement 
chevaucher lès sections de mesures FE cellules 
simples. 

Les résultats des mesures et des calculs sont 


donnés dans le tableau I. Nous obtenons 
TABLEAU I 
A6, A6 (1 pv)? 
1278 + 0064 1,63 2 0,160 5.873 :10 
1,080040,054 4,17. + 0,12 4,214 .10—3 
0,992 + 0,050 0,984 + 0,10 319010 
0,914 + 0,046 0,835 + 0,084  2,503.10—5 
0,802 + 0,043 0,643 + 0,069 2,061 .10—® 
TABLEAU I (suite) 
po M, VA a 
2,608.10—5% 1,94 + 0,14 3,60 25,5 
1,874 .10 1,44 + 0,10 2,66 29,95 
1,399.10—? 1,47 + 0,10 2,76 26,7 
1,111.10—5 1,28 + 0,09 2,38 25,85 
0,915.10—2 0,99 + 0,07 1,81 21,65 


K, = 24,4 Æ 1,0ete — 0032 ; ces valeurs sont en 
accord avec les valeurs que nous avons obtenues 
précédemment par la méthode graphique. Pour un 
nombre total de cellules doubles indépendantes 
de 960, l’erreur sur X, due à la fluctuation statis- 
tique est égale à 0,7. 


Rapport de la différence seconde et de l’angle de 
la diffusion multiple. 
calculé la différence seconde A? 6 des 6. La diffé- 
rence seconde est particulièrement utile pour éli- 
miner l'influence des déformations de l’émulsion. 
Nous avons des relations d’une forme analogue à 
celles des différences premières avec des coefficients 
différents, c’est-à-dire 3/2 et Y : 


A20 = K? y2(3/2) (1 [po)? + €’? 
A°0 — Kè %(3/2) 210 Jpu)? + e?J4. 


MESURE DE LA DIFFUSION MULTIPLE DES MÉSONS x 19 


Le coefficient 3/2 résulte du calcul de la diffé- 
rence des A6, mais le coefficient y est une cons- 
tante empirique légèrement supérieure à 4. La 
moyenne pondérée de y est 1,05 + 0,03, valeur en 
accord avec celle donnée par Di Corato et al. [12], 
qui est de 1,06 + 0,01. Les résultats des calculs 
sont donnés dans le tableau IL. 


TABLEAU II 
nr, 2 2\ 1/2 

A26,A26, (4A262—A261)  w{1]pr pv) Re 
1,590 2,172 4,042 5,107.10—2 23,4 0,918 
1,420 1,884 3,490 9010 0255120506 
1,299 4,802 3,362 3,783.10—2 26,25 0,983 
1,209 1,664 3,101 3,370.10—? 27,2 1,052 
1,077 1,413 2,613 3,061.10—2 25,25 1,166 


Conclusion. — La constante de la diffusion mul- 
üple déduite de l’étude des mésons IT de faible 
énergie ne diffère pas de plus de 1 % de la valeur 
théorique. Il en résulte que la méthode de la tan- 
gente est d’une utilisation particulièrement recom- 
mandée dans le domaine des très faibles énergies de 
particules. 

Nous avons établi des formules d'élimination qui 
se sont trouvées vérifiées par cette expérience et 
par les calculs théoriques. Il est à remarquer éga- 
lement que le bruit de fond des cellules doubles est 
inférieur à celui des cellules simples. Il n’est pas 
excessivement grand, de sorte qu'avec un miero- 
scope plus précis, il est possible d'utiliser des cel- 
lules plus courtes. 

Je tiens à exprimer ici ma gratitude à M. le 
Professeur M. Morand pour l’accueil que je reçois 
dans son laboratoire. Je remercie également 
M. Tsai-Chü de l’aide qu’il m’a apportée dans ces 
recherches et M. C. Guerrier qui a veillé au bon 
fonctionnement du microscope. Ces recherches ont 
été rendues possibles par une allocution du Centre 
National de la Recherche Scientifique. 


Manuscrit recu le 16 octobre 1958. 
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MESURE DE LA RÉPARTITION ANGULAIRE DE LA DIFFUSION 
AVEC ÉCHANGE DE CHARGE PAR DES ATOMES D’ARGON 
D'’'IONS POSITIFS D’ARGON 40, DE PALLADIUM 108 
ET D’URANIUM 238 D'ÉNERGIE ÉGALE A 30 keV 


Par C. CASSIGNOL, F. DANIEL, G. DELMAS et G. RANC, 


Section des Réactions Nucléaires à Basse Énergie, G. E. N. Saclay. 


Résumé. — Description de la méthode employée, du dispositif expérimental, du mode opéra- 
toire. Les résultats sont présentés sur des diagrammes et critiqués. 


Abstract. — Description of the method used and of the experimental arrangement and pro- 


cedure. The results are presented in diagrams, 


Dans le but de contrôler un mécanisme impor- 
tant de contagion isotopique dans les séparateurs 
électromagnétiques à secteur, nous avons été 
amenés à étudier la diffusion à petit angle accom- 
pagnant l’échange de charge des ions de quelques 
dizaines de kilo-électron-volts d'énergie cinétique et 
de charge positive unitaire. Lorsqu'un tel ion ren- 
contre un atome gazeux, l'interaction peut revêtir 
différentes formes ; la plus fréquente dans ce 
domaine d’énergie est la capture électronique ou 
échange de charge. L’ion rapide est neutralisé et 
continue son chemin avec une faible variation 
d’énergie et l’atome-cible devient un ion positiflent. 

Il n'existe aucune théorie satisfaisante des colli- 
sions des particules lourdes aux énergies susdites. 
Il n’entre pas dans nos possibilités, ni dans nos 
intentions, d'entreprendre une étude de longue 
haleine sur ce sujet. 

Les sections efficaces correspondant à plusieurs 
types d'interactions ont été mesurées dans plusieurs 
cas particuliers [1]. 

Dans le cadre d’un vaste travail sur les collisions 
des particules lourdes, Fedorenko [2] a, pour la 
première fois, étudié les répartitions angulaires des 
particules rapides. Le domaine d’énergie qu'il à 
exploré déborde largement le nôtre et ses recherches 
sont aussi plus approfondies. Par contre, si nous ne 
nous sommes intéressés qu'aux répartitions angu- 
laires de l’échange de charge, nous avons dû étudier 
des ions provenant d’éléments inhabituels dans 
cette classe de recherche. C’est ainsi que, à côté de 
l’argon 40, nous avons expérimenté sur des ions de 
palladium 108 et d'uranium 238 à charge unitaire. 

Le principe de la méthode consiste à former un 
pinceau d'ions très délié et bien défini géomé- 
triquement et, au moyen d’une surpression locale 
d’un gaz donné, provoquer un certain nombre 
d'échanges de chaiges. Le résidu du faisceau pri- 
maire et les particules chargées provenant des 
autres types d'interaction sont déviés par un champ 


électrique transversal alors que les neutres rapides 
seuls continuent leur chemin dans une direction 
moyenne qui est la direction du pinceau avant la 
déviation. La distribution angulaire des neutres est 
relevée au moyen d’un collecteur mobile fonction- 
nant sur le principe de l’émission électronique 
secondaire. Aussi longtemps que l’on ne s’intéresse 
pas aux valeurs des sections efficaces, toutes les 
mesures sont des mesures relatives. Il n’est donc 
pas nécessaire de connaître la valeur exacte de la 
surpression ni le taux d’émission des électrons 
secondaires au collecteur. | 

Le dispositif expérimental comprend un collima- 
teur formant le faisceau délié primaire, une cham- 
bre de surpression, un condensateur déflecteur, un 
collecteur de mesure. L’infrastructure du dispositif 
est le Séparateur d’Isotopes de Saclay [3] (fig. 1). 


- 


COLLIMATEUR. — Il comprend : 

— un diaphragme carré À de 1 mm de côté (cons- 
titué en réalité par la partie commune de deux 
fentes croisées portées sur des plaques parallèles et 
voisines) (fig. 2) ; un collecteur de centrage B en 
bout de l’appareil. Le trou de définition a un dia- 
mètre de ? mm ; un diaphragme auxiliaire C placé 
avant le champ magnétique est destiné à compléter 
l’action de À. Il comporte un système de réglage de 
la largeur des fentes (de 1 à 30 mm). En réalité, ce 
dernier diaphragme se compose de quatre plaques 
en forme de demi-lune solidaires de deux paires de 
joints de Wilson (une paire horizontale et une verti- 
cale) se faisant face deux à deux. Un système de 
butées règle l’écartement minimum des bords recti- 
lignes des plaques (fig. 3). 

Diverses chicanes D sont disposées dans la 
chambre du séparateur dans le but d’interdire le 
passage aux 10ns parasites. Ces chicanes consti- 
tuant des résistances de pompage, l’enceinte de: 
mesure est évacuée par une pompe indépendante 
(Edwards 403). 


PAEAMFL! 


Fic. 2. — Diaphragme A et chicanes, 
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CHAMBRE DE SURPRESSION E. — C’est un 
cylindre métallique de 5 em de longueur et de 2 cm 
de diamètre que le pinceau d’ions traverse dans le 
sens de la longueur, en passant par deux trous 
de 5 mm de diamètre pratiqués au centre des bases. 
Cette chambre est solidaire d’un tube qui sert à 
amener le gaz ; ce tube coulisse dans un joint de 
Wilson porté par la paroi de l’enceinte en per- 
mettant de régler le centrage de la chambre et 
aussi de l’éclipser complètement. 

Le condensateur-déflecteur F est situé immédia- 
tement après la chambre de surpression. La tension 
appliquée est de 3 000 volts. Elle est plus que 
suffisante, compte tenu de l’écartement et de la 
longueur des électrodes, pour arrêter complètement 
toutes les composantes chargées du faisceau (fig. 4). 


COLLECTEUR DE MESURE G. — Il est constitué 
par une boîte, parallélépipède rectangle en feuille 
d’inox, dont la face d’entrée porte en son centre un 
trou circulaire de 1 mm de diamètre. Derrière ce 
trou, à l’intérieur de la boîte, se trouve une élec- 
trode de mesure. La boîte et l’électrode sont soli- 
daires d’un tube glissant dans un joint de Wilson 
centré dans une flasque latérale. La connexion de 
l’électrode de mesure est sortie dans l’axe du tube 
au moyen d’un passage étanche à haut isolement en 
plexiglas. La boîte est polarisée à quelques dizaines 
de volts au-dessus de la masse dans le but d’accé- 
lérer vers la boîte les électrons créés lors de l’impact 
des neutres rapides sur l’électrode de mesure et 
d'empêcher ceux qui sont créés sur les bords du 
trou de parvenir à l’électrode de mesure (fig. 5). 

Les translations du collecteur sont repérées parun 
index et une règle. Comme le trou de la face d’entrée 
n’est pas sur l’axe du tube et que ce dernier a la 
faculté de tourner dans le joint de Wilson, on doit 
faire les mesures à azunut constant, après avoir 


déterminé l’azimut qui permet d’atteindre le centre 


du phénomène. La constance de l’azimut est assurée 
de la façon la plus commode en contrôlant (au 
moyen d’un niveau à bulle) l’inclinaison d’une tige 
métallique solidaire de l’axe du tube et perpen- 
diculaire à celui-ci. 


Mode opératoire. — Le mode opératoire com- 
porte les étapes suivantes : ; 

1) Réalisation du pinceau délié primaire. 2) Éta- 
blissement de la surpression et élimination des 
composantes chargées. 3) Relevé de la distribution 
angulaire des neutres. 


RÉALISATION DU PINCEAU PRIMAIRE. — a) Le 
séparateur d’isotopes étant plein d’air, on enlève la 
flasque terminale dans laquelle le collecteur B est 
centré, et on la remplace par un réticule. On’aligne 
par visée optique le centre du réticule, le trou du 
diaphragme À, les centres des faces de la chambre 
de surpression, On règle les électrodes du conden- 


No 
sateur. Pendant cette opération et pendant la sui- 
vante, le collecteur de mesure est éclipsé et la boîte 
mise à la masse. 

b) On remet en place le collecteur B. On met 
l'appareil sous vide et on extrait un faisceau 
d’argon sous 30 kV. Le diaphragme C étant ouvert 
au maximum et les électrodes du condensateur F 


étant mises à la masse, on règle la valeur du champ 


magnétique jusqu’à obtention du courant maxi- 
mum au collecteur B. 

c) On rapproche par tâtonnements successifs 
jusqu’à la distance minimum les demi-lunes de C en 
veillant à conserver le maximum de courant sur B. 


ÉTABLISSEMENT DE LA SURPRESSION ET ÉLIMI- 
NATION DES COMPOSANTES CHARGÉES. — Àu moyen 
d’un robinet à fuites, on introduit dans la chambre 
de surpression de l’argon jusqu’à ce que le courant 
au collecteur B ait diminué d’environ 20 %. On 
met le condensateur sous tension, ce’qui a pour 
résultat d'annuler le courant sur B.Â 


RELEVÉ DE LA DISTRIBUTION ANGULAIRE. — Îlse 
fait après détermination de l’azimut du collecteur 
en mesurant le courant d’électrons secondaires du 
collecteur pour différentes positions de celui-ci, 
espacées de millimètre en millimètre. La distance 
du trou de G au centre de la chambre de surpres- 
sion est de 600 mm ce qui, compte tenu de la dimen- 


sion des diaphragmes, introduit un élargissement 


instrumental de quelques minutes d’angle. 

# Les lacunes expérimentales sont assez nom- 
breuses mais faciles à caractériser et à corriger. On 
pourrait envisager de commander par vis tous les 
réglages mécaniques, de réduire les dimensions des 
diaphragmes, d'introduire un troisième diaphragme 
de collimation analogue à C entre À et C, d’aug- 
menter la sensibilité des mesures, etc. 

La stabilité de la source a quelque peu laissé à 
désirer dans le cas de l’uranium. Plusieurs séries de 
mesures ont de ce fait été inutilisables. Une élec- 
trode de normalisation et l’intégration des mesures 
de courant auraient permis de surmonter cette 
difficulté. 


Résultats. — Les répartitions angulaires sont 
reportées sur les figures. Sur la figure 6 on notera la 
variation de l’élargissement en fonction de la masse 
du projectile, à énergie constante et la répartition 
fourchue de l’uranium. Sur la figure 7, on peut 
constater le rapprochement des deux maxima de 
l'uranium quand l’énergie s'élève. 

La quasi-concordance au voisinage du centre des 
distributions de l’argon et du palladium peut être 
attribuée à l’élargissement instrumental que nous 
avons estimé plus haut. Les distributions réelles 
sont done encore plus pointues et les maxima plus 
accentués que sur nos diagrammes. 

Dans.sa forme peu élaborée, notre étude nous 


Fre. 4. —— Chambre de surpression, condensateur-déflecteur et volet d’occultation. 
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H1G. 5. — Collecteur de mesures. 
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permet, par des recoupements, de caractériser la L'ensemble mécanique de mesure a été réalisé 
neutralisation des ions sur l'atmosphère résiduelle par MM. Bouyer et Knittner. 
d’un séparateur comme une cause importante de NS 
contagion isotopique [4]. Manuscrit reçu le 10 octobre 1958. 
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I. ÉTUDE DU RAYONNEMENT DE FREINAGE INTERNE, 
DE L’'AUTOIONISATION ET DES ÉLECTRONS ÉMIS SIMULTANÉMENT 
AVEC LE RAYVONNEMENT & DU PHOSPHORE 32 


Par F: SUZOR et G. CHARPARK, 
Laboratoire de Synthèse Atomique, CG. N. R. S. 


. Résumé. — On étudie, par analyse des impulsions en coïncidences de deux compteurs propor- 
tionnels, les spectres de photons et d'électrons de basse énergie (1 à 25 keV) émis en même temps que 


les rayons B, par ®?P. 


Le spectre de photons comprend la raie K de 2,3 keV provenant de l’autoionisation de l’atome 
résiduel d’intensité 5,5.10—4 par désintégration, et le spectre continu du rayonnement de freinage 
interne. Ce dernier concorde parfaitement avec les prévisions théoriques de 1 à 6 keV puis monte 


fortement au-dessus. 


Le spectre électronique comprend la raie Auger de 2 keV et un spectre continu important, 
d'intensité supérieure à celle prévue par la théorie de l’auto-ionisation. Le désaccord devient consi- 
dérable de 5 à 25 keV, où l’on trouve 6,6.10—% électrons associés par désintégration. 


. Abstract. — By analysis of the coincident impulses in two proportional counters study is made 
of the low energy photon and electron spectrum (1 keV to 25 keV) emitted,simultaneously with 


beta rays, by %?P. 


. The Photon spectrum comprises 2,3 keV K-line following autoionisation of the residual atom, of 
intensity 5,5. X 10—4 per disintegration, and the internal bremsstrahlung continuous spectrum ; 
the latter agrees perfectly with the theoretical prediction from 1keV to 6 keV and thenitrises strongly 


above ihe theoretical values. 


Theelectron spectrum comprises the 2 keV Auger line and an important continuousspectrum. The 
intensity of the latter is higher than predicted by autoionisation theory. The discrepancy becomes 
considerable from 5 keV to 25 keV where 6,6.10—% associated electrons are found per disintegration. 


I. Introduction. — Le spectrographe à double 
compteur proportionnel que nous avons cons- 
truit [1] a été utilisé pour l’étude des photons et 
des électrons de basse énergie émis lors de la désin- 
tégration du phosphore 32 simultanément avec le 
rayonnement 8. L’exposé de ces expériences mon- 
_trera qu’il est possible d’éliminer les phénomènes 
parasites, et que l’on obtient entre 1 et 12 keV le 
_ spectre des photons et entre 1 et 25 keV celui des 
électrons avec une bonne précision, et seulement 
un ordre de grandeur du phénomène pour des 
énergies supérieures. Les résultats seront comparés 
avec les théories du freinage interne et de l’auto- 
ionisation, faisant apparaître, surtout dans le cas 
des électrons, un désaccord important. 
_ La grande proximité des 2 compteurs est res- 
ponsable- de l'existence de phénomènes parasites 
qu’il est nécessaire d'éliminer ou tout au moins de 
réduire et d'évaluer avec précision. Les impulsions 
en coïncidence ainsi obtenues sont de deux sortes : 

a) un rayon Ê enregistré dans un compteur pro- 
jette dans le support de source, les écrans ou le gaz 
environnant un électron qui est enregistré dans 
Pautre compteur ; : 

b) un rayon 6 après un très court parcours de 
l’ordre de un ou plusieurs millimètres est rétro- 
diffusé par le gaz du compteur ou le bord de la 
rondelle sur laquelle est fixée la membrane sup- 
portant la source ou toute manière environnante et 
est ensuite enregistré dans l’autre compteur. 


Le premier de ces phénomènes parasites a été 
systématiquement étudié [2]par l’addition d'écrans 
d'épaisseur variable à des distances variables de la 
source. Les résultats en accord avec les prévisions 
théoriques montrent que l'effet peut être considé- 
rablement diminué et rendu souvent négligeable 
à la condition d'employer des sources radioactives 
obtenues par vaporisation sur un support très 
mince en matière plastique et de ne pas employer 
des écrans en contact immédiat de la source mais à 
des distances d’au moins un demi-millimètre envi- 
ron. 

Le deuxième.phénomène parasite dû à la rétro- 
diffusion du rayonnement 8 est beaucoup plus 
facilement contrôlable. On peut modifier la dispo- 
sition géométrique de la matière entourant la 
source ; c’est ainsi qu'il est toujours préférable de 
disposer la membrane porte-source sur un support 
métallique percé d’un trou de petit diamètre. Nous 
avons choisi, dans la plupart des cas, comme sup- 
port de source une membrane en matière plastique 
aluminisée d’une épaisseur de 40 microgrammes 
environ par centimètre carré, tendue sur une ron- 
delle en cuivre ou en béryllium de 0,1 à 0,2 milli- 
mètre d'épaisseur avec un trou circulaire de 3 milli- 
mètres de diamètre. Le meilleur moyen pour con- 
trôler et mesurer les coïncidences parasites de ce 
type consiste à diminuer la pression du gaz rem- 
plissant les compteurs. Pour toute rétrodiffusion 
sur une matière autre que le gaz, les conditions 
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géométriques restant les mêmes, la grandeur de 
impulsion correspondant au parcours dans le 
premier compteur avant rétrodiffusion diminue 
proportionnellement à la pression du gaz ; pour la 
rétrodiffusion sur les molécules du gaz, il s’ajoute à 
cette même raison de diminution celle due à la 
diminution de la probabilité de diffusion. 


IT. Spectre des photons. — Le domaine d’énergie 
des photons émis simultanément avec le rayon- 
nement f que nous puissions étudier avec notre 
spectrographe est limité du côté des grandes éner- 
gies par l'efficacité du compteur et du côté des 
basses énergies par l’écran qu'il est nécessaire de 
mettre devant la source pour arrêter les électrons 
et diminuer les effets parasites. Cette limite infé- 
rieure est de 1 keV. Pour descendre plus bas, il 
faudrait laisser la source nue, sans écran, et on est 
alors considérablement gêné par les électrons émis 
simultanément avec le rayonnement 8 et les effets 
parasites précédemment décrits. Par contre pour 
l’étude d’une désintégration par capture XÆ, ces 
raisons n'existent plus et rien ne s'oppose, la source 
étant nue dans l’un des compteurs, à la mesure de 
photons d’énergie inférieure à 1 keV en coïnci- 
dence avec des rayons X ou des électrons Auger. 

La figure 1 représente la disposition des divers 
supports et écrans autour de la source dans les 2 cas 
choisis pour l'étude du spectre des photons. 
Ces 2 dispositifs ont l’avantage de rendre prati- 
quement nulles, dans le domaine d’énergie étudié, 
les coïncidences parasites par rétrodiffusion. Ils 
fournissent par contre des impulsions en coïnci- 
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Dispositif A : Les écrans E, et E, sont en béryllium 
de 3,3 mg/em”. La membrane mince supportant la 
source est tendue sur une plaque de béryllium de 0,2 mm 
d'épaisseur percée d’un trou de 5 mm de diamètre. 

Dispositif B : L'écran E, est constitué par une mem- 
brane de matière plastique aluminisée de l’ordre de 
40 microgrammes/cm? tendue sur une plaque de cuivre 


de 0,2 mm d'épaisseur percée d’un trou de 10 mm de 
diamètre. La membrane mince supportant la source est 
tendue sur une plaque de cuivre de 0,2 mm d’épaisseur 
percée d’un {trou de 3 mm de diamètre. 


N°1 
dence dues à des électrons émis simultanément avec 
le rayonnement $ et franchissant les écrans, ainsi 
que des coïncidences parasites dues à des électrons 
projetés. Nous ne tenterons pas de distinguer ces 
impulsions dues à des électrons de deux prove- 
nances différentes ; il nous suffira de remarquer 
qu’un changement de la nature ou de la pression du 
gaz remplissant les compteurs fournira toujours le 
même spectre d’impulsions dues aux électrons. En 
fait, nous obtiendrons le spectre des photons en 
comparant les résultats obtenus avec les trois 
remplissages suivants des compteurs : 

a) argon : 60 em de mercure + propagane 6 em ; 

b) propane : 57 cm de mercure ; 

c) propane : 8 cm de mercure. 

En négligeant l'efficacité des compteurs par effet 
paroi, il est facile de calculer efficacité en fonction 
de l’énergie pour chacun de ces remplissages, 
compte tenu des dispositifs À et B représentés sur 
la figure { et de la forme connue des compteurs ; 
ces courbes sont représentées sur la figure 2. 


Efficacité d'un compteur 
aux photons 
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L’étalonnage en énergie des compteurs est 
obtenu à l’aide des rayonnements X (X-L) de 
l'aluminium —1,5 keV, (K-L) de Pargon —2,9 keV, 
(X-L) du manganèse — 5,9 keV, un mélange (K-L) 
et (X-M) du cuivre — 8,5 keV, (L-M) du bis- 
muth = 13 keV, (L-N) du bismuth — 16 keV et la 
raie de 47 keV du RaD. à 

La figure 3 donne les résultats bruts en coïnci- 
dence après soustraction du mouvement propre et 
des coïncidences fortuites ; ils sont obtenus dans le 
sélecteur à 50 canaux avec le dispositif B et le 
remplissage b (propane : 57 cm de mercure) pour 
des gains de l’amplificateur de la voie 1 ayant les 
valeurs 2 500 et 5 000, et pour la voie 2 un gain 
de 580 et une discrimination ne laissant passer que 


No 
les impulsions corre 
rieure à 14keV. 


La figure A représente les résultats complets dans 
Jes mêmes conditions, les abscisses étant gradués 


pondant à une énergie supé- 
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en keV et les ordonnées représentant le nombre 
d’impulsions comptées dans un intervalle de 4 keV. 
De plus, les courbes ont été ramenées à un taux de 
comptage de 10% dans le compteur 2 (NW, — 10 000 
par unité de temps, pour une énergie supérieure 
à 14 keV). Ce qui nous intéresse est de connaître le 
spectre des impulsions dues au rayonnement enre- 
gistré dans le compteur 1 et associé aux 10 000 
rayons $ du compteur 2 ; il faut donc soustraire 
de la courbe de la figure 4 les coïncidences dues aux 
rayons $ enregistrés dans le compteur 1 et associés 
au rayonnement d'énergie supérieure à 14 keV du 
com pteur 2. Ces coïncidences représentées en poin- 
tillé sur la figure 4 sont calculées moyennant cer- 
taines hypothèses : on suppose que le rayonnement 
émis simultanément au rayonnement & est indé- 
pendant de l’énergie de ce dernier, ce qui donne 
pour la courbe pointillée le même spectre à une 
multiplication près des ordonnées (7/104.N,) que 
le spectre V, du rayonnement 8 enregistré dans le 
compteur 1 ; » est le nombre total d’impulsions 
d'énergie supérieure à°14 keV dans le compteur 2 
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émis simultanément avec 10* rayons $ enregistrés 
dans le compteur 1. Tenant compte des spectres 
intégrés de N, et des coïncidences pour des énergies 
supérieures à 14 keV, et connaissant les angles 
solides sous lesquels sont vus, de la source, chacun 
des compteurs, on en déduit pour n la valeur 92 
On voit que la courbe pointillée ainsi obtenue et 
venant en déduction des coïncidences, représente 
une correction très faible pour des énergies infé- 


rieures à une dizaine de keV. Ce résultat est obtenu . 


. grâce aux grandes dimensions des compteurs ren- 
dant les valeurs de !V, très petites pour des énergies 
inférieures à une dizaine de keV environ. 

A côté des résultats représentés sur la figure 4, 
d’autres résultats semblables ont été obtenus avec 
chacun des remplissages a, b, et c du compteur ou 
avec chacun des dispositifs Aet B (fig. 1). Les in pul- 
sions en coïncidence dues aux électrons étant ind é- 
pendantes du remplissage, il est facile de les éli- 
miner par différences entre les courbes obtenues 


pour deux remplissages différents. Tenant compte 
des efficacités aux photons du compteur (fig. 2) 
nous avons obtenu le spectre des photons qui est 
représenté sur la figure 5. La courbe en trait plein 
représente le spectre des photons venant de la 
source, à l’exception de la bosse située autour 
de 9 keV qui est due à des photons d’énergie supé- 
rieure enregistrés dans le compteur par effet paroi 
(raie À du cuivre). Cette partie de la courbe a 
d’ailleurs été obtenue avec le dispositif À (fig. 1) 
où la plaque supportant la membrane porte-source 
est en béryllium. Avec le dispositif B où cette 
même plaque est en cuivre les résultats obtenus 
sont les mêmes, sauf autour de 9 keV où la bosse est 
considérablement augmentée, ce qui est dû à la dif- 
fusion des rayons 8 par les bords du trou et exci- 
tation de la couche X du cuivre. La comparaison 
de ces résultats avec plaque de cuivre et plaque de 
béryllium, ainsi que la prévision théorique per- 
mettent d'affirmer que la contribution due au frei- 
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nage externe du rayonnement f est absolument 
négligeable. La courbe pointillée de la figure 5 se 
confond avec la courbe en trait plein à l'exception 
de 2 raies, l’une à 1,5 keV et l’autre à 2,9 keV. La 
première a été obtenue lorsque la membrane en 
matière plastique supportant la source est recou- 
verte par vaporisation sous vide d’une couche 
d'aluminium d’une dizaine de microgrammes par 
centimètre carré pour la rendre conductrice. L’exci- 
tation de la couche XÆ de l'aluminium par le rayon- 
nement 8 est alors responsable de cette raie (celle-ci 
est visible sur les figures 3 et 4). La deuxième raie 
est obtenue avec le remplissage du compteur par de 
l’argon et est due à l’excitation par le rayon- 
nement $ de la couche X de ce gaz. Dans ces 
bandes d’énergie, la courbe en trait plein provient 
de mesures faites avec une source sans aluminium 
et un remplissage de propane. 

En résumé, à l'exception de cette légère bosse 


autour de 9 keV, la courbe en trait plein de la 
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figure 5 fournit le spectre des photons émis simul- 
tanément dans un angle solide 4x avec 10 000 
rayons 6. Ce spectre se compose d’un fond continu 
d'intensité décroissante en fonction de l’énergie dû 
au rayonnement de freinage interne et d’une raie 
de 2,3 keV qui est la raie À de l’atome résiduel de 
soufre et qui est due au phénomène d’autoioni- 
sation. La séparation du fond continu et de la raie 
présente un certain arbitraire qui est considéra- 
blement réduit par l'obligation d’obtenir après 
séparation une raie symétrique ayant la largeur 
convenable imposée par le pouvoir de résolution 
connu du compteur. Le résultat, après séparation, 
est donné sur la figure 6 ; la surface de la raie 
correspond pour le phénomène d’autoionisation à 
l'émission de 5,5 photons À pour 10 000 rayons £. 
La figure 6 reproduit aussi en trait pointillé le 
résultat prévu pour le rayonnement de freinage 
interne par la théorie de Knipp et Uhlenbeck [3] 
dans notre dispositif expérimental. Cette courbe a 


No1 


été calculée en tenant compte de l'intensité prévue 
par la théorie en fonction de l’énergie et aussi de la 
corrélation angulaire prévue par cette même théorie 


entre rayon 8 et photon de freinage. Cette corré- 


lation angulaire présentant un maximum vers 
l'avant, il en résulte que, avec notre dispositif A 


+ En pointillé prévision théorique 
| pour le rayonnement de frei- 
+ ; À VE 
nage interne dans le dispositif 
expérimental 


pour. 10 “rayons 


Nombre de photons par keV 


D—e 0 
ne TC —e- 2 —6- —0— + —o- 


AU dE — TR D TC 
CE E 10 e 
FIG. 6: 


(fig. 1), le nombre de photons enregistrés en coïnci- 
dence ne représente que 41 % de l'intensité que 
l’on obtiendrait avec le rayonnement isotrope. Avec 
le dispositif B, cette p'oportion est de 37 %,. (Ces 


valeurs ont été obtenues par intégration gra- 


phique.) | 
On voit sur la figureT6 qu’il existe un accord 


excellent entre la théorie de freinage interne et 


notre résultat expérimental entre 1 et 5 keV. De 5 à 
12 keV la courbe expérimentale s’élève progres- 
sivement et rapidement au-dessus de la courbe 
déduite de la théorie. Ce fait est à rapprocher des 
résultats obtenus par certains auteurs [4]. H. Lan- 
gevin-Joliot étudiant le spectre des photons de 
freinage interne du phosphore 32 jusqu’à 15 keV a 


trouvé pour cette énergie un excès de 50 % par 


rapport à ce que prévoit la théorie. Nos expériences 
indiquent à 12 keV l'existence d’un nombre de 
photons presque 3 fois supérieur à la prévision 
théorique. Si l’on veut chercher un accord entre 
nos résultats et ceux de H. Langevin-Joliot, il 
faudrait admettre que l’excès de photons constatés 


dans les deux cas ne suit pas la même loi de corré- 
‘lation angulaire que celle prévue par la théorie de 
freinage interne, mais est émis d’une façon prati-. 
quement is0trope par rapport à la direction d’émis- 


sion de la particule £. Indiquons pour terminer que 
nos résultats ne sont pas en accord avec les expé- 
riences de G. À. Renard [5] sur l’autoionisation 
et le rayonnement de freinage interne entre 1 et 
5 keV du phosphore 32. Sauf erreur systématique 
dans une des expériences, le désaccord pourrait 
provenir d’une loi de corrélation angulaire’photon- 
rayon $ différente de la prévision théorique. 
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IT. Spectre des électrons. — L'étude du spectre 
des électrons émis simultanément avec le rayon- 
nement $ lors de la désintégration du phosphore 32 
a été effectuée avec le dispositif C comme porte- 
source (fig. 7) dans lequel aucun écran n’est inter- 


Compteur 1 Laiton 


Compteur 2 


Dispositif C Plexiglas 


, F1G. 7.— La membrane mince d’environ 40 micro- 
grammes par cm?supportant la source est tendue sur une 
plaque de cuivre de 0,2 mm d’épaisseur percée d’un trou 

F de 3 mm de diamètre s 


posé entre la source et le gaz du compteur 4. Le 
dispositif C’ est celui dans lequel un écran est 
placé directement contre la source du côté du 
compteur Î. Avec le dispositif C, nous enregis- 


trerons dans le compteur 1 des photons, des élec- : 


trons et les coïncidences parasites précédemment 
mentionnées. Les résultats que nous donnons dans 
la suite de cet article ont été obtenus après sous- 
traction des coïncidences dues aux photons ; ces 
coïncidences ont été calculées en partant du spectre 
représenté sur la figure 6 et de l’efficacité connue du 
compteur (fig. 2). 

Il existe pour des énergies faibles un spectre 
rapidement décroissant de coïncidences parasites 
dues à la rétrodiffusion du rayonnement £. L’impor- 

tance de ce spectre avec un remplissage du comp- 
teur en argon rend impossible l’étude des électrons 
par ce procédé. Nous avons utilisé les remplis- 
sages b (propane : 57 cm de mercure) et € (pro- 
. pane : 8 cm de mercure). Avec le remplissage b ce 
spectre parasite est nul pour des énergies supé- 
rieures à 2 keV ; avec le remplissage c, il disparaît 
complètement pour des énergies supérieures à 
1 keV, qui représente la limite inférieure de nos 
résultats. L'emploi du dispositif C’ supprime égale- 
ment ce spectre parasite, avec le remphssage b, 
montrant ainsi qu'il est dû à la rétrodiffusion du 
rayonnement 6 sur les bords de la rondelle sup- 
‘portant la membrane sur laquelle est déposée la 
source, et non pas à une rétrodiffusion sur les 
molécules du gaz. 

Le spectre d'électrons que reproduit la figure 9 se 
rapporte à 10 000 rayons 6. Les coïncidences, 
d’ailleurs en petit nombre, dues aux rayons $ enre- 
gistrés dans le compteur 1 ont été soustraites après 
_ calcul par une méthode analogue à celle décrite 

précédemment. Comme nous l’avons dit, les coïnci- 
dences dues aux photons, d’un ordre de grandeur 
inférieur à celles dues aux électrons, ont été égale- 
ment soustraites. Les coïncidences parasites dues à 
la rétrodiffusion n’apparaissent pas car le spectre a 
été obtenu avec le remplissage b pour une énergie 

supérieure à 2 keV et avec le remplissage c pour une 
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énergie supérieure à 1 keV. (Au-dessus de 2 keV, 
les deux résultats sont identiques, aux erreurs 
d'expérience près.) La figure 9 donne le spectre des 
électrons dans le dispositif C’ après traversés d’un 
écran de 0,04 mg/cm? ou de 0,2 mg/em? montrant 
que l’absorption est conforme à ce que l’on peut 
attendre. 

S'il existe des coïncidences parasites dues à la 
projection d’électrons par le rayonnement £, ce 
phénomène ne peut se produire que dans la matière 
environnant la source (fig. 7). Cette matière est 
constituée par environ 2 à 3 millimètres d’épais- 
seur de gaz, par la membrane support de source et 
par la membrane faisant écran à la partie inférieure 
du dispositif C au contact du gaz du compteur 2. 
Ces deux membranes sont en matière plastique 
aluminisée d’une épaisseur totale de 0,04 mg/em?. 
Un calcul basé sur la probabilité de projection d’un 
électron et sur l’absorption de cet électron projeté 
nous a donné une valeur négligeable devant l’inten- 
sité du spectre de la figure 8. Ce point de vue est 
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confirmé expérimentalement par l'étude analogue 
faite sur le soufre 35, que nous développons dans un 
autre article de ce même journal. 

Sur Ja figure 8, nous avons reproduit en pointillé 
une décomposition du spectre des électrons donnant 
d’une part une raie de 2,0 keV et d’autre part, un 
fond continu décroissant avec l’énergie. Cette raie 
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Source recouverte par un écran 
c'e formvar aluminise 0,04 mg/cm° (courbe 1) 
d'aluminium de 0,2 mg/cm? (courbe 2) 
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10 keV 


est due aux électrons Auger consécutifs au réar- 
rangement du cortège électronique de l’atome de 
soufre. Lorsqu'il y a autoionisation de l’atome, 
l’électron Auger est émis simultanément avec l’élec- 
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tron éjecté de la couche X, et cet électron Auger 
apparaît dans la raie de 2,0 keV seulement dans le 
cas où la direction d'émission de l’électron éjecté ne 
permet pas qu'il soit enregistré dans le même 
compteur. Par contre lorsqu'il y a autoexcitation 
de l’atome et que l’électron éjecté de la couche X 
est lié, l’électron Auger est toujours compté dans un 
cas sur deux, le compteur 1 étant vu de la source 
sous un angle solide 27. La proportion inconnue des 
phénomènes d’autoionisation et d’autoexcitation 
ne permet pas de déduire de la surface de la raie 
d'électrons de 2,0 keV le nombre d'électrons Auger 
effectivement émis par la source. Il est important 
de remarquer que les électrons enregistrés dans 
notre expérience sont non seulement ceux émis 
dans un angle solide 2x, mais également ceux rétro- 
diffusés par la membrane supportant la source. 
Cette rétrodiffusion s'accompagne d’une dégra- 
dation de l’énergie ; il est pratiquement impossible 
de séparer la partie relativement faible du spectre 
qui correspond aux électrons rétrodiffusés ; indi- 
quons seulement que la fraction d’électrons rétro- 
diffusés dans ce domaine d’énergie par une épaisseur 
saturante d'aluminium est de 15 %, [7]. 

En résumé, nous trouvons des électrons mono- 
cinétiques de 2,0 keV et des électrons d’énergie 
variable dont le spectre est représenté approxi- 
mativement par une loi en fonction inverse de 
l’énergie avec une intensité égale à 156.104 entre 
Let 25 keV. 

La théorie de l’autoionisation [8] prévoit pour 
le phosphore 32 l’émission d’un spectre continu 
d'électrons que nous avons représenté en pointillé 


No“ 


sur la figure 10, notre spectre expérimental, qui est 
le même que celui de la figure 8, y étant représenté 
en trait plein. On constate vers { et 2 keV environ 
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3 fois plus d'électrons que n’en prévoit la théorie ; : 
à partir de 4 à 5 keV, le désaccord devient plus 
considérable, puisque l'intensité du spectre théo- 
rique pour des énergies supérieures à 5 keV est 
de 2,10-4#et que notre valeur expérimentale est de 
66.10—4+ entre 5 et 25 keV. Ce fait confirme et 
précise des expériences plus anciennes effectuées à 
laide de 2 compteurs Geiger-Muller fonctionnant 
en coïncidence [9]. 
Manuscrit reçu le 29 octobre 1958. 
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Par G. CHARPAK et F. SUZOR, 


Laboratoire de Synthèse Alomique, C. N. R. S. 


Résumé. — L’analyse des impulsions en coïncidences dans deux compteurs proportionnels a 
permis l’étude du spectre de photons et d’électrons de basse énergie (1 à 20 keV) émis en même 


temps que les rayons 8 par %5$. 


Le spectre de photons comprend la raie X de 2,7 keV provenant de l’autoionisation de l’atome 
résiduel, d'intensité 1,8.10—* par désintégration et le spectre continu du rayonnement de frei- 
nage interne, dont l'intensité est du même ordre que celle prévue par la théorie, de 4 à 10 keV. 

Le spectre électronique comprend la raie Auger de 2 keV environ, et un spectre continu, qui 
comme pour **P, mais avec une intensité moindre, est plus intense que celui prévu par la théorie. 
De 5 à 20 keV, on trouve 3,6.10—% électrons associés par désintégration. 


Abstract. — By analysis of the coincident impulses in two proportional counters, study is 
made of the low energy photon and electron spectrum (1 keV to 20 keV) emitted simultaneously 


with beta rays by 558. 


The photon spectrum comprises 2,7 keV K:-line following auto-ionisation of the residual atom, of 
intensity 1:8 X 10—* per disintegration and the internal bremsstrahlung continuous spectrum of 
intensity of the same order of magnitude as predicted by theory, from 4 keV to 10 keV. 

The electron spectrum comprises the Auger line of about 2 keV and a continuous spectrum 
which, as in *?P, although with smaller intensity, is much higher than predicted by autoionisation 
theory. From 5 keV to 20 keV 3,6.10—% associated electrons are found per disintegration. 


I. Spectre des photons. — Une étude analogue à 
celle décrite dans l’article précédent sur le phos- 
phore 32 a été entreprise pour le soufre 55. La diffé- 


rence de comportement de ces 2 radioéléments dans: 


notre spectrographe, due aux énergies très diffé- 
rentes des rayonnements f nous a contraints à 
employer d’autres dispositifs expérimentaux. 
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Laiton 


2 
Rene Dispositif E 
PR  ——— 
F1G, 4: S 
Dispositif D.— L'écran E,; est en béryllium de 
3,3 mg/em?. La membrane mince supportant la source 
est tendue sur une plaque de laiton de 0,2 mm d’épaisseur 
percée d’un trou de 3 mm de diamètre. ns 
Dispositif D’. — Identique au dispositif D sauf 
écran E, qui est en béryllium de 6,6 mg/cm?. 
Dispositif E. — La source du côté du compteur est 2 
= est déposée sur une plaque de béryllium de 6,6 mg/cm 
et de diamètre 18 mm. 


C’est ainsi que le dispositif À (partie I, figure 1) 
ne pouvait convenir, le rayonnement 8 du soufre 35 
étant trop absorbé par l'écran de 3,3 mg/em? de 
béryllium situé du côté du compteur 2. Nous avons 
employé le dispositif D (fig. 1) ne présentant pas 
cet inconvénient. 

Les mêmes remplissages en gaz des compteurs 
que pour l’étude du phosphore 32 ont été utilisés ; 
les efficacités des compteurs dans ces différents cas 
ont été données dans la figure 2 de la partie I 
précédente. 

La figure 2 représente les coïncidences obtenues 


3 Argon 60 cm de mercure 
= Propane 6cm ,, 5 
< Dispositif D 
e L ? 
is N, = 10000 
8 
5 
= ei se 
Q | 
= | ; — 

| (2 

(ERtere are 

Ci 6 SMS AIDE 2AkeV 

(IRTG. 2: 


après soustraction du mouvement propre et des 
coïncidences fortuites dans le dispositif D avec le 
remplissage a {argon 60 cm de mercure — pro- 
pane 6 cm) ; la raie très importante:de 3,0 keV due 
à l’ionisation de l’argon par le rayonnement f rend 
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impossible toute étude en dessous de 4 keV avec ce 
ren plissage des compteurs. 
La figure 3 représente le spectre des photons 


En paintillé prévision théo- 
\ rique pour le rayonnement 
de freinage interne dans le 
dispositi 


experimental 


= 
mm me 


Nombre de photons par keV 
pour 10* rayons 


entre 4 et 10 keV compte tenu des angles solides, de 
l'efficacité et du pouvoir de résolution des comp- 
teurs. Ce résultat a été obtenu avec le remplissage a 
d’argon et de propane et les dispositifs D et D’. Le 
dispositif E qui favorise le freinage externe dans le 
gaz à cause de la surface beaucoup plus grande de 
l’écran séparant les compteurs, et également dans le 
support de béryllium à cause du contact direct de 
celui-ci avec la source fournit, dans le même 
domaine d’énergie, un nombre environ double de 
photons. Il résulte de cette comparaison que les 
photons de freinage externe dans le gaz doivent 
avoir une influence négligeable dans les dispo- 
sitifs D et D’. L'identité des résultats obtenus avec 
le dispositif D (E, — 3,3 mg/em? de Be) et avec le 
dispositif D’ (EE, — 6,6 mg/cm? de Be) permet de 


Propane 57 cm de mercure 


Dispositif D 


W, = 10000 


Coincidences par keV 
a 
(éa] 


10 keV 


penser que l’influence des photons de freinage dans 
le béryllium est aussi négligeable. Étant donné le 
nombre très faible des photons du spectre repré- 
senté sur la figure 3 et l’importance relative très 
grande des corrections dues au mouvement; propre 
et aux coïncidences fortuites, il ne faut considérer 
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celui-ci que comme un ordre de grandeur d’une 
précision approximative de 30 % en bon accord 
avec la prévision théorique représentée en poin- 
tillé ; comme dans le cas du phosphore, nous avons 
tenu compte pour cette prévision théorique non 
seulement de l'intensité prévue par la théorie de 
Knipp et Uhlenbeck mais aussi de la corrélation 
angulaire entre rayon Ê et photon de freinage. 

La figure 4 représente les coïncidences obtenues 
(courbe 4) avec le dispositif D et le remplissage b 
(propane : 57 em de mercure) ; dans le compteur 2 
ne sont comptées que les impulsions d'énergie supé- 
rieure à 14 keV. Cette expérience nous permet 


d’atteindre avec précision la raie d’autoionisation 


de 2,7 keV de l’atome résiduel de chlore. Pour cela 
nous devons soustraire de la courbe 1 les coïnci- 
dences dues à d’autres phénomènes et illustrées par 
les courbes 2, 3 et 4. La courbe 2 donne la contri- 
bution du freinage interne tel qu’.l est représenté 
sur la figure 3 entre 4 et 10 keV, l’extrapolation en 


dessous de 4 keV étant faite conformément à la 


théorie. La courbe 3 représente les coïncidences 
dues au rayonnement $ compté dans le compteur 1 
et associé à des électrons ou à des photons émis 
simultanément dans le compteur 2 ; cette courbe 3 


= — HS 
a [éa] [æ) 


=: 
(Cal 


Nombre de protons par keV pour 10“ rayons p 


1 
5 keV 


F16. 5.— En ordonnées, lire photons, et non protons. 


est dessinée homothétique au spectre du rayon- 
nement $ fourni directement dans le compteur 1 et 
telle que la somme des courbes 2 et 3 redonne la 
courbe 1 dans le domaine compris entre 7 et 10 keV. 


La courbe 4 représente la raie X de l’aluminium 


(1,5 keV) due à l’ionisation par le rayonnement B de 
la couche d’aluminium déposée sur la membrane 
supportant la source ; cette courbe 4 a été dessinée 


compte tenu de la nécessité d'obtenir pour celle-ci 


et pour la raie d’autoionisation des courbes symé- 
triques présentant une largeur compatible avec la 
résolution connue du compteur. 


as Diamant antbmafton! 


noté min sénnnimetaititoumests its ambre aps 
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Tenant compte en outre de l’angle solide, de 
l'absorption par l’écran et de l'efficacité du comp- 
teur, nous obtenons la raie d’autoionisation repré- 
sentée sur la figure 5. La surface de celle-ci donne 
pour les photons À du chlore une probabilité 
d'émission de 1.8.10-4, Ce résultat est un peu infé- 
rieur à celui trouvé par Howland et Rubinson [1] 
qui est de (2,47 + 0,6) 104 ; il faut remarquer que 
ces auteurs devaient soustraire du résultat brut 
qu’ils obtenaient un fond continu dont l’impor- 
tance relative était grande, alors que dans notre 
cas en ne soutrayant aucune des corrections figurées 
par les courbes 2, 3 et 4 et en attribuant ainsi à la 
raie d’auto-ionisation la totalité des photons dispo- 
nibles, on arriverait seulement à 2,1.10-4, 

Il est important de remarquer que le résultat 
obtenu pour les photons d’autoionisation dans le 
cas du phosphore 32 est de 5,5.10-4et qu'il n’est 
que de 1,8.1074 pour le soufre 35. Ces deux radio- 
éléments ayant des numéros atomiques différant 
seulement d’une unité la théorie [2] prévoit des 
résultats sensiblement analogues : celui pour le 
phosphore étant supérieur de 10 à 20 %. 


IT.. Spectre des électrons. — L'étude du spectre 
des électrons émis simultanément avec le rayon- 
nement £ du soufre 35 a été effectuée avec le même 
* dispositif C (partie I, figure 7) déjà employé pour 
Pétude semblable faite sur le phosphore 32. Pour 
_ éviter les coïncidences dues à la rétrodiffusion du 
rayonnement $ nous avons seulement employé le 
remplissage C du compteur (propane : 8 cm de mer- 
cure) et par diminution progressive de la pression du 
propane dans le compteur, nous avons vérifié que 
cet effet était alors nul pour des énergies supérieures 
à 1 keV. Dans le compteur 2, le rayonnement f, 
donnant des impulsions d’énergie supérieure à 
12 keV, était seul enregistré, rendant négligeable la 
probabilité d'obtenir une coïncidence entre un 
rayon 8 compté dans le compteur { et un électron 
dans le compteur 2. Rappelons que le temps de 
résolution du sélecteur à coïncidences est connu 
avec une précision du centième, ce qui permet de 
soustraire parfaitement les coïincidences fortuites. 
Il n’est pas nécessaire de tenir compte des coïnci- 
dences dues aux photons et étudiées au début de cet 
article car celles-ci sont en quantité négligeable 
devant celles dues aux électrons. 

Il reste un dernier effet à soustraire du résultat 
ainsi obtenu, c’est l'effet parasite provenant des 
. électrons projetés par le rayonnement f dans toute 
matière environnant la source. La section efficace de 
projection d’un électron (DÉnelRLe comprise entre Q 
et Q + dOest égale à _. a d’après la formule de 
 Bohr, v étant la vitesse de l’électron incident (dans 
le cas où Q est faible devant l’énergie de Pélectron 
incident, il est inutile d’utiliser la formule plus com- 
plexe tenant compte des spins). On tire de là que 
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la probabilité de projection d’un électron par un 
rayon Ê traversant 1 micro-gramme par em? de 


de L Lee: 
matière est égale à Es a prenant la valeur 4,7 


pour le rayonnement 6 du soufre 35 dans nos condi- 
tions expérimentales et la valeur 1,2 pour le phos- 
phore 32. I] est important de remarquer que dans 
des conditions expérimentales identiques le nombre 
d'électrons projetés dans le cas du soufre 35 sera 
ainsi 4 fois plus grand que dans le cas du phos- 
phore 32. Il est pratiquement impossible à partir de 
cette probabilité de projection, des conditions géo- 
métriques de l’expérience et des lois d’absorp- 
tion [3] des électrons projetés de déduire avec pré- 
cision le spectre que l’on doit obtenir ; on peut 
seulement prévoir des limites supérieures qui con- 
duiraient pour le spectre des électrons projetés du 
soufre 35 associés dans un angle solide 4x à 104 
rayons $ à des valeurs de l’ordre de 10 vers 2 keV, 
de 7 vers 5 keV, de 3 vers 10keV et de 1 vers 17 keV, 
la réalité étant très probablement inférieure à ces 
limites supérieures. Pour le phosphore 32 ces 
valeurs devant être divisées par 4 sont comme nous 
l’avons dit négligeables devant le spectre des élec- 
trons associés au rayonnement £. 

À côté de ces prévisions théoriques nous avons 
également étudié expérimentalement cet effet des 
électrons projetés dans le domaine des énergies 
supérieures à 7 ou 10 keV ; dans ce cas en effet 
l'absorption des feuilles de 100 ou 200 micro- 


-Photons K de J’AL 


Electrons prajetes 


One 


grammes par em? est assez faible pour que le 
nombre des électrons projetés soit proportionnel à 
l'épaisseur de la feuille. Nous avons alors utilisé le 
dispositif C” qui se distingue du dispositif C unique- 
ment par le fait qu’un écran d’aluminium de 
230 microgrammes par cm? est placé derrière la 
membrane supportant la source et en contact direct 
avec celui-ci. Dans ces conditions le dispositif C 
doit nous donner le spectre des électrons émis simul- 
tanément avec le rayonnement 8 et celui des élec- 
trons projetés dans le support de 40 microgrammes 
par em? ou dans le gaz avoisinant ; le dispositif C 
doit donner en plus le spectre des électrons projetés 
par le rayonnement $ dans la feuille de 230 micro- 
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grammes par em? d'aluminium. La différence entre 
les spectres obtenus avec C” et avec C est repro- 
duite sur la figure 6 et représente l’effet dû à 
l'aluminium. Ce spectre ramené à l'émission dans 
un angle solide 4x et correspondant à 104rayons 8 
enregistrés comprend une raie d'énergie 1,5 keV 
qui correspond aux photons X d’ionisation de l’alu- 
minium et un spectre très étalé dû aux électrons 
projetés. 

Même en tenant compte d’une incertitude sur 
l'épaisseur du support de source employé, on peut 
admettre que les électrons projetés d’une énergie 
supérieure à 7 ou 10 keV seront dans le dispositif C 
en nombre au plus 3 fois inférieur à ceux repré- 
_sentés sur le spectre de la figure 6, ces valeurs étant 
d’ailleurs inférieures aux limites supérieures tirées 
des prévisions théoriques et précédemment men- 
tionnées. 


35S 


pour 10* rayons f 


Nombre d'électrons par keV 


ETGE 


Nous avons tracé sur la figure 7 le spectre des 
électrons émis dans 4x et associés à 104rayons 8 du 
soufre 35 après soustraction pour des énergies supé- 
rieures à 7 ou 10 keV du tiers du spectre des élec- 
trons projetés de la figure 6. Comme dans le cas du 
phosphore 32 ce spectre se décompose en un fond 
continu d'intensité décroissante avec l’énergie et 
une raie voisine de 2 keV. Celle-ci est due aux élec- 
trons Auger d'énergie 2,3 keV consécutifs au réar- 
rangement du cortège électronique de l’atome de 
chlore et a une surface approximative de 5 pour 104 
rayons £. 

Nous avons représenté à nouveau sur la figure 8 
le spectre continu d’électrons associés au rayon- 
nement $ du soufre 35. Celui-ci présente le même 
aspect que celui obtenu dans le cas du phosphore 32, 
avec cependant une intensité inférieure pour des 
énergies supérieures à environ 10 keV. L’intensité 
entre À et 5 keV est égale à 81 et celle entre 5 et 
20 keV est égale à 36 pour 104 rayons 6. Nous 
avons tracé en pointillé sur la figure 8 le spectre 
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continu d'électrons prévu par la théorie de l’auto- 
ionisation [2]. Comme dans le cas du phosphore 32, 
il existe un désaccord flagrant surtout pour les éner- 
gies élevées de ce spectre, résultat confirmant les 
expériences de Duquesne [4], effectuées également 
par la méthode des coïncidences. 


En pointillé 
prévision théorique 


Nombre d'électrons par keV pour 10* rayons $ 


— 
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L'existence de ce spectre anomal d'électrons pour 
les deux radioéléments 8 *P et ‘55 nous a incités à 
effectuer la même recherche sur le fer 55 se désin- 
tégrant par capture Æ. Le spectre théorique des 
électrons d’autoionisation variant dans ce cas de 
15 à 1.105 par keV dans le domaine de 1 à 20 keV 
est beaucoup trop faible pour être mesurable dans 
notre spectrographe. Nous n’avons trouvé aucun 
électron associé pour le fer 55 à la raie K de 5,9 keV 
avec comme limite supérieure la valeur 3.10-4dans 
la bande de 1 à 7 keV, celle de 0,5.10-4 pour la 
bande de 10 à 14 keV et 0,5.10-4 pour la bande 
de 15 à 20 keV, ce qui nous permet de conclure 
qu’un spectre anomal analogue n'existe pas pour le 
fer 95. 

Nous sommes très reconnaissants au regretté 
PT F. Joliot pour l'intérêt avec lequel il a suivi nos 
expériences et pour les discussions fructueuses que 
nous avons eues avec lui. Nous remercions 
M. J. Merinis pour la fabrication des sources radio- 
actives. 
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ÉTUDE THÉORIQUE DE LA DISTRIBUTION ÉLECTRONIQUE 
DANS UN PLASMA LORENTZIEN HÉTÉROGÈNE ET ANISOTROPE 


Par Raymonp JANCEL et THéo KAHAN, 
Institut Henri-Poincaré, C. N. R. $.. Paris. 


… Résumé. — Les auteurs étudient un développement général de la fonction de distribution 
électronique en fonctions sphériques, en vue d’établir une solution de l’équation de transfert rela- 
tive à un plasma faiblement ionisé, hétérogène (existence d’un gradient de densité) et anisotrope 
(présence de champs électromagnétiques). On établit ainsi un système d’équations permettant de 
calculer les deux premières approximations de cette fonction de distribution. 


.Abstract."— The authors use a general expansion in spherical harmonics for the electronic 
distribution function to obtain a solution of the transport equation for a weakly ionized inhomo- 
geneous (existence of density gradient) and anisotropic plasma (in the presence of electric and 


magnetic fields). 


This solution gives rise to a system of differential equations, and the solution 


of this system of equations gives the first two approximations to the distribution function. 


On sait que le développement de Chapman-Enskog 
appliqué à un plasma faiblement ionisé soumis à 
l’action d’un champ électrique cesse d’être valable 
dès que le champ électronique devient relativement 
intense ; en effet, les électrons sont fortement accé- 
lérés par le champ durant leur libre parcours alors 
que l’énergie perdue lors d’un choc contre une parti- 
cule lourde sera de l’ordre de 2 fois le rapport des 
masses, c’est-à-dire de l’ordre de 10%. Les élec- 
trons libres gagnent donc plus d’énergie qu’ils n’en 
perdent au cours des chocs ; il en résulte un accrois- 
sement sensible de la température électronique et 
la fonction de distribution des vitesses électroni- 
ques s’écartera notablement de la répartition max- 
wellienne. Le calcul de cette fonction de distribu- 
tion électronique ne relèvera donc plus de la théo- 
rie de Chapman-Enskog qui n’est valable que pour 
un système s’écartant faiblement de léquilibre 
thermique et qui consiste à développer la fonc- 
tion de distribution en fonction d’un paramètre À: 


f= PR + P + XP 4 (1.1) 


2. Développement en fonctions <phériques. — 
Pour un champ électrique d’intensité quelconque 
appliqué à un plasma, la symétrie sphérique du 
problème suggère de faire appel à un déveioppe- 
ment de la fonction de distribution électronique j, 
en harmoniques sphériques dans l’espace des vites- 
ses sous la forme : 

+oo +1 


f= » 2 
1=0 Mm=—1 
où 0 et w sont les angles de v en coordonnées sphé- 
riques dans l’espace des vitesses. Ce développement 
tout en permettant de mettre en évidence les 
anisotropies liées à la présence des champs élec- 
trique et magnétique, ne préjuge en rien du carac- 


aj(r, V,t) 1 (V,0,0) (2.1) 


tère maxwellien ou non maxwellien de la partie 
isotrope aÿ(r, V, {) de la fonction de distribution 
électronique. Nous utiliserons pour Y? les fonctions 
sphériques sous la forme : 


YT = CAVERNE 


4T(l + m)|! 
sin” 0 oo ENT Pi(cos 0) eimp (2.2) 
avec la relation 
Vrai Tr (2.3) 
Les conditions de réalité exigent : 
a = (— 1) o ?. (2.4) 


Supposons maintenant que notre plasma soit 
soumis à un champ magnétique extérieur cons- 
tant B qu’on prendra pour débuter parallèle à Oz, 
ce qui n’amènera aucune perte de généralité. Les 
électrons sont alors soumis à une force spécifique : 


Ts = (VA B); (2.5) 
en posant 
Be (2.6) 
m 
on à : 
Te = wB /A\ Ve. (2.7) 
Si 
TER (2.8) 


m 


est l’accélération subie par l’électron du fait du 
champ électrique E(r, t), l'équation de Boltzmann 
prendra la forme 

Djelr, V, t) 


D FR) 


= J(fe) 
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et où D /Ot est l'opération boltzmannien : 


o) 
| ; + Ve.Vr + IF + (wm A vi].Vr (2.10) 


Mt D 
où J est l'opérateur de choc de Boltzmann. Si le 
gaz est lorentzien, J(f) ne comporte pas les chocs 
élastiques électron-électron. Il ne comprendra que 
les chocs élastiques électron-particules lourdes ; 
dans ce cas c’est un opérateur qui s'écrit : 


Je) = 2 f [ (fe fe — fe fo) gen b db as (241) 


où fe et f, sont les fonctions de distribution pour les 
électrons et les particules lourdes (ions, molécules) 
respectivement. En nous référant aux propriétés 
des systèmes d'équations de Boltzmann couplées 
il est possible de décrire un gaz lorentzien avec 
la seule équation (2.9) relative à la composante 
électronique, la composante lourde ayant une 
distribution maxwellienne : 


3/2 —m, v2/2XT 
REA Lin BP, (2.19) 
Si le gaz est considéré comme parfaitement 
lorentzien, c’est-à-dire si M/m, — 0, le terme de 
chocs (2.11) se simplifie en remarquant que f, —f ; 
on peut donc décomposer les intégrales en écri- 
vant pour l’opérateur J, parfaitement lorentzien : 


Ji(fe) — or [fs avs (fé 1e) 8er db. db (2/18) 


— 27 np [ue — fe) Le» b db 


avec 


gep — |(Ve — Vp)|. (2.14) 


Comme mem —0, on pourra remplacer ge 
par PV, ; d’où pour (2.13) 
Jn(fel = 27 m9 Ve [lle — fe) b dé. (215) 
Nous voici ramenés alors à un opérateur J indé- 
pendant de f, et agissant sur la fonction f, ; il est 
manifestement linéaire ses fonctions propres sont 
précisément les fonctions sphériques Y} qui ser- 
vent de base à notre développement. Si l’on dés- 
igne par v(v) la valeur propre associée aux fonc- 
tions propres Ÿ?, on a les relations (voir appen- 


dice) 
Ji(XT) = — (Ve) YF, (2.16) 


où les v, sont définis par les relations : 
CO 
TE noVe [ 1 — Pilcos x] bdb. (217) 
La quantité v, s’interprète comme la fréquence 


de collision électron-particule lourde pour les élec- 
trons de vitesse v. Elle dépend essentiellement, par 


No 1 


l'intermédiaire de la déviation x, de la loi d’inter- 
action électron-particule lourde. 

Dans ce cas particulier, la valeur propre corres- 
pondant à / — 0 est 


y NU 


Dans le cas où m./m, est petit mais différent de 
zéro, -ce qui est le cas physique, les résultats précé- 
dents ne sont qu’approchés. En portant le dévelop- 
pement de f, dans (2.11), l'opérateur de choc J(fe) 
se décompose en : 

dr). 


/ © 
Jo Y5) + 4 D 
1=1m=—1 

Dans le développement effectif, nous ferons 
l'hypothèse explicite que le premier terme J(a«ÿ Y6) 
se rapporte à un plasma imparfaitement lorentzien 
(I. L.) et que les termes avec ! > 0 décrivent un 
plasma parfaitement lorentzien (L). D’après ce qui 
précède, l’opérateur J pour ! > 0 peut donc être 
représenté par le J, de (2.15) qui, appliqué aux Y?7 
donne les résultats (2.16). On aura donc : 


(2.18) 


Jai VE) = — v(Ve) à YT (2.19) 
et par conséquent (2.18) s’écrira : 
3 CO +1 Fe ñ 
Jir(aÿ Yo) — 2 ZX wlVe) ai (Ve) Y1: (2.20) 
1=1 Mm—=—1 


Reste à définir ce que nous entendons par le 
terme « imparfaitement lorentzien. » J;r : on le 
calcule à partir de la formule exacte (2.11) en se 
bornant aux termes du 1 ordre en m,./m, ; ceci 
revient à tenir compte d’une distribution des 
vitesses v, des particules lourdes et des échanges 
d'énergie cinétique entre les électrons et ces parti- 
cules lourdes au cours des collisions (en effet, si 
on ne tenait pas compte du terme du 1% ordre 
en Me/M», le terme de choc donnerait une contri- 
bution nulle pour le 17 terme du développement). 
Dans cette manière de procéder les termes 
avec 1 > 0 mettraient en jeu des puissances 
de me/m, plus grandes que l’unité, si on apportait 
à J, la correction que l’on fait dans Ji1.. 

On effectue ce calcul en tenant compte sépa- 
rément du recul de la particule lourde et du chan- 
gement de f. du fait de l'existence d’une distribution 
des vitesses pour la particule lourde. Pour évaluer 
ces deux contributions, partons de la formule 
reliant les vitesses initiales et finales des 2 parti- 
cules au cours d’une collision, qui peut se mettre 
sous la forme suivante 


F(Ve) + RME 


On a pour le terme en J;5, de f?, d’après des résul- 
tats classiques [1], [2]: 


(2.21) 


ST Ne eut cou ©; 2AT AT 


7 à PF LARES 
Me + My r. ave LE (té Mme d(Vä) 
(2.22) 


N°1 


D’après (2.11) et (2.20), on peut donc écrire le 


terme de collision sous la forme ÿ, a 


Ji) = Fruit US À 


1=1 m=—1 
(| Me ie Curl 2kT d(axÿ) 
CUT) 
VrMe + M Pure CAN LE Sen ET 
CO ot 
T2 >) je Gi -Yre (205) 
=im—-— 


3. Application de l’opérateur ( /Ot sur le déve- 
loppement. — Le développement du terme de 
collision étant ainsi achevé, 1l nous reste à calculer 
ne de l’opérateur 7/1 sur le développement 
(2.1). 


On a sucessivement : 


dfe nr V a da. m 
9 Re Fe 
CO +1 _ 
b) Ve-Vrk = 2 YlVe.Vroi] (3.2) 
=1m—=— 
puisque Ÿ}° est indépendant de r. 
CO +1 
c) 11 vele = r.2 D  Vrelar Nr) (3.3) 
etenfin 
CO +1 
d) (ws AVe).Vrke = (8 A Ve). Z 2 Vreli YF). 
=im—— 
(3.4) 


On voit, d’après les expressions précédentes, que 
l’on est amené à calculer d’abord l’expression V, Y} 
en fonction des Y?. 


3-1. Cazcuz Des V, Y}. — On a, en utilisant les 
notations et formules de l’appendice 


1+1 


HARREE _—- sin 0.07 ei cos © 
T 
piti 
_ Te sin 0.@7 [eitm+1e L eitm—1)e] 
PAV2T 
DAS : 
_ = d— GT C(, m) + OÙ" D(I, m)] eitm+ De 
2 V527T 
+ [OŸit E(I, m) — OÙ F(I, m)] eitm—ne } (3.5) 
soit, d’après la définition des Y? 
m+i 2 wym+l 
Pan XP = — TEE OU, m) + EE DU, 
M—1 M—1 
* Hi TRUE 22 FU, m) (3.6) 
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On obtient, de la même manière 


Fe  . Ve FE 
Vs A Cr to te D{I, m) 
ire Ve Yi 
on mA, M) RÉ F(l,m) (3.7) 
et 
Vez YT = Yis1 AU, m) + VE Yi B(l, m). (8.8) 


3-2. CALCUL DE Vr, Y}. pour évaluer les expres- 
sions (3.2)-(3.4) nous avons également besoin de 
calculer Ve, Ÿ}* en fonction des Y*. En nous 
reportant toujours à l’appendice, on a succes- 
sivement 


” vas 
= ee Yr — = [ OT sin 0 ee cos © 
d O0 . SU ME 
vmp = CY10 
ce CS 0 cos æ Un = ge + sin gl 
pit : ï lo Hd m 
= 2 aient ne [ @j" sin 0 + do cos 0 — m di | 
2V9r d0 sin 0 


m Un 
+ eitm—1)e [ Of sin 0 + D cos 0 + m Si | L (3.9) 
d0 sin O1 \ 


Posons pour abréger 


| PER d O7 (OH 
a): all, m) = [ ®7 sin 9 + 50 cos 8 MS 


] (3.10) 


M mr 
b): bl, m) = [ @7" sin 0 + LES cos 0 + 1 sl: (3.11) 


Ces expressions se calculent en utilisant les for- 
mules de l’appendice ; il vient ainsi : 


; m nm d 7 
aller ( sin Ô ca (10 Sos 0 
n 2 
== (: sin 0 - ne m 1) @ 
sin 0 sin 0 


L VE m) (+ m +1) *! 


= (1 — m) sin0 O7 


+ VI m) ( + m +1) O7 cos 8 


cos 0 


= (1—m) [— OÙ CU, m) + OT DU, m)] 
LV m) ( Em +1) [074 AU, m + 1) 
(3.12) 


HO BL me At 

En utilisant les formules de définition des A(/, m), 

B{(1, m), ete. l'expression précédente se simplifie et 
l’on a : 


all, m) = (2 + 1) D(, m) OF". (3.13) 


Par un calcul en tous points analogue, on trouve 
également : 


b{t, m) = — (2 + 1) F(, m) OÙ x. (3.14) 
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En portant (3.13) et (3.14) dans (3.9), il vient : 


Pre 
ù zu Ve m+l ; 
JA TU — — (27 1) D(I, O7 eitm+l)e 
de 0 re + 1) D(E, m) Or 
— (21 + 1) F{I, m) OT? eitm—De] 
21 +1 


A Ce 0 ne MEN) 


(3.16) 
2 { = 
= Y — = re à, [D(E, m) YÉA! + F(, m) Yi] 
d 


ve LS (2 . 1) B(L, m) Ep (OA) 


d Vez 


3-3. CALCUL DE V..V. fe. — On a, par définition 


CO +1 ( Dar Doi ou 
me { >: 
Ve Vr le — em ne Yi | Ver dx Vey dy Vez 32 ) 


et, en vertu de (3.6) et (3.7), il vient : 


Ve Vrfe — DE ; (Es 2) 


(—= YEN CU m) + VE ST DU, m)| 


+ (RE + LS) CHER EU m) — Ve ES FA, m))] 


Se 2 (Yir1 AU, m) + VE Yi Bll, m)) : (3.18) 


En changeant les indices de sommation, on a 
enfin : 


Ve .Vr fe 


m \ 1 | Do" dat x À 
2 : CURE 3 
2 Lil tr 1) À _ He 
nt bn 
+ Di +1,m 1) ve( Se 1 du ) 
: mit | wEs) 
+ EU 4, m + 1) | 5x t 7 
f , 2 dote re 
F(l+1,m +1) ra ( RU . 
dort 2 dot+1 
+ Afl—1, hr + Vé B(I +1, (Or a 5410) 


8-4: CALCULE DE LV pe fe — On a, par définition : 


RAT PAS SN Re 2 (YE Vre 0 + of Vre YŸ). (3.20) 
_ Comme, d’après le développement, «7° est fonc- 
tion de V, seulement, on a : | 


Ve do 


ve Æ ee — 
Vreu — AA (8.21) 


d’où pour (3.20) : 
T:WVrefe — > (re LV + TVreYf) 


Ve dV4 

= Na Ps . (8.22) 

avec : 
Pr. m [Lx Ver + Ty Vey ue l',. «) es 3.99 
W: => [ri ( LA dPe ; (3. ) 

D dE dYr | 

RE — Ur LED D LERTe ; DE 
w.=E [4 (È ee on (22 


En tenant compte des formules (3.6), (3.7) et 
(3.8), on trouve pour W, : 


HR APN Cm) ET ED 
ER 0 A à 2 el tn) 
1 do M 2 m 
LE D (Vie Am) VAT ERP NERO) 
Ve d Ve : 


et, en regroupant les indices de sommation, il 
vient finalement : 


x S YT äl ; do" 
Le F. “hu cu ( 2 Cl — 1, m—1) 
m—1 À 
L pa PEL D(I +14, m 1) 
(e e 
Dore 1 
+ (lo + il) | ELA mA) 
v2 date FI = : 
8 Sy. (SEE EE) 
Dates 2 dA+1 } 
LT A(I pl : 
(SEE AU 4m) + EE BG + 1, m)) 


(3.26) 


Il reste à calculer W,; on utilise à cet effet les 
formules (3.15), (3.16) et (3.17) qui, portées 
dans W,, donnent après des transformations du 
même genre que précédemment 


RAT AD ee 
W:; =Ù (2 +3) tot +4, m1) ia 


1,m 2 
—— FU +1,m +1) «ri 

+ T, BÜ F1, m) al. (3.27 

3-5. CALCUL DE (wz /\ W).Vrefe — Nous 


ferons ce calcul en deux étapes ; nous commen- 
cerons par supposer que le champ magnétique B 
est parallèle à l’axe des z. Nous obtiendrons ainsi 
des formules de validité restreinte qui permettent 
néanmoins d'obtenir au premier ordre des résultats 
valables pour une orientation quelconque du champ 


. 
£ 


No 1 


magnétique ; nous confirmerons par la suite la 
validité au premier ordre de ces formules en faisant 
le calcul pour une orientation quelconque du champ 
magnétique. 


3-5.1. Cas où B est parallèle à Oz. — On a alors : 
(wB /\ Ve) Vye Î 
= À { (w8 A Ve) [PP (Vre d') + af Vre YF] }. (8.28) 


. D’après (3.21), le premier terme de cette expres- 
sion est nul et l’on a : 


(wB A Ve) .Vre f — 2 d'(wr /\ Ve) Vre Y?. (3.29) 


D’après (3.15), (3.16) et (3.17), on a (en posant 


Op — |oz|) : 


(wB8 /\ Ve .Vre fe — oB » a ( Veu 
= BE TC DU, m) ÉÉIVe + Von 
— F1, m) FE (Ver — Vey)] (8:30) 
d’où l’on tire, en raison de (3.6) et (3.7) : 
(wB /\ Ve).Vre fe — io8 2 moÿ Yi. (3.31) 


3-5.2. Cas général : orientation de B quelconque. 
— Comme précédemment, on doit calculer l’expres- 
sion (3.29) qui s’écrit maintenant 


= 
(avec Op — (oz, Qy; &,)) 


> ai (8 A Ve) .VPre Yi 
,m 


(oy Vez — 


En vertu des formules (3.15), (3.16) et (3.17), 
(3.32) s’écrit encore : 


1m 
Ë [D(, m) YÉA! — FI, m) YE) (oy Ver — 4 Ve) 


—i{D(, m) YÉT + F(, m) VE.) (we Vez — ©x Vez)] 
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Pme, Fe; to, Vex) 
=D +1 x 
1,m 


D(L, ; 
re PP on + 107) — ol Vey + iVox)] 
FS FI, m) 


2 Hola — 10z) — oz Vey — 1Vex)] 


DOTE Gt 0, ee (3.33) 


En vertu des formules (3.6) et (3.7), (3:33) s’écrit 
encore : 


DL, m) 
D ne en) lo to 
1m ? 
PCA Er Ben 
F(I, 2e 
re A RU 2 4 


A 


Ye Pl lsmedli) 


7m 
ne 1 | à E(l—1,m+1) 
ve 
1 pale FU — 1, m + 1)| 
L 
F(I Fo 
Fe on | : CUS hr ET) (3.34) 
M 
Se = ENT m—1)] 
+ B(l, m) ox X 
pri Va YP+ 
[- 7: CL —1,m) + — ni D{I — 1, m) 
Vs Vera | 
Re E(l—1,m) + TE F(l—1, m) 
— Pl, m) oy X 
m+1l 2 vuMm+l 
[- LE C(l—1,m) + ere = D(i—1,m) 
mM—1 2 M—1 
VEUT eee) Pere = FUN, | , 


soit en regroupant les indices de sommation 


(21 + 14) F(E, m + 4) re aa) PR Al —1, m) 


(3.35) 


L vP} (2 Te A 1 : 
1,m 2 
2 ci S pti 2,m-1)loy + io) Véoite BI +1,m) 
2 à F(L + 9, m + 1) (oy — io) Ve es B(I +1, m) 


+ me PET) ip, m) El — 1, m + 1) + F(, m) CU — 1, m—1)): 
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Lo PERS om, pa DIE + 2, m) FU + 1, m + 1) FU + 2, m6) DU Fm 4) 
21 


Bll;m—1){(21+1) No ro 


BU Le 20m 1) ET VE feet DIT Em 1e te) 


PAS ot : 
— BI, m + nr DAS © + 4) (oy —1ox) 


LS] 


(27 
9 


4 


à aol F(I + 1, m + 1) (oy — io) | 


L Bil +2,m +1) 


= vr'| (toz + oy) FE . 1) af l(D(I, m—1) AU —1,m) + BI, m—1)C(l—1,m—1)) 


Le 245 ee ATAU + 2, m—1) B(I + 1, m) BI 2; m 1) DIT + À, m—1))] 


2 


SERRE) É . Lip mm ce 1) BU mt UNE PU, mc 4) AU =, m) 


ae S mti Vê(E(L + 2, m + 1) BI + 1, m) — BI + 2, m + 1) FI + 1, m + 1))] 
RO Ne | | ; ee 
De + on Ui (D(I, m) EI — 1, m +1) F{l, m) C(l— 1, m—1)) (3.36) 
io ES os VE D(I + 2, m) FU + 1, m + 1) + FIL + 2, m) DL +4, m—1) . 


Il nous faut maintenant évaluer les six coef- il vient, en utihsant les définitions données dans 
ficients en A(l, m), etc. dans les parenthèses ;  l’appendice : 


D{(I, m— 1) AU —1,m) + B(I, m—1) CH —1,m 1) TETE (t + m) (Ü— m + 1), 

D({I + 2, m—1) BI + 1, m) — B(I + 2, m—1) D{I +1,m—1) = 0 

E(T— 1; m +4) B(, m+1) Æ F{li1m + 1) A1; m)— 5 = 7 VAT m)(l+m ve 4); 

F(L + 2,m + 1) B(L LA, m) — BI +2, m + 1} FU 4, m +1) =, (3.37) 
5 : D _2m “3 

D(t, m) E(l —1,m + 1) — F{I, m) C( — 1, m—1) Hi 


F{(I + 2, m) D(I + 1, m—1) — D(I + 2, m) F(I +1, m +1) = 0. 


Avec ces résultats (3.36) s’écrit alors : 


: m—1 
(wz À Ve) .Vpye fe — 2 Ven (Lox + Oy) _ AU + m) (l m +1) 
m+1 k 
— (oy — 10%) ; V(E— m) (1 + m +1) + iomad |. (3.38) 
4. Équations générales pour les a. — Ayant générales en utilisant les résultats de (2.41), (3.1), 


ainsi effectué le développement de chaque terme de (3,19). (3.26), (3.27) et (3.31) (dans le cas où B est 
l'équation de Boltzmann en fonction des Y?, ja parallèle à Oz ou (3.38) (pour B quelconque). 
fonction f.(r, V:, t) sera connue si l’on peut calculer 

les divers coefficients &"(r, V., t) de son dévelop- . 4-1. ÉQUATION RELATIVE AUX x} DANS LE CAS 
pement. Les équations auxquelles doivent obéir ou B EST PARALLÈLE A Oz. — Il vient alors en 
les &” sont obtenues en identifiant terme à terme utilisant (3.31) : 
les coefficients des Y}. On obtient les équations 


ee ce fe a Eponge né jo ts 
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do 1 % dar Dar ; À m—1 à m—1 
. CM re se ie ) ED +1,m—1) pa ( — ES _ 
; ! œ y 
: og es d\ S.Mm+1 PALNS 
+ E( —1, m + nl er x) — (I + 2e OUAIS DAT 
| { = l dy (+4, m + 1)pÉ à LE 5 
nm : m 
+ Al —1, m) _ + Ve BI +1, mn) PE | 
03 dz 
ii dar dati. \ 
+7 fem SA C(l—1,m—1) + VE NA Dame TA)) (4.1) 
Dr * M +1 
ET EU EM» 1 Re Le 
(Lx + il y) (E m + {) 2 Ve F(l +1, m +1) NA 
3 ; ù LU ; m 
+ DA —1, m) A + VEBU + 1, m) | 
Se + ir 
+ er +3)| = 5 D(I +1,m 1) «+1 . : TFU+L1,m- 1) «47 + TV, B(I + 1,m) dl] 


Lim oB a + ma — 0. 


Il y a lieu de faire deux remarques concernant 
cette équation : Pour ! — 0,il faut remplacer le 
terme de choc vw «œ* par l'expression donnée 
en (2.22) qui tient compte de l’effet de recul de la 
composante lourde maxwelliennement distribuée 
sur la fonction électronique. D’autre part, pour une 
valeur de / donnée, on a autant d’équations qu’il ya 
de valeurs possibles pour m. Toutefois lorsque 
indice supérieurm’ de «est plus grand en valeur 
absolue que l’ (c’est-à-dire lorsque [m'| > l'), le 
coefficient «° est identiquement nul d’après la 
définition même du développement en fonctions 
sphériques. 


4-2. ÉQUATION GÉNÉRALE DES COEFFICIENTS œ}!. 
— (Orientation de B quelconque.) On l’obtient 
immédiatement à partir de (4.1) : il suffit de rem- 
placer dans (4.1) le terme dû au champ magné- 
tique im oz a} par l'expression (3.38). 


5. Développement de f, sous forme tensorielle 
(irréductible). Considérons à nouveau le dévelop- 
pement général : 

Ce +1 
Re os r) (5.1) 


1=0 \to = 4 


et remarquons que les Ÿ? constituent, pour un 
donné, un tenseur sphérique irréductible de rang l, 
dont les 2! + 1 composantes se transforment sous 
une représentation à 2] + 1 dimensions du groupe 
de rotation À selon la formule : 


RYŸ R— = DE Dim (a, B, Y) VE, (5.2) 
Uo m 
où l'opérateur À défini par : 
R = exp [— 10 n.L] (5.3) 


est l’opérateur de rotation, L l’opérateur moment 
cinétique et n le vecteur autour duquel s'effectue 
la rotation d'amplitude 0, 

Comme le premier membre de (5.1) est un sca- 


laire, il en est de même du second membre, done 

AUS 

> - Les «/’ constituent donc à leur 
M =—1 

tour un tenseur sphérique de rang / en vertu de la 

loi du quotient (1). En désignant ces tenseurs par 

une lettre majuscule ronde (5.1) peut s’écrire sous 


la forme : 


de la somme 


ae Ÿ A ; YO, 


1=0 


(5.4) 


le symbole : désignant le produit contracté des 
deux tenseurs défini par : 


+l 
tt) : YU) = >> a Ve 


m=—}| 


(5.5) 


On peut écrire ce développement en termes de 
tenseurs réels, en prenant des combinaisons 
linéaires des Ÿ} qui s'expriment d’ailleurs sous la 
forme de polynômes homogènes de degré L (véri- 
fiant l'équation de Laplace) en Ve, Vey, Ve. Nous 
désignerons donc le tenseur sphérique de rang / 
exprimé sous cette forme par 4 et le dévelop- 
pement s’écrira dès lors : 

CO 
LR = 2 VU: F0, (5.6) 
k 1=0 
où #Ÿ a pour composantes des combinaisons 
linéaires des «*. Dans le cas L = 0, 0 et #0) sont 
des scalaires ; pour ! —1, 9 et #® sont des 
vecteurs ; on a d’ailleurs 1 = V, Dans le 


() Cf., par exemple: A. Licanerowicz, Calcul tensoriel 
(Armand Colin). : 
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RUD — (0(Y, V.) 


LÉ 
sans divergence et de même pour 2% —°f? (le 
symbole <-—> indiquant un tenseur et l’indice (° en 
haut à gauche désignant un vecteur de divergence 
nulle (vecteur solénoïdal). 

Il nous faut donc exprimer le terme général du 
développement (5.6) en fonction du terme général 
de (5.4) et (5.1), en tenant compte des relations 
entre fonctions sphériques. Ce terme général s’écrit 
à partir de (5.1) : 


cas L —2, V® est une dyade : 


Fe 
mn M 
DMC a 


m=—|l 


at) : VU — RA(1) : UN) - 


L 
0 0; m yM — —M 
= a PE D (et di oo PA): 
m = 1 


(5.7) 


En utilisant la définition des YŸ?7, la somme du 
second membre s’écrit : 


12< {(o° + (— 1) «y ”)® COS Mo 


+ of — (— 1) 07”) OÙ sin mp}. (5.8) 


Comme par ailleurs on a : 


m 


RAS 12 NT0 EPR 
—— Ve ® cos mp — (eimp —Æ er im) 
Var VIT 
1 + 
= +), 
1 A @!! 
rs y O}"'sin mo = AC eimp — e— imp) 
27 2iV2T 
1 mr — M 
——(Y : 1) Y 
sl L { ) (4 jh 


(5.8) peut s’écrire encore : 


l 
DS fe tir) (Etes) 


—n) ee — LS Fa) | 


ot — (— 1)" ; ) 


En vertu de (2.3) et de (2.4),les deux premières 
parenthèses de l’accolade sont réelles, et le deu- 
xième terme est réel pour les mêmes raisons. En 
rapprochant (5.9) et (5.7) on voit que les tenseurs 
sphériques réels V® et #0 ont pour composantes : 


_ 


(5.9) 


m ? M yM | 2 
DL) (re, DS mn DR an di 12 y 
2 27 
(5.10) 


FO (xÿ, © ie ET 1) de (— 1)" ao) 


(5.11) 


Considérons maintenant le développement (5.6) 
en nous arrêtant au terme en L — 1 ; le terme 
d'ordre zéro est un produit de scalaires et le terme 
d'ordre Î est un produit scalaire de vecteurs ; nous 
écrirons donc : 


+ — 
le = FO VO + Va Fa = {9 LV, fO (512) 


N°14 


avec 
= FO V0) 


— — : 
VD FD = OV, fo 

Quelles sont les composantes de /(? Elles s’ex- 
priment à partir des of, «l et x; ! en remplaçant 
d'autre part les Y9, Yi et Y;:!en fonction des 
composantes de la vitesse. 

Commençons par f®; par définition, à 4 cor- 


è 1 k 
respond Y$ = ——0@; d'où 
/ 27 
(OE cs O9 — 208 (5.13) 
V9 2Vr 

(d’après lappendice). 

De même pour V(® on a, d’après (5.10) : 

1 1 JE HU 

0, CENT e 2 6 — — LVar7 
He Fi 2 HR 5 /3 ez 
Yi — Ze, Leip — re 3 Y, sin 0 eiv (5.14) 

VOr 2 DC 

1 Free tre Œx : ; 

Y1 ge 17 ei VS Ve sin 0 ete. 
DÉREEN  nt ie 
En formant === 1t"# te vient : 
v) 21 

ARS € Pre ES ei MOVE 

Sn rt . sin 0 cos o — Da S. Ve 
NT En V. LS 3 
a Se 3, sin 0 sin ?=—i 3x Veu. 

On a donc 
nee 1e SG rev 1 
VU) ,g ee Vrai Vez 3 ox ) Ver 


ex.) Fer NA RTE) 


1 Se 
NÉLCE 


Il s'ensuit que les composantes de {® seront par 
définition : 


(1) CRT 1 1 

cu À 1) = —= (a + 4”) 

(1) LRAES ET TA 

fes rt) = (ai ko ) — —VÈ (a —— «}*) 
CDR 3 o 

lez — 2 ls (5.18) 


Dans la seconde partie de ce travail, nous don- 
nerons les équations du premier ordre que véri- 
fient fo et 19°, et nous traiterons le cas du champ 
électrique alternatif. 


Manuscrit reçu le 21 mars 1958. 
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DANS QUELLE MESURE LA MÉCANIQUE CLASSIQUE PEUT-ELLE PRÉDIRE LES TRAJECTOIRES ? (1) 


Par Max BORN. 


Résumé. — Étant donnée l'incertitude sur les conditions initiales — due à l’imprécision inévi- 
table des mesures — les prévisions de la mécanique classique ne se rapportent jamais à une trajec- 
toire précise, mais à un ensemble de trajectoires donné par une distribution de probabilités. La 
tâche réelle de la mécanique est de prédire l’évolution de cette distribution au cours du temps. 

On montre que, dans tous les cas, sauf celui des oscillateurs harmoniques, cette distribution 
s'étale et que la prévisibilité s’estompe. Dans le cas des mouvements périodiques et quasi-pério- 
diques, toute prédiction sur la configuration du système considéré devient impossible à partir 
d’un temps critique calculable en fonction des conditions initiales. Si l’on décrit l’état de ce système 
par des variables angulaires o et leurs variables conjuguées d’action J, on atteint au bout d’un 
temps assez long un état-limite où les J sont presque rigoureusement determinés tandis que les ® 
restent totalement indéterminés. Ceci correspond à l’état stationnaire de la théorie quantique de 
Bohr, mais ici les J ne sont pas nécessairement multiples entiers de h. 

Ces considérations ne doivent nullement effacer la différence fondamentale entre mécaniques 
classique et quantique ; leur but unique est de montrer que l’idée de détermination en mécanique 
classique repose sur une hypothèse dénuée de sens et que l’absence de détermination en mécanique 
quantique ne peut servir d’objection à cette théorie. 


Abstract. — In view of the uncertainties concerning the initial conditions, due to inevitable 
experimental imprecision, the predictions of classical mechanics never describe a unique trajectory 
but rather a group of trajectories defined by a probability distribution. The true task of mecha- 
nics is to predict the evolution of this distribution with time. It is shown that for all cases, except 
for the harmonic oscillator, this distribution flattens out. 

In the case of periodic and quasi-periodic systems, any predictions concerning the configuration 
of the system considered become impossible beyond a critical value of the time. This critical 
value is easy to calculate as a function of the initial uncertainty. If the state of the system is 
described by angular variables ® and their conjugated action variables J, then, after a certain 
time a stationary limiting state is reached for which the J’s are nearly exactly determined whilst 
the os are totally in determined This corresponds to the stationary state of the Bohr quantum 
theory, but here the J’s may take any values and not necessarily integral values of h. 

These considerations do not affect the fundamental differences between classical and quantum 


mechanics. 


Their intention is solely to show that the idea of determination in classical mechanics 


is based on a hypothesis which is void of meaning and that the absence of determination in quan- 
tum mechanics cannot serve as an objection to this theory. 


Introduction. — La question de l’influence d’une 
petite indétermination dans les conditions initiales 
sur le calcul de l’évolution future des systèmes de la 
mécanique classique a été examinée dans plusieurs 
travaux (?). Je reviens sur ce problème, d’abord 
parce que les raisonnements que nous avons publiés 
D. J. Hooton et moi, sont assez compliqués et ne 
laissent pas apparaître assez clairement la racine 
intuitive du résultat principal. 

Mais la raison essentielle d’un nouvel exposé est 
d'ordre principiel, philosophique. Il existe encore 
des groupes de physiciens et de philosophes qui 


(2) Manuscrit remis simultanément au J. Physique Rad., 
pour lequel M. Edmond Bauer a bien voulu accepter dele 
traduire, et aux Z. Phystk, où il vient de paraître en 1958, 
153, 372. 

(2) Born (M.), Acad. Copenhague Math. Phys., 1955, 
30, n° 2. Born (M.) et HooTen (D. J.), Statistische Dyna- 
mik mehrfach periodischer systeme. Z. Physik, 1955, 142, 
201. Statistical dynamics of multiply periodic systems. 
Proc. Cambridge Phil. Soc., 1956, 52, 287. Le P' W. A. Fock 
de l’U. R. S. $. m’a appris que son élève N.S. KryLow a 
publié des considérations analogues dans le but de fonder 
sur de nouvelles bases la mécanique statistique. 


mettent en doute l’interprétation statistique de la 
mécanique quantique et les idées de complémen- 
tarité de Bohr. 

Je ne vais pas ici énumérer systématiquement 
toutes ces tendances mais seulement en citer deux 
exemples appartenant à des camps tout à fait 
opposés. 

Le premier exemple est celui du philosophe 
catholique Van Laer (*). 

Celui-ci ne distingue pas seulement causalité et 
déterminisme (point sur lequel beaucoup de phy- 
siciens — dont je suis — seront d’accord, bien que 
pour de tout autres raisons), mais aussi détermi- 
nisme et prévisibilité (ou possibilité de prédire). 
D’après lui, le déterminisme est l’expression de la 
nécessité causale dans l’abstrait. Il est une condi- 
tion nécessaire, mais non suflisante de la prévi- 
sibilté, car, pour prédire à l’avance, il faut une 


(3) Van Laer (P. Henry), Philosophico-scientific pro- 
blems (Pittsburg Pa, The Duquesne Univ. Press, 1953, 
pp. 115-132). Voir aussi Koren (Henry J.), Readings in the 
Dore of Nature. (The Newman Press, Westiminter 
Maryland, pp. 388-401.) 
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connaissance suffisante de l’état initial. Van Laer 
reproche aux physiciens de méianger ces deux 
« concepts. Pour moi, cette distinction me semble une 
coupure de cheveux en quatre, en contradiction 
d’ailleurs avec un principe heuristique qui est cons- 
tamment appliqué avec succès en physique 
moderne : les concepts auxquels ne correspond 
aucune situation de fait, au moins en principe, 
n’ont pas de place en physique. Le déterminisme, 
au sens de Van Laer, est un de ces concepts super- 
flus, car il ne peut être vérifié que par une pré- 
diction. D'ailleurs Van Laer semble croire que sa 
distinction permettrait de sauver le déterminisme 
en mécanique quantique, car celle-ci ne nierait que 
la possibilité de prévoir et non la détermination en 
principe. Mais cela est une erreur de fait, car la 
mécanique quantique ne permet, d’un point de vue 
purement conceptuel, la fixation rigoureuse, ni de 
l’état initial ni d’aucun autre état ({). 

C'est pour éviter toute discussion oiseuse de ter- 
minologie que j’ai évité dans le titre du présent 
travail le mot de déterminisme. 

Le deuxième exemple de l’opposition à l’inter- 
prétation statistique de la mécanique quantique 
est un petit volume paru récemment en français (°) 
qui contient 9 articles de chercheurs à tendance 
marxistes. Leurs arguments principaux sont de 
deux ordres : d’une part ils s'appuient sur des auto- 
rités telles que Planck, Einstein, Langevin, De 
Broglie, Schrôdinger, etc. qui, ayant eu une part 
active importante à la construction de la méca- 
nique quantique, ont néanmoins rejeté l’inter- 
prétation statistique de la théorie ; d’autre part, ils 
se fondent sur la philosophie marxiste dont un des 
dogmes est le déterminisme des phénomènes natu- 
rels. 

La beauté et la cohérence logique de la physique 
classique, la simplicité de la philosophie sous- 
Jacente possèdent un pouvoir d’attraction puissant 
auquel ne peuvent se soustraire même certains 
grands esprits révolutionnaires. Planck a expliqué 
dans une esquisse d’autobiographie combien le 
fascinait l’idée de lois rigoureuses de la nature, et 
l’on trouve chez Einstein des déclarations ana- 
logues. Je crois cependant que ces convictions 
reposent sur une illusion provenant de l’identi- 
fication, courante au siècle dernier, du monde phy- 
sique et du continu mathématique. En effet, en 
mécanique classique, un mouvement n’est déter- 


(*) Max PLanckx voulait sauver le déterminisme en 
appliquant au comportement de la fonction d’onde qui 
est effectivement déterminée par les conditions initiales et 
aux limites. Mais comme le carré de la fonction d’onde ne 
représente qu’une probabilité, rien n’est gagné ainsi pour 
la mécanique du point. Voir par exemple PLancr (M.), 
€ La Causalité dans la nature » dans le livre: Weg zur 
Physikalischen Erkenntniss (S. Hirzel, Leipzig 1933), p. 233. 

(5) Recherches internationales à la lumière du Marxisme. 
Physique : Quelques problèmes philosophiques, 1957, n° 4 
(septembre-octobre). 
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miné pour tous les temps par les conditions initiales 
que si celles-ci sont données avec une précision 
absolue. Les mathématiciens ont créé le concept 
raffiné de continu dans un but purement logique, 
à savoir: établir sans contradiction l’analyse, la 
théorie des fonctions et la géométrie analytique. 
Pour cela il est indispensable de définir le concept . 
de nombre réel par des passages à la limite, de telle : 
sorte qu’il possède une exactitude absolue et que, . 
par exemple, en représentation décimale, deux. 
nombres sont considérés comme différents si, leurs » 
premières décimales coïncidant, les (nr + 1) %% 
diffèrent — même si nr — 10%. En physique, un 
tel concept n’a aucun sens. Il y a des années déjà 
que Felix Klein a demandé que l’on créat une 
mathématique des approximations adaptée à la 
physique. Plus tard, c’est la physique elle-même 
qui a pris ce problème en mains et l’a résolu dans : 
une certaine mesure en introduisant le concept de 
probabilité et la statistique. 

Le déterminisme mécanique, comme on l’entend 
d'habitude n’a de sens que dans l’espace abstrait 
mathématique, où la position d’un point peut être 
donnée avec une précision absolue par des nombres 
réels, ainsi que la vitesse et d’autres grandeurs. 
Dans l’espace réel de la physique le déterminisme 
n’a pas de sens. Car il postule que l’on peut consi- 
dérer comme distinctes des choses qu'aucun appa- 
reil de mesure, si raffiné soit-il, ne saurait dis- 
tinguer. De tels concepts n’ont, d’après notre prin- 
cipe heuristique, aucune place en physique. 

On peut cependant donner du déterminisme une 
définition qui a un sens physique. Tout état phy- 
sique implique une imprécision petite, mais finie. 
Il ne sera donc jamais représenté par un nombre, 
mais par une distribution de probabilités et la tâche 
de la mécanique est de prédire, à partir d’une dis- 
tribution initiale, les distributions futures, On peut 
ensuite faire tendre vers zéro l’incertitude initiale. 
Si alors le domaine d’incertitude qui, par suite du 
mouvement, résulte au temps { du domaine initial, 
tend vers zéro dans toutes ses dimensions (quand le 
système est dynamiquenrent stable), l’état ultérieur 
est prédictible et la’ théorie peut être considérée 
comme déterministe. 

Cette définition du déterminisme se distingue de 
la définition habituelle par l’ordre de succession des 
opérations dans un double passage à la limite. La 
méthode habituelle qui consiste à rassembler 
d’abord en un point la dispersion initiale, à calculer 
ensuite la trajectoire et la dispersion finale (en un 
temps quelconque, en particulier au temps t), n’a 
pas de sens physique et doit être remplacée par la 
règle suivante : Calculer d’abord les trajectoires et 
la dispersion finale en un temps { quelconque à 
partir d’une distribution initiale donnée, puis res- 
serrer la dispersion initiale et déterminer comment 
se comporte la dispersion au temps £. Le problème 
est alors de savoir s’il existe des systèmes stables 
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en ce sens, C'est-à-dire tels que l’on puisse calculer 
rigoureusement à l’avance leur évolution. 

Dans la suite, nous montrerons que la mécanique 
classique n’a pas cette propriété, sauf dans un cas, 
celui des systèmes d’oscillateurs harmoniques 
idéaux qui, par suite de l'absence d'interactions, ne 
sont guère utilisables. | 

En outre, nous ramènerons le cas du système le 
plus général au plus simple, celui du mouvement 
Hinéaire périodique d’un point matériel. Il semble 
donc peu vraisemblable que l’on puisse imaginer 
d’autres lois de mouvement qui se comportent de 
manière déterminée. 

Mais alors est retirée toute base aux efforts qui 
tentent de fonder la mécanique quantique sur des 
h ypothèses de coordonnées cachées obéissant à des 
lois déterministes (comme le propose par exemple 
J. P. Vigier (f). 

Bien que j'aie grandi moi-même dans les idées du 
siècle dernier, il m’a toujours semblé étrange qu'il 
dût exister un domaine, celui de la mécanique et de 
la physique classique construite sur son modèle, 
où tout soit absolument exact et libéré des incer- 
titudes qui, par ailleurs, dominent la vie et la 
pensée humaines. C’est pourquoi je pense que l’éli- 
mination de cette rigueur fictive par la physique 
moderne est un progrès dans la tendance vers une 
conception unitaire du monde. 

Il faudrait déjà tenir compte dans l’enseignement 
de la thèse que nous soutenons ici. À mon avis, des 
considérations analogues à celles qui vont suivre, 
notamment celles du paragraphe 1, pourraient 
trouver leur place dans des cours élémentaires de 
mécanique, afin d'éviter dès l’abord la formation du 
préjugé déterministe. 


1. Mouvement apériodique d’un point matériel 
non soumis à des forces. — Le cas le plus général 
peut se ramener au plus simple, celui du mou- 
vement libre apériodique d’un point matériel. Nous 
allons traiter celui-ci d’abord de manière tout élé- 
mentaire et intuitive, puis analytiquement, du 
point de vue de la dynamique statistique. 

La loi du mouvement est la même pour toute 
trajectoire à courbure quelconque sur laquelle se 
meut sans frottement le point matériel. On peut 
donc se borner à considérer une trajectoire recti- 
ligne, comme dans le mouvement libre. 

Soient s la distance à un point fixe O, v la vitesse. 
On a alors pour le point matériel considéré : 


S = So + Ut D = (4) 


Considérons un ensemble de mouvements défini 
par une densité de probabilité (D. P.), celle de l’état 
initial P(s0, Vo: 0) — (80, 29). Comme le nombre des 


{s) Voir par exemple l’article de J. P. Vicrer à la 
page 130 du livre cité p. 2, et le volume : Structure des 
. micro-objets dans l'interprétation causale de la théorie des 
quanta, Paris, Gauthier-Villars, 1956. 
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particules est invariable, la D. P. au temps { sera : 
PDT DEN Dr 0). (2) 
Nous supposerons en particulier dans la suite 


que la distribution initiale est la distribution nor- 
male de Gauss : 


exp [— : Q(s, v)] (3) 


ee 


f est normé à l’unité, $ et v sont les valeurs 
moyennes initiales de s et v, 6 et + sont leurs dis- 
persions (écarts quadratiques moyens). 

En fonction des variables réduites 


BE S "ND 


où 


ET 


FT EUIO ES CN NUE (5) 
la D. P. prend la forme 
Jr, 0 NCAA UE D exp F3 Gate, 4200)| 0) 
Si l’on pose : 
=) fo (7) 


Q, peut s’écrire sous les trois formes équivalentes : 
l (x — Miy)? + y? 
2? — 2Mixy + (1 + M?) y? 


Mit 2 (8) 
(1 + M?) er F as) Si 


+ 


Qilx, y; t) ee 


x? 
Décrivons d’abord de manière intuitive et quali- 
tative la propagation de la D. P. Dans le plan sv, 
(1) représente un écoulement laminaire parallèle à 
l'axe s, où chaque domaine subit une déformation 


IR TIR 


de glissement (fig. 1) : un carré devient un paral- 
lélogramme dont les côtés parallèles à l’axe des x 
(ou des s) restent invariables (AB — A’B”, mais 
BC  B'C') ; un cercle se change en ellipse dont 
les axes diffèrent peu des diagonales du parallé- 
logramme. Quant { augmente, le parallélogramme 
s’étire évidemment en devenant plus étroit, et il en 
est de même de l’ellipse. Les aires restent cons- 
tantes, conformément au théorème de Liouville de 
la mécanique. Quant t devient très grand, la direc- 
tion du grand axe Ë de l’ellipse se rapproche de la 
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direction Ox, celle du petit axe n de la direction Oy. 
Le domaine de probabilité notable s’étend de plus 
en plus, pour { croissant, dans la direction Ë, tandis 
qu’il se rétrécit dans la direction 1. Mais on ne peut 
dire d'emblée qu’il en est de même pour les direc- 
tions Ox et Oy, bien que les deux systèmes d’axes 
se rapprochent de plus en plus. Les choses sont un 
peu compliquées et ne peuvent s’éclairer que par 
un traitement analytique. 

Tout d’abord, de la première et de la troisième 
expression de @;, l’on déduit immédiatement que 
la D. P., en fonction de v (ou y), est indépendante 
de la’ valeur de s (ou x) : 

(9a) 


Poly, t) = 1 PACAURONIAUE Ver ÉnaUel 


et qu’inversement la D. P., en fonction de s (ou x) 
est indépendante de v (ou y) : 


+00 
Pix, t) Le 


1 me 
Var + M? ae | 21 + M? x) 


Les deux distributions sont gaussiennes, mais, 
tandis que celle en y est invariable, celle en x est 
fonction du temps, de 1 + M? 1°. Il en résulte 
que la dispersion de s s’écrit : 


ot) = os VA _. M? 2 = Vo? + ee 


P,(x, y. t) dy 


(98) 


(10) 
Soit : 
x = Ë cos y — nsin y 
— Ésin y + n cos y 
la transformation aux axes principaux de l’ellipse 
Q,(x, y,t) — constante. 
Un calcul élémentaire montre que l’angle y entre 
l’axe Ë et celui des x est donné par : 


(11) 


tg 2x — 1/2 Mt (12) 
(11) transforme Q,(x, y, t) en 
QE, n, t) = (ËJa)? + (n1b)? (13) 


où a?et b? sont les inverses des racines de l’équation 
séculaire de Q,, que l’on peut écrire en tenant 
compte de (12) et de l'identité 


1 +2 ctg? 2% = (1/2) (te? x + cte? y) 
22 — à (te? x + ctg? y) + 1 = 0, (14) 


équation du second degré réciproque dont le 
racines sont À, —tg? y, À, —ctg? y. On a donc : 


MIE cie x, 


is CAO (15) 
b—1/VXR = te y, 
À Q, correspond la D. P. 
PAE, 0 = exp QE mal (16) 


La constance du produit a.b est l'expression du 
théorème de Liouville : Paire de l’ellipse Q, = r? 
(r — const.) est indépendant du temps. Si le 
produit a.b —1, c'est parce que, d’après (8), 
l'équation de l’ellipse devient, pour { = 0, celle du 
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cercle x? + y? —r?, et que l’on a par conséquent 
n(ra){rb) = rt soi Gb 
Le carré de la dispersion en Ë& s'écrit : 


E se PA(E, y, 0) dE dy — a. (17) 

On obtient de même pour la dispersion en net le 
coefficient de corrélation entre ëet n 

m = b? (18) 


D'où l’on peut tirer à l’aide de (11) les grandeurs 
correspondantes pour «et y 


En = 0. 


22 1=%a2 008 F7 DE SIRET 
y? — a? sin? y + b? cos x (19) 
ay — (a? — b?) cos y sin x. 

Si l’on introduit ici les expressions (15) de &et b, 
on obtient par un calcul élémentaire 


22 —1 + 4 ctg29y = 1 + M?12 = [o(t)/o]? 


y? =1 (20) 


my = 2 ctg 2% — Mi. 

Les deux premières valeurs coïncident avec ce 
que donne le calcul direct (formules (94) et (9b)). 
La troisième montre que les distributions en x et y 
ne sont pas indépendantes et donne une mesure de 
leur corrélation : celle-ci disparaît, bien entendu, 
pour # > 0, et tend vers l'infini avec t. De (11) et 
(12) on déduit que 

(21) 
4 > rla,z—(1/2)(ê— 1), y —(1/2)(6 + n) 
pour t— co: y 0, x—£# JET 

Les axes principaux &, »n sont donc d’abord — 
comme on le voit sur la figure 1 — les diagonales 
du carré OABC construit sur les axes des x et des y, 
puis, quand t croît, ils se rapprochent de ces axes. 

Néanmoins, les composantes x, y et 6,1 du mou- 
vement se comportent de manière très différente : si 
la dispersion en x, et celle en &, tendent toutes deux 
vers l’infini avec £, celle en y reste constante, tandis 
que celle en n, inversement de celle en Ë, tend vers 
zéro. La raison en est évidemment que les compo- 
santes £, n sont indépendantes et-non les compo- 
santes %, y. 

Nous allons rattacher à ce point une remarque 
qui anticipe sur les résultats des paragraphes sui- 
vants. [Il sera démontré que les mêmes raison- 
nements exactement restent valables pour des sys- 
tèmes périodiques, non seulement à un degré de 
liberté mais à un nombre quelconque de ceux-ci. La 
coordonnée s et la vitesse v sont alors remplacées 
par les paires de variables angulaires et d’action +, 
J dont s’est servi à Bohr dans sa théorie quasi- 
classique des quanta. Si alors la dispersion d’une 
variable angulaire devient grande par rapport à sa 
période, l’état, en tant que défini par cette variable, 
devient totalement indéterminé. | 

Si l’on se rappelle maintenant la théorie de 
Bohr, on voit que son trait essentiel est l’introduc- 


pour t— 0 : 
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tion d’états stationnaires distincts définis par des 
valeurs discrètes des variables d'action (J = nA, 

— 1, 2,...). Les variables angulaires ne jouent 
aucun rôle, elles restent indéterminées. 

Mais cela est exactement le comportement de 
l’ensemble (gaussien) classique de trajectoires pour 4 
grand quand on le décrit par les variables &, 1 : 
l’état en £ devient totalement indéterminé, l’état 
en » rigoureusement défini. En outre, on ne peut 
plus distinguer &, n de x, y, x étant l’écart des 
variables angulaires à partir de leur valeur 
moyenne, y le même écart pour les variables 
d’action. 

Ce résultat éclaire la méthode de quantification 
de Bohr d’un point de vue nouveau. Elle ne se 
rapporte évidemment pas à des trajectoires uniques 
mais à un ensemble de trajectoires dans l’état final 
stationnaire { — co, 

Bien entendu, on n’explique nullement ainsi la 
condition que seules sont permises les valeurs de J 
multiples entiers de h. Ceci est du ressort exclusif 
de la mécanique quantique. 


2. Mouvement périodique libre de toute force. — 
Soit une particulé se mouvant, libre de toute force 
(sans frottement), sur une ligne fermée de lon- 
gueur /. N’ont alors de sens que les fonctions pério- 
diques de s ayant la période L. Or la distribution 
gaussienne (3) n’est pas périodique. [] faut la rem- 
placer par une fonction périodique analogue Le 
plus simple est d’écrire 


HISD) MD UT IS EL" D) 
È k=—0co 
f étant la fonction définie par (3) 
La D. P. au temps test alors 


(22) 


+00 
DRE SEE PAIN D) (25) 
k=—00 
Il est facile de voir que, si f(s, v) est normalisé 
dans le domaine —co <s< +oo,—co<v< +00 
F(s, v) et par conséquent P(s, v, t) seront norma- 
lisés dans le domaine 0 £<s < {, —oo < v +00. 
La démonstration est fondée sur le théorème 
suivant : 
Si G(s) est périodique 


P(s, v,t) = F(s — nt, v) — 


eo (24) 
on a 
Je Ge Al). ds == ae &(s) f{s) ds (25) 
En effet 
+00 
Je GA Fisids => Ë G{s) fls + A1) ds 
A Cure los 
k=—00 Vkl 
+ co (R+1)1 
— 2. JE G{s) f(s) ds 


et la sommation donne (25). 
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En particulier, pour G —1,ona 


ile Fs) ds sl f{s)- ds 


ce qui démontre notre proposition sur la norma- 
lisation de F. 

Il s’en suit en outre que la valeur probable de 
toute grandeur physique, c’est-à-dire périodique 
en s, G(s, v, t) calculée sur la période / à l’aide 
de P(s, v,t) — F(s — ut, v), est identique à la 
valeur probable de G calculée sur tout le domaine 
— COS L + coavec f(s — vé, v). 

Par conséquent on a le droit d’étendre les résul- 
tats obtenus dans le paragraphe précédent au cas 
périodique. Il suffit simplement de considérer 
comme identiques tous les points dont les valeurs 
de s coïncident modulo ! (c’est-à-dire dont la diffé- 
rence est un multiple entier de l). 

On admettra naturellement pour la dispersion 
initiale en s 


(26) 


Cl: (27) 


La période en æest alors //6 > 1. 

Si alors o(t) > {, aucune prédiction quant à la 
position ne sera plus possible. C’est ce qui arrivera, 
d’après (10), quand ot) = 6 Mt = rt > 1, c’est-à- 
dire quand 

DUC RON 


ONE (28) 


S1 la précision sur la vitesse est grande, si + est 
petit, l'intervalle critique de temps & est grand. 
La figure 2 illustre géométriquement ce processus. 


Vi 


0 Lo 


FrG:22 


Le parallélogramme de la figure 1 y est (pour # 
grand) étiré en un ruban long et mince qui s’étend 
sur un grand nombre de périodes //6 ; les sections 
de ce ruban qui se trouvent dans les périodes supé- 
rieures sont ramenées dans la première. Il en est de 
même exactement de l’ellipse inscrite dans le paral- 
lélogramme. Le cercle initial engendre ainsi une 
série de bandes parallèles, qui recouvrent unifor- 
mément, quand t est assez grand, tout le domaine 
de la première période. 

Ce cas contient en particulier celui du rotateur 
libre, si l’on remplace respectivement s et v, par 
l'angle de rotation 6 et la vitesse angulaire © ; la 
période est alors 27. 

Il en est de même du cas traité autrefois (7) d’une 


(7) Voir la note (?) de la page 1. 
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particule qui oscille sur une droite entre deux parois 
où elle rebondit élastiquement ens —0ets —1. I 
suffit, pour étendre à ce cas les considérations pré- 
cédentes de remplacer la fonction Æ(s, v) de l’équa- 
tion (22) par 
- F(s, dv) (29) 


= Due f(QRI + æ — ot, v) + H(2KI — x — v1, — v) } 
et de remplacer le mouvement oscillant par deux 
mouvements de vitesses égales et opposées s’effec- 
tuant sur des droites illimitées : il faut alors 
ramener tous les points extérieurs à l’intervalle 
0<s< 1 à l’intérieur de ce domaine par réduc- 
tion modulo 2/. 


3. Mouvement général à un degré de liberté. — 
Il est possible de ramener aux exemples traités 
dans les paragraphes précédents le cas général d’un 
mouvement à un degré de liberté sous l’action de 
forces conservatrices. 

Soient g, p la coordonnée et le moment, Æ(g, p) la 
fonction de Hamilton. Le problème est résolu si 
l’on peut, à l’aide d’une transformation canonique 
(g, p) (+, J), introduire de nouvelles variables 
telles que A(q, p) — E(J) indépendant de +. Si le 
mouvement reste dans un domaine fini, il est 
périodique et l’on a 


L 


= Ÿ pdq 


(30) 


intégrale étendue à une période. On donne alors à J 
le nom de variable d’action. à © celui de variable 
angulaire. On norme généralement © de manière 
que la période soit 1. Le mouvement est représenté 


par 


où 


PPT OM ER NES (31) 
dE 
AY 47 
est une fonction de J, sauf dans le cas de l’oscil- 
lateur linéaire de fréquence v,, où £ = v, J. 
Abstraction faite de ce cas exceptionnel, nous 
pouvons étendre nos résultats antérieurs au cas 
général, car l'intervalle en J dans lequel la D. P. 
P(®, J, t) présente des valeurs notables est petit. 
Considérons les développements 


(32) 


E(J) = E + 5{J 


(33) 


VAR En PE AS 
ne tenons compte que des termes écrits : nous 
sommes ramenés exactement au cas où il n’y a pas 
de forces ; les résultats que nous avons établis 
montrent que la dispersion en J (correspondant à v) 
ne croit pas. Il en résulte que cette approximation 
est valable pour tous les £. 
En particulier on a pour la D. P. 


P(9, J,t) = (1/2xo+) exp ee Q(o, J,1t)] (34) 


4 


Oo Tue (ES) Be ie (33) 


(o} 


où 


Si l’on pose 


x = (®— © — Vt)o, y = (J —J)/r (36) 
les formules (3) restent valables pour Q;, avec 
M — (fo) — (rio) (4? EldJ?);. GA 


L'époque critique où toute prévision devient 
impossible a maintenant pour valeur 

te INTRO EU (38) 

Toutes les autres considérations restent sans 

changement. Dans le cas de l’oscillateur linéaire et 

dans celui-là seul on a, d’après (37) M =0,u —0 

et, par suite, & — co: la distribution initiale se 

conserve. 


4. Système à n degrés de liberté. — Il reste à 
examiner le cas d’un système conservatif à n degrés 
de liberté. Soient 4, px (k — 1, 2 ... n) les coor- 
données généralisées et les moments, Æ/(q, p) la 
fonction de Hamilton. On supposera que l’on peut 
introduire de nouvelles variables canoniques ®x, Jy 
pour lesquelles 7 se transforme en une fonction 
des J seuls 

H{g.p) — E(J). (39) 

Les mouvements pour lesquels les g, p sont 
fonctions périodiques des ® sont dits multipério- 
diques. Les w et les J sont alors appelés variables 
angulaires et d’action. Comme d’habitude les 
périodes en les sont prises égales à 1. 

Pour les mêmes raisons qu’au paragraphe pré- 
cédent il suffit de ne tenir compte que des termes 
écrits dans les développements 


E(J) =E + 2 (Jr — Jr) 
+ 1/2 > ue(Jr 5. Jx) (7 Re Ji) A 0) 
wlJ) = + D ef = Ji) + 2. 


Les équations canoniques s’écrivent, pour les +, 


AE 
do ea dJx us 
+ VE, 1e 0. (41). 
Elles ont les solutions 
= px + vw, Te = J9; (42) 


Dans l’espace g, p, la D. P. P(q, p, t) satisfait à 
Péquation différentielle de Louville 


dP DES ss 
dd 

dt k \0x px  dPr dgr (3 

qui, dans l’espace +, J prend la forme 

2P dP 
Me en Fi 
Sa solution satisfaisant à la condition initiale 

P(?, J, 0) = f(p, J) (45) 


s'écrit 


P(9, dt) = flqi— vit, ga vat,.… ; Ji, ae. (6): 


NOT 


Choisissons pour f une fonction de Gauss normée 
(6, J) = HIT)" 0, 


On T1 - - Tr] 

exp [5 Qle, 1 (47) 
où 
Qe, 2) = DE Lee — ou) lou + [Ur — Hi) } (8) 


x, Gx Sont les valeurs moyennes initiales et les dis- 
persions des x ; J4et 7 celles des J4. 
Si l’on introduit les variables réduites sans dimen- 
sions 
(49) 
1 NU 1Tx) (Jy — Jr). 


La solution (46) prend ia forme 


2x = (1/0x) (9e — 9x — Vi), 


Pile, y, 0) = (1/2m)% exp [— 1/2 Qi(x, y, 1] (50) 
où 

Que, v,9 = Zf(m—1Z May +uè} (51) 

Mu = Mu = pui (sxlo) + (rlo) À. (52) 


Servons-nous maintenant de la notation matri- 
cielle : considérons les x, et y;comme éléments des 
matrices à une colonne x et y ; x et y seront leurs 
matrices transposées à une ligne. M est la matrice 


symétrique dont les éléments sont Mat M}, 
Nous avons alors : 


Qu = (à — tÿM) (x — My) + ÿy 
= àx — 2iMy + y(1 + 1 M?) y, 


(53) 


où il faut noter que * Myet y Mxsont deux expres- 
sions de la forme bilinéaire À My Ty Y. 


ke! , 
Faisons subir aux 4 et y, la même transformation 
orthogonale 
z = Ur’ y = Uy (54) 
U étant la matrice orthogonale qui diagonalise la 
matrice symétrique M 


ÜU = UU =1 (55) 
m; 0 0 
gmMu=m=| ° me... (56) 


U diagonalisera de même toute puissance de M, 
donc tout polynôme en M, avec des valeurs propres 
qui sont représentées par le polynôme corres- 
pondant en m. 


Üf(M) U = f{m). (57) : 
Par la transformation (54) Q, se change en : 
Qx = za — 21% my" + Y'A + m°) y’, (58) 


soit développé : 


On = 24/3 — 2imja y + (1 + Emity 5} (59) 
î 
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On peut transcrire l’expression entre accolades 
sous l’une des deux autres formes données en (8). 

La DP. (50) est ainsi décomposée en un produit 
de facteurs dont chacun a la forme admise dans les 
paragraphes précédents. Toutes les conclusions 
antérieures restent valables ; il faut seulement rem- 
placer les variables réduites par des combinaisons 
linéaires bien définies, correspondant à une même 
rotation dans les es, aces des x et des y. 

Quand £ > 1/m,, on peut remplacer les disper- 
sions en Z’ et y’ par les premiers termes de leurs 
développements en mn; t. 

Les périodes du systèmes sont égales à 1 pour 
les @, à 1/04 pour les z3. Si l’on pose 6? =Y of, 

k 


on peut considérer 6 comme l'extension du sys- 
tème ; et comme D'aË —=Y 1, © est aussi l’exten- 
k 


k 
sion dans l’espace x’. 
Si m est le plus petit des nombres m1, m, ... my 
et si l’on pose 


te = 1/mo, (60) 


on est certain que, pour t > &,, l’on a atteint l’état 
limite où tous les J, sont, pratiquement, égaux à 
leurs valeurs moyennes J,, tandis que les @, restent 
tout à fait indéterminés (au sens expliqué au $ 1). 

(60) est la généralisation rationnelle de (38) au 
cas de nombreux degrés de liberté. 

L'état limite correspond exactement à l’état 
stationnaire de la théorie quantique de Bohr qui 
est décrit par des valeurs précises des J (multiples 
entiers de À) sans que les valeurs des @ y inter- 
viennent. 

Le seul cas d'exception est celui d’un système 
d’oscillateurs. M et tous les m; s’annulent alors ; 
tk est infini; une D. P. initiale quelconque se 
conserve. Toutes les théories pures de champs sans 
couplage sont équivalentes à des systèmes d’oscil- 
lateurs. Elles appartiennent donc à ce cas excep- 
tionnel. Mais, si l’on introduit des termes de cou- 
plage (non linéaires) entre les champs de différentes 
espèces, le phénomène que nous venons d’étudier 
doit apparaître : les variables angulaires s’estom- 
pent, les variables d’action se précisent. Il serait 
intéressant d'examiner ce sujet de plus près. 


5. Exemple : le pendule. — Afin d'évaluer au 
moins sur un exemple concret la grandeur du temps 
critique, nous considérerons le pendule pesant, 
avec la fonction de Hamilton (masse 77, bras de 
levier du centre de gravité /, moment d'inertie À) : 


H = pê[2A + mglli — cos 0) = p6[2A 
+ (1/2) mgl(0? — (1/12) 04 +...) (61) 


C’est un oscillateur anharmonique. Ecrivons la 
fonction de Hamilton sous la forme normée 


H = 1/2[p? + (2rvo)? q°] + ag° + bgé +... (62) 
in 
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(62) découle de (64) par la substitution 
g—=VA0 p = pelVA 
Vo = (1/27) Vmgl]A b — — 7° 64. 


(63) 


GR: ==) 


Le calcul de À en fonction des variables d’ac- 
tion ($) donne en toute généralité 


H— E = vo, J + (1/2) u J? + .… (64) 


où 


u = (15/2) a2](27)6 v8 + 36[(2x)t ve. (65) 


Si l’on tient compte des valeurs (63) de a et b, 
on obtient 


u = —1/(32r°À). (66) 


Soit 7, — "5! la période d’oscillation harmo- 
nique, On à : 


= No IE 2 Er ANT IE (67) 


Pour rendre intuitif le sens de cette formule, 
considérons le cas où le pendule est écarté de sa 
position d'équilibre jusqu’à son élongation maxi- 
mum 0, et où celle-ci est donnée avec une dis- 
persion relative y — d0,,/0%. 


Pour te 0 = 0e VA 0, 


E = H = (1/2) (2rv)° qm — (1/2)(2xv0)? A 0%. (68) 
Or on a 
(69) 
nn LS dE 
m VU — J)? = I — on = (2x)? À Vo On 0m 
et d’après (67) 
telTo — 80m Y- (70) 


Si l’élongation initiale est, par exemple, 1/20 
détour 00 = 710 4/5 04119 /#et#hon 
obtient &/T, = 72]y. Si la précision en 6, est, 
par exemple, de 1/1 000, &/T, — 72 000 et, pour 
T, —1 sec, & — 20 heures. Lorsque l’amplitude 
initiale diminue, le temps critique croît en raison 
inverse de son carré. 

Pratiquement cette limite de détermination 
n'entre pas en ligne de compte, car d’autres in- 
fluences (frottement, résistance de l'air, etc...) 
jouent bien plus tôt. Le cas des horloges est diffé- 
rent du précédent, car leur pendule n’oscille pas 
librement, mais avec apport d'énergie extérieure. 


6. Conclusion. -— Bien entendu, les considé- 
rations qui précèdent ne doivent nullement effacer 
la différence fondamentale entre mécaniques clas- 
sique et quantique. Leur but unique est de montrer 


($) C£., par ex. : Born (M.), Atomtheorie des festen Zus- 
tandes (Teubner, 1923), pp. 671-674. 


NU! 


que lidée de détermination en mécanique clas- 
sique repose sur une hypothèse dénuée de sens et 
que l’absence de détermination en mécanique quan- 
tique ne peut servir d’objection à cette théorie. 
Bien plus, il est également raisonnable de formuler 
dès l’abord la mécanique classique de manière 
statistique. On obtient ainsi un avantage supplé- 
mentaire : certaines objections à la mécanique 
quantique ne peuvent plus être soulevées, celles qui 
se rapportent à l’impossibilité de séparer les points 
de vue de l’objet et du sujet. 

En effet si, en mécanique classique aussi, toute 
affirmation est de probabilité, elle représente, 
comme en théorie quantique, une combinaison de 
données objectives (l’état initial) et de traits sub- 
jectifs (l'attente de l’état final). Et c’est cela surtout 
ce qu'on appelle la « réduction de la fonction 
d’ondes » par l’observation : on calcule par la 
mécanique quantique la propagation d’une fonc- 
tion d'ondes et trouve ainsi la distribution de proba- 
bilité d’un phénomène atomique. Aussitôt que ce 
phénomène est observé (par exemple l’apparition 
d’un électron en un lieu déterminé), toute la fonc- 
tion d’ondes devient sans valeur — elle s’évanouit 
pour ainsi dire ; et il faut en calculer une nouvelle 
sur les nouvelles données. 

Einstein, Schrôdinger et d’autres ont interprété 
cela comme impliquant pour l’acte d'observation, 
une « action à distance à vitesse de propagation 
infinie » et l’ont rejeté comme dénué de sens. La 
même chose exactement apparaît déjà en méca- 
nique classique, dès qu’on la formule statisti- 
quement et que l’on considère la D. P. comme étant 
la grandeur à déterminer. Alors, comme en méca- 
nique quantique, le but de la théorie n’est pas de 
décrire une actualité, mais une potentialité (dans la 
terminologie d’Aristote), et tout acte d'observation 
change la situation de fond en comble. Ce qui fait 
la véritable originalité de la mécanique quantique 
ce n’est pas ce qui précède, mais le fait que la 
variable primaire n’est pas la D. P. elle-même, 
mais une fonction d’ondes, dont la carré donne la 
D. P. De là viennent d’abord les relations d’incer- 
titude de Heisenberg qui limitent les dispersions de 
couples de variables conjuguées, ensuite les effets 
d'interférences par superposition des champs 
d’ondes. Ce sont ces effets qui sont déterminés par 
la constante de Planck. Leur existence entraîne 
nécessairement l'interprétation statistique. En 
théorie classique au contraire celle-ci dépend du 
rejet ou de l'acceptation de l'hypothèse que le 
concept de précision infinie a un sens physique. 
Je pense que ce concept est inacceptable. 


Manuscrit reçu le 20 décembre 1958, 
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ÉTUDE DES RAYONS Y PRODUITS PAR LA DIFFUSION INÉLASTIQUE DES NEUTRONS RAPIDES 


Par Boris MILMAN, GrorGes AMSEL et Marte-CLaupe LOYAU, 


Laboratoire de Physique de l’École Normale Supérieure. 


Résumé. 


! à ON a étudié les spectres des rayons y produits par la réaction (x, n°) dans les noyaux 
suivants : Fe, ?7A1, 1°%F, 55Mn, Cu, 5*Co, #1P, As, avec des neutrons de 2,7 MeV. 


Les transitions 


mises en évidence permettent de compléter les schémas de désexcitation de ces noyaux. 


Abstract. 


The y-ray spectra resulting from the (x, n’) reaction have been investigated for 


the following nuclei : Fe, ?7A], 29F, 55Mn, Cu, 5°Co, #P and As, for a primary neutron energy 


of 2,7 MeV. 


1. Introduction. — La diffusion inélastique des 
neutrons rapides à été utilisée par de nombreux 
chercheurs pour linvestigation des niveaux nu- 
cléaires. Nous nous sommes attachés à l’étude des 
spectres de rayons y émis par certains noyaux 
excités par des neutrons de 2,7 MeV. 

Ces expériences présentent l'intérêt d’exciter des 
niveaux jusqu’à maintenant inconnus ou de vérifier 
d'une manière nouvelle les schémas obtenus par 
d’autres réactions : (p, p'), (p, y), ete... 

Nous avons choisi des éléments isotopiquement 
purs, Ou ne présentant que deux isotopes stables, 
et qui n’ont été que peu ou pas étudiés par la dif- 
fusion inélastique : Cu, Co, #1P, 75As. 
=. Nous avons étudié, pour vérifier notre méthode, 
le Fe, *’AI, #F et 55Mn, dont les spectres sont bien 
connus. 


Il. Technique expérimentale. --- 1. LA sourcx 
DE NEUTRONS. -— Les neutrons étaient produits par 
la réaction d-d. Le faisceau de deutons était pro- 
duit par l’accélérateur de 100 KV de l'Ecole Nor- 
male Supérieure. On opérait sous une tension 
de 65 kV avec une intensité de 400 wA. La cible 
était formée par locclusion du deutérium apporté 
par le faisceau dans une feuille d'aluminium de 
0,2 mm d'épaisseur, refroidie par un jet d'air. Cette 
disposition présentait Pavantage de réduire la quan- 
tité de matière qui formait la cible ; on produisait 
ainsi un minimum de rayons y parasites. Dans ces 
conditions, nous obtenions un flux total de 106 neu- 
trons par seconde. 

Les neutrons produits dans une cible à ocelusion 
ne sont pas mono-cinétiques. On peut calculer leur 
spectre en tenant compte : 

de la densité des deutons occlus dans Palumi- 
nium à une profondeur donnée. Cette densité à 
été récemment étudiée par Fiebiger [1] ; 

de la variation du pouvoir d’arrêt de l'aluminium 
avec l'énergie des deutons : en extrapolant les 
données de Warshaw [2], et d’après d’après Jorgen- 
sen [3], nous avons trouvé, pour E < 100 keV : 


——— = 90 ?/E keV/mg/cm?° ; 
dx V 


The observed transitions allow the level schemes of these nuclei to be completed. 


de la section efficace (Æ) de la réaction D(d,n) 
He, donnée par [4]. re 
Le spectre obtenu présente un pic à 2,7 MeV 


2700 


2800 
E, keV 


F16. 1. -- Spectre des neutrons d’une cible à occlusion à 00 
( 


Ha 65 keV: 


avec une largeur à mi-hauteur de 100 keV, soit 
avec une dispersion de # %, (fig. 1). 


2. APPAREILLAGE. — Le détecteur consistait en 
un cristal de NalÏ(TD monté sur un photomulti- 
plicateur DUMONT 6292. Les impulsions passaient 
par un amplificateur ECKO 5688 et étaient 
envoyées sur un analyseur à 100 canaux du 
type SUNVIC. Des transformateurs d’impulsions 
permettaient l’adaptation d'impédance du long 
câble de 50 ohms qui reliait ces deux derniers appa- 
reils. Un groupe d’échelles rapides permettait 
d'effectuer les corrections de pertes de taux de 
comptage de l’analyseur SUNVIC, qui étaient 


‘importantes, vu le temps mort élevé de cet appa- 


reil : 800 us. 

La stabilité et la linéarité de l’ensemble étaient 
contrôlées après chaque enregistrement, à laide 
des sources suivantes : *%Hg, °?Na, 1%7Cs, 5*Mn, 
60Co et Rth. Tout spectre dont la stabilité n’était 
pas assurée à mieux que 1 % à été rejeté. La réso- 
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lution du détecteur était de 9 % pour la raie 
de 662 keV du 1#7Cs pour chacun de nos cristaux. 


3. GÉOMÉTRIE. —— Les neutrons produisaient 
dans le détecteur un bruit de fond important ; 
pour rendre optimum le rapport signal/bruit, nous 
avons utilisé : 

un blindage de tungstène de 20 cm qui protégeait 
le détecteur contre le faisceau direct des neutrons ; 
ce blindage réduisait le bruit de fond d’un facteur 
de cinq, et rendait son spectre plus uniforme ; 

des cristaux dont l'épaisseur était comparable à 
un parcours dans le Nal des rayons y de l’énergie 
moyenne de la bande explorée ; 

des diffuseurs épais également d’un parcours, 
dans la direction du cristal. 

Des expériences préliminaires nous ont permis 
de déterminer ces deux dernières conditions. 

Les diffuseurs étaient appliqués contre la face du 
cristal et avaient le même diamètre que ce dernier. 
Cette disposition est représentée figure 2 ; elle nous 


50 mm be 
DÉAMEREE 


200mm : 


| 


FE x 
j Nal(Tl | jui 
UD, 1/82 
F1. 2, — Géométrie avec le grand cristal. 


permettait d'obtenir un rapport signal/bruit de 
l’ordre de 1 avec des masses de diffuseur de l’ordre 
de 100 gr. 

Les cristaux choisis avaient un diamètre 
de 45 mm et une épaisseur de 50, 25 et 2 mm. Le 
cristal de 2 mm était employé à des énergies infé- 
rieures à 200 keV. Dans cette région, l'efficacité de 
détection du cristal pour les rayons y est bien supé- 
rieure à celle pour les neutrons et le blindage 
devient inutile. Pour évaluer le bruit de fond, on 
filtrait les rayons y par une feuille de plomb 
de 2 mm, selon une méthode dite «en sandwhich », 
due à Rothnmian [2]. En gardant la quantité de 
matière placée dans le faisceau constante, comme 


No 


l'explique la figure 3, on élimine pratiquement 
l'effet des neutrons dans les spectres. 


| ESS -— Diffuseur 
V2 +— Pb 
SSSSS + Na (TU) 

Position pour le bruit 

- de fond r 


| us +— Pb 
UTP — Diffuseur 
SSSS<— Na i (TD) 


Position pour les 
rayons Y 


Fic. 3. — Géométrie avec le cristal de 2 mm. 


4. L'EXPÉRIENCE, -— Pour obtenir les spectres 
des rayons y avec les cristaux épais, on effectuait 
des mesures avec et sans diffuseur, et on calculait 


la différence : 
ÿe —= A Ne B 


où y représente la contribution du diffuseur, À et B 
les spectres avec et sans le diffuseur, pris avec le 


.même nombre total de coups, et « un facteur de 


normalisation. «B représente le spectre de bruit 
normalisé au flux de neutrons utilisé dans la 
mesure, compte tenu des pertes de taux de comp- 
tage de l’analyseur. Les taux de comptage étaient 
de l’ordre de 200 coups/s. On opérait d’une manière 
analogue pour le cristal mince. 


5. DÉPOUILLEMENT. — Un spectre ainsi obtenu 
contient le spectre cherché des rayons y superposé 


au spectre produit dans le détecteur par les neu- 


trons diffusés élastiquement et inélastiquement par 
le diffuseur. Ce bruit de fond inhérent au phéno- 
mène est faible, vu l'efficacité de détection réduite 
du cristal pour les neutrons, comparée à celle pour 
les rayons y. Néanmoins, les pics caractéristiques 
de l’I, du Na et de l’AI (qui enveloppait les cris- 
taux) pourraient apparaître dans la différence. A 


chaque fois que nous observions un pic dont l’ori- 


gine était en doute, nous faisions des expériences de 
vérification. Quatre méthodes ont été utuisées : 

a) Comparaison des spectres de plusieurs élé- 
ments de poids atomiques voisins. Si la raie en 
question était retrouvée dans tous ces spectres avec 
la même intensité. elle était attribuée aux neutrons. 
Nous avons préféré cette méthode à celle employée 
par d’autres chercheurs (Day [6], etc...) qui esti- 
maient l'effet des neutrons en soustrayant la contri- 
bution d’un diffuseur en graphite. En effet, ni 
l’énergie, ni la distribution angulaire (Remund [7]) 
des neutrons diffusés élastiquement ne sont compa- 
rables d’un corps à l’autre, si leur poids atomique 
diffère sensiblement. 

b) A des énergies mférieures à 300 keV, une expé- 
rience d’absorpuon qui consistait à intercaler des 
épaisseurs croissantes de Pb entre le diffuseur et le 
cnstal. Un pie dû aux neutrons n’était prati- 
quement pas modifié, pour des épaisseurs de Pb 
absorbant complètement des rayons y jusqu’à 
300 keV (voir le cas de 75As). 
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c) Des rayons y d'énergie supérieure à 4,5 MeV 
peuvent facilement être détectés au moyen d’un 
scintillateur plastique, sans bruit de fond dû aux 
neutrons. La précision d’une telle mesure est faible, 
mais le doute quant à l'existence même du rayon 
peut ainsi être levé (voir le cas du ?1P). 

d) Nous avons remarqué que le quotient des 
spectres avec et sans diffuseur, tout en mettant en 
évidence les pics dus aux rayons y, atténuait ceux 
dus aux neutrons. En effet, le spectre dû aux neu- 
trons est sensiblement proportionnel au bruit de 
fond, dû également aux neutrons (voir le cas 
du 55Co). Chaque spectre est la somme de quatre 
spectres enregistrés après vérification de la stabi- 
lité. Deux expériences indépendantes ont été exé- 
cutées, au moins, dans chaque cas. Leurs résultats 
ont toujours été concordants. 

La pureté chimique des échantillons était assurée 
à mieux que 99,5 %. Ils étaient contenus dans de 
minces boîtes en papier. 

Tous les calculs, ainsi que le passage du système 
de numération binaire, utilisé par le SUNVIC, au 
système décimal, ont été exécutés par l’ordi- 
nateur 704 de I. B. M., dont l’emploi nous a permis 
d'économiser beaucoup de temps. 
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IIT. Résultats. — 1. ÉLÉMENTS ÉTUDIÉS À TITRE 
DE VÉRIFICATION. — Les spectres du Fe, AI, F, Mn 
montrent la présence de rayons y aux énergies 
suivantes, exprimées en keV : 


Fe 27 A] 107 55Mn 
840 + 20 (56Fe) 860 +20 110 +& 128 +5 
1230 + 25 (56Fe) . 1 035 +-20 200 + 6 
1 420 + 25 (54Fe) 2 240 + 30 


Fe et Al ont été étudiés avec le cristal de 50 mm, 
F et Mn avec celui de 2 mm. 

Toutes ces valeurs sont en bon accord avec les 
résultats de Beghian, Hicks et Milman [8]et Day [6] 
qui ont étudié ces éléments avec la réaction 
(n,n', y) les neutrons ayant une énergie de 2,6 MeV. 


2. CUIVRE. — Les rayons y détectés ont les 
énergies suivantes, exprimées en keV : 

E Iso- CRISTAL E Iso- CRISTAL 
TOPE MM TOPE MM 
360 +10 65 2 14925 +95 68 50 
660 +15 63 25—50 1540 +25 68 50 
770 420.65 12550. 1.860-+ 30. 63 50 
965 +20 63 50 2 060 + 30 50 
1420 +20 65 50 2 300 + 30 50 
4 340 +25 63 90 2 460 + 35 20 
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Fic. 4. — Spectre du Gu, cristal de 50 mm, 27 keV/canal. 


a) Région au-dessus de 500 keV (fig. 4). — Le 
diffuseur avait une épaisseur de 11 g/em? jusqu’à 
1 600 keV/canal de 27 g/em? au-dessus. Dans la 
région de recouvrement, vers 1 400 keV, on peut 


observer que la résolution avec le diffuseur mince 
est meilleure qu'avec le diffuseur épais, conformé- 
ment à nos conelusions quant aux épaisseurs 
optima des diffuseurs, 


QX 
RS 


Tous les pics observés sont clairement dus aux 
rayons y du Cu. Le pic de 660 keV se décompose en 
deux pics : 660 15 keV et 770 + 20 keV. 


b) Région au-dessous de 500 keV (fig. 5). — Le 
diffuseur avait une épaisseur de 8 g/em?. 


+-10000 
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F1G. 5. — Spectre du Cu, cristal de 25 mm, 11 keV/canal. 


À 360 keV, on trouve un pic assez faible, mais 
confirmée par plusieurs expériences. 

Le spectre de la figure 8 contient encore un pic 
de 210 keV, que l’on retrouve dans tous les spectres 
et que nous avons attribué aux neutrons, ainsi que 
des pics de 410 et 500 keV qui sont probablement 
les pics Compton des rayons de 660 et 770 keV. 
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Fi. 6. — Schéma du Cu. 
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Le cristal de 2 mm n’a permis la détection 
d'aucun rayon Y. 

c) Schéma (fig. 6). — L'attribution des rayons y 
aux deux isotopes du Cu ($%Cu 69 %, 55Cu 31 %) 
a été facilitée par les études sur isotopes séparés de 
Sinclair [9]et Mazari [10]. Ceci a conduit au schéma 
de la figure 9, où les valeurs des niveaux sont celles 
déterminées avec précision par Mazari [10] avec la 
réaction (pp’) et les valeurs des transitions sont 
celles trouvées par nous. 

Les résultats de Day [6](n,n’, y)et Cranberg [11] 
(n, n’) confirment les nôtres. 

Seul le rayon y de 360 keV n’est pas une tran- 
sition fondamentale. Parmi les cascades possibles, 
nous proposons la cascade entre les niveaux 
de 1 482 et 1 114 keV du 655Cu. En effet, cette 
transition a déjà été observée par Siegbahn [12] 
dans la décroissance du S5Ni, en compétition avec 
la transition fondamentale de 1 482 keV, qui a la 
même probabilité. Le rayon y de 1 482 keV n’a pas 
pu être observé, son intensité étant très faible. 
L'efficacité élevée du cristal aux basses énergies a 
permis l'observation du rayon y de 360 keV. Celui-ci 
est si peu intense que l’on peut être sûr que le 
niveau de 1 114 keV a été surtout directement 
excité par les neutrons. 


3. COBALT. — Des rayons y ont été détectés aux 
énergies suivantes. exprimées en keV : 


E CRISTAL MM E CRISTAL MM 
260 + 10 20 1 480 + 25 50 
380 + 10 25 1 750 + 30 50 
590 +15 25 50 2 100 + 30 50 
685 +15 25 — 50 2 500 + 35 50 
1 200 + 20 0 
a) Région au-dessus de 500 keV (fig. 7). — 
L’épaisseur du diffuseur était de 10 g/em?. 
pan 
5 | es 
FES 
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sut Gil 
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10 ® “30 Canaux ‘50 à ad à 
Fi, 7, — Spectre du Co, cristal de 50 mm, 27 keV/canal. 
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L'existence des rayons y de 550 et 685 keV a été 
confirmée par toutes nos expériences, avec les deux 
cristaux (fig. 7 et 8). 


— 


Nombre de coups canal 


UE 
a 
S 
=) 
= 


Quotient des Zspectres 
avec et sans 
diffuseur 

2€0 380 
l % 
: cg 
MLRRTS 


res 10 20 50 
" . L- 1 + 


Canaux 40 50 
ñ ENarERe 


1G. 8. — Spectre du Co, cristal de 25 mm, 11 keV/canal. 


La forme du pic de 1 200 keV suggère qu’il est 
triple. Une décomposition a été effectuée sur un 
autre spectre, avec 10 keV/canal, enregistré pour 
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Fi. 9. — Schéma du Co. 


avoir plus de détails. Elle donne trois rayons y de 
1 200 + 30, 1 200 + 20 et 1 300 + 30 keV (fig. 7). 
Le pie de 1 480 keV est dû à un seul rayon Y, les 
niveaux voisins de cette énergie dont l’existence 
est connue par ailleurs, ne sont que faiblement 
excités (-< 10 %). Les pies de 2 100 et 2 500 keV 
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sont multiples ; leur décomposition est impossible, 
vues les mauvaises statistiques. 


b) Région au-dessous de 500 keV (fig. 11). — Le 
diffuseur était épais de 7 g/cm?. 

Le pic de 210 keV est dû aux neutrons : il dis- 
paraît d’ailleurs dans le quotient (fig. 11). 

Le cristal de 2 mm n’a permis la détection 
d’aucun rayon y. | 


c) Schéma (fig. 9). — Alors que dans la décrois- 
sance P— du Fe, seuls les niveaux de 4 100 et 
1 290 keV sont excités [13], nous avons excité la 
plupart des niveaux déterminés par Mazari [14] 
(p, p'), qui sont représentés sur la figure 12. 

Pour le triplet, vers 1 400 keV, nous avons 
conclu que seul le niveau de 1 479 keV est excité 
de manière détectable. En effet : cette raie est 
unique (voir plus haut); la transition de 
380 + 10 keV ne peut être qu’une cascade entre 
les niveaux de 1 479 et 1 097 keV, ce qui prouve 
que c’est le niveau de 1 479 keV qui est excité, et 
non celui de 1 458 keV. Ceci permet également de 
situer le rayon de 260 + 10 entre les niveaux 
de 1 473 et 1 479 keV, et non entre ceux de 1 458 
et de 1 189 keV. La cascade entre les niveaux 
de 1 289 et 1 097 keV vue dans le Fe [13] fournit 
un rayon y de 191 keV ; il a été masqué par le pic 
de l’I, contenu dans le bruit de fond. Enfin, les 
rayons de 550 et 685 keV ne pouvaient être placés 
que comme indiqué figure 9. 


d) Excitation directe des niveaux. — On peut se 
demander si les rayons y de 1120, 1 200 et 
1 480 keV proviennent de niveaux directement 
excités par les neutrons, ou s’ils sont dus seulement 
aux cascades partant des niveaux supérieurs et qui 
fournissent les rayons y de basse énergie que nous 
avons observés. Ces derniers sont peu intenses, et 
on peut, dans chaque cas, être sûr que les rayons y 
en question proviennent surtout de niveaux direc- 
tement excités par les neutrons. 


4. PHosPHORE. — L'étude du spectre de %1P, 
faite avec un diffuseur de 2,5 g/em?° (la faible densité 
du phosphore, (0,8), nous empêchait d'en accu- 
muler davantage utilement), et le cristal de 50 mm, 
met en évidence deux rayons y d’une énergie de 
4 270 + 25et 2 230 + 30 keV. 


a) Discussion du spectre (fig. 10). — Le pic de 
1 270 keV est clairement dû au *P. Le pic de 
2 230 keV se superpose exactement au 2 240 keV 
de l’AI (fig. 14a). Pour nous assurer que ce pic est 
bien dû au #:P, nous avons fait une expérience avec 
un scintillateur organique NE 102, 40 X 50 mm 
(Nuclear Enterprise, Édimburgh), non contenu 
dans de l'aluminium. Tout aluminium a été éliminé 
du faisceau. Le spectre de la figure 116 confirme 
l'existence de ce rayon y. Le bruit de fond est ic 
très faible, car le spectre dû aux neutrons 
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F1G. 10. — Spectre du P, cristal de 50 mm, 27 keV/canal. 
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de 2,7 MeV dans ce scintillateur a la forme d’un 
spectre dû à des rayons y de 4 MeV. 

Remarquons que le niveau de 2 230 keV appa- 
raît très bien sur le quotient (fig. 12). Le pic 
vers 4 700 keV est le pie «paire + 1 » du 2 230 keV. 
Les deux pics vers 600 keV (fig. 13) semblent pro- 
venir des neutrons, l’allure des courbes de la 
figure 14 dans cette région, ainsi que le quotient 
(fig. 12), le suggèrent nettement. 

Nous avons recherché avec une attention parti- 
culière la cascade entre les niveaux de 2 230 et 
4 270 keV qui donnerait un rayon y vers 1 000 keV. 
On voit (fe. 11) le pic important de l’Al vers 
1 000 keV qui se répète dans les spectres avec et 
sans diffuseur. Nous avons d’autre part comparé 
le spectre du rayon de 1 270 keV (fig. 10) à celui 
de la source ??2Na qui émet un rayon y d’énergie 
voisine : 1 280 keV. Ce spectre, également repré- 
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Fic. 11, — a) Spectres avec et sans diffuseur du P. 
b) Spectre du P avec un scintillateur organique. 


senté figure 10, montre que le pic observé vers 
1000 keV est le pic Compton du rayon de 
1 270 keV. Nous en avons conclu que la cascade a 
une intensité inférieure à 5 % de celle de la tran- 
sition fondamentale de 2 230 keV. 


b) Schéma. — L'accord est satisfaisant avec les 
expériences de (p, p') [15], (n, n’) [16] (qui n’a 
mis en évidence que le pie de 1270 keV 
#S1(p, y) [17] et #1S1(8—) [18] et infirme les résul- 
tats de Scherrer [19]. 

L’absence de la cascade est également confirmée 


par d’autres auteurs, comme le montre le schéma 
publié par Endt et Braams [17]. 
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Ces résultats sont en contradiction avec ceux 
de Boley, etc. [20] (qui ont paru après l’achè- 
vement de nos expérience). Ces auteurs travail- 
laient sans blindage et le rapport signal/bruit s’en 


trouvait réduit. Il semblerait que ce soit là la 
cause de notre désaccord. 
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Fic. 43. — Spectre de l’As, cristal de 50 mm, 27 keV/canal. 


5. ARSENIC. — Le’5As étudié par les trois cris- 
taux et avec des épaisseurs de 6,5 et 4 g/cm?, a 
donné les spectres des figures 13, 14 et 15. Voici les 
rayons y trouvés, leur énergie étant exprimée en 


keV : 


ÉTUDE DES RAYONS y 57 
No E CRSTAL MM N: E CRISTAL M 
VIS DNCERRS 2 de 575 +15 25 —50 
AE SOC AP EE RTE EE) 50 
DD ST TON D 100 Ye A 0/0 NE 20 50 
0010 460 25 50, y, 43504 25 50 


a) Discussion des spectres. — Les pics de 135, 
410, 575, 1 070 et 1 350 keV sont clairement dus 
à 75As, ainsi que ceux de 275 et 820 keV qui sont 
doubles. 

Le pic de 820 keV se décompose en deux pics 
de Y7 — 810 Æ 25 keV et y, — 870 + 25 keV, 
comme le montre aussi le quotient (fig. 13). 

L'énergie du pic de 275 keV était connu à 2 keV 
près par comparaison au pic du *#Hpg, de 279 keV. 
Connaissant l’existence des niveaux de 265 et de 
280 keV de 75As [21, 221, nous avons pu décom- 
poser ce pic et attribuer à ces rayons les intensités 
respectives de 2/3 et 1/3 (fig. 14). 
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Fic. 14. — Spectre de l’As, cristal de 25 mm, 11 keV/canal. 
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Fic.,15. — Spectre de l’As, cristal de 2 mm, 4 keV/canal 
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Le pic à 200 keV se superpose au pic de 210 keV 
de VI et apparaît dans tous les spectres. Comme il 
était très intense, nous avons pensé qu’il provenait 
partiellement de 75As. Pour le prouver, nous avons 
effectué une expérience d'absorption par le Pb avec 
un eristal de 12 mm. Sur la figure 19 (échelle 6 keV/ 
canal), il apparaît que le pic de 200 keV, contrai- 
rement à celui de 275 keV, ne s’atténue pas expo- 
nentiellement. On voit que ce pic est la superpo- 
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F1G. 16. — Spectre d'absorption de l’As par Pb, 
6 keV/canal. 


sition de deux pics d’intensités comparables, l’un 
dû aux neutrons, à 210 keV, l’autre dû aux rayons y 
de 200 keV. Ceci est confirmé par le décalage pro- 
gressif du pic de 200 keV vers 210 keV avec les 
épaisseurs croissantes de Pb. La constance du gain 
est démontrée par le pic de 275 keV dont la posi- 
tion est restée fixe. 


b) Schéma (fig. 17). — Le capture X du 75Se [21] 
confirme la présence de Y:, Y2, Ya: Ya Ys- 

Par excitation coulombienne [23], ÿ4 Ye et Y: 
ont été trouvés. 

Les rayons Vs, Yo et Y sont inconnus. Nous 
suggérons que le rayon de 870 keV provient d’une 
cascade du niveau de 1 070 à celui de 200 keV. 

Remarquons que les niveaux de 470 et 620 keV 
mis en évidence par l’étude de 7°Ge et le niveau 
de 304 keV connu par 75%$e [21, 22] n’ont pas été 


ou très peu excités. 


N°1 


Les valeurs des niveaux indiquées sur le schéma 
de la figure 17 sont tirées de Van den Bold [21]. 


Fic. 149. — Schéma de l’As. 


c) Excitation directe des niveaux. — Van den 
Bold, etc. [21] donnent les intensités relatives des 
cascades qui partent du niveau de 401 keV de 75As. 
Nous avons calculé les intensités relatives dans 
notre cas avec assez de précision pour pouvoir 
faire des comparaisons : ù 


E VALEURS DE NOS VALEURS PROPORTION 
VAN DEN Born D’EXCITATION 
DIRECTE 
401 1 1 1 
136 97 5 CA 
265 27 6 1/2 
280 1 3 2/3 


Nous n’avons pas pu étudier complètement la 
cascade entre les niveaux de 401 keV et de 280 keV, 
car le rayon y de 121 keV correspondant n’a pas 
pu être séparé du rayon y de 136 keV. 
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LETTRES A LA RÉDACTION 


ÉTUDE 
PAR RÉSONANCE MAGNÉTIQUE NUCLÉAIRE 
DE L'EAU FIXÉE 
SUR UN GEL DE SILICE 


Par M. Kan-ichi Kamiyosui, 
Université de Tohoku (Japon) (1) 


Nous avons étudié la résonance magnétique de l’eau 
adsorbée sur un gel de silice (actigel Prolabo). La 
figure 1 représente la largeur AH (en milligauss) de la 
raie, en fonction de la teneur en eau (masse d’eau par 
g de gel de silice) qui est déterminé par séchage 
6 heures à 110 °C. Le poids d’échantillon utilisé pour 
la détermination de AX est d’environ 0,15 g. 

La largeur est d’environ 200 milligauss pour 10 % 
d’eau ; elle diminue graduellement jusqu’à moins de 
10 milligauss vers 70 % d’eau. Rappelons que la 
largeur pour l’eau liquide et la glace sont respecti- 
vement de l’ordre de 0,7 milligauss et 16 gauss. 


60 70 
Teneur en eau (g d'eau pour 1g) x 110 
F1G. 1. — Résonance magnétique nucléaire de l’eau fixée 


sur un gel de silice et sur le verre en poudre (largeur, 
AH, de la raie en fonction de la teneur en eau). 


Sur la figure 1 sont indiqués également les résultats 
relatifs à l’eau fixée sur le verre en poudre ; pour la 
même teneur en eau la largeur est plus faible que celle 
du gel de silice. 


(2) Travail effectué au 
Recherches Physiques à 
hertzienne ; PT FREYMANN). 
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COEFFICIENTS DE CONVERSION INTERNE 
DANS LES NIVEAUX L,, Li Lu 


Par Marcel Frizzey et Manuel VALADARES, 


Laboratoire Curie, Paris 
et Laboratoire du Grand-Aimant Permanent, Bellevue. 


Nous avons publié en 1957 [1] des résultats expéri- 
mentaux concernant les intensités relatives des raies de 
conversion dans les niveaux Lx, Zar et Larr du rayon- 
nement y de 46,50 keV émis au cours de la trans- 
mutation RaD — RaE. - 

Pour les comparer avec la théorie, nous n'avions 
alors que les calculs de Rose et Goertzel [2] qui négli- 
geaient les effets des dimensions nucléaires mis en 
relief par Sliv. 

Après la publication de notre mémoire, Rose et ses 
collaborateurs ont effectué de nouveaux calculs pour 
les couches Zx et Lrr [3], en tenant compte des dimen- 


- sions nucléaires et ont admis que celles-ci n’ont pas 


d'influence sur la conversion dans la couche Zrrr. Nous 
venons de recevoir les tables de Sliv et Band [4] con- 
cernant ces mêmes calculs pour les trois couches. 

Il nous a paru intéressant de comparer dans le 
tableau suivant les résultats de l'expérience et de ia 
nouvelle théorie. 


THÉORIQUES 
Noyau [Noyau 
s PONCTUEL NON PONCTUEL 
RAPPORTS EXPÉRIMENTAUX Rose ar Boss 
GŒRTZEL Ban 
LlLn 9,5 + 0,4 11,2 10,9 , 9,3 
Lu /Lrrx 10,5 + 0,9 14,7 11,6 41,7 


Les valeurs de Rose sont doublement interpolées : 
pour les Z entre 75 et 85 ; pour les énergies, entre 0,05 
et 0,10 m, c?. Les valeurs de Sliv et Band sont données 
pour Z — 83 et sont interpolées seulement pour les 
énergies, entre 0,05 et 0,10 m, c?. 

L’accord entre l’expérience et les calculs théoriques 
de Sliv et Band est excellent, en ce qui concerne le 
rapport Lx/Lar. 

Pour le rapport Lrr /Lrr, le désaccord d’environ 40% 
qui existait a disparu, mais la valeur théorique reste 
légèrement supérieure à la valeur expérimentale. Il 
serait pourtant prématuré de parler de désaccord, en 
effet, ce rapport est beaucoup plus difficile à déter- 
miner expérimentalement que le premier, pour deux 
raisons : la raie Lux est très faible par rapport au 
fond continu et l'intervalle énergétique Lir, Larx est 
quatre fois plus grand que celui de Lx à Lar. Il est 
possible, par conséquent, que nous ayons sous-estimé 
l’erreur commise. 


Remarquons, d’ailleurs, que ces calculs sont faits en 


N° 1 


admettant que le rayonnement + de 46,50 keV est un 
dipôle magnétique pur. En supposant qu’avec le M1 
un rayonnement E2 serait mélangé dans la proportion 
en photons de 1/1 000 seulement, les rapports Zx [Lux 
et Lu /Lurx deviendraient, d’après les tables de Sliv et 
Band, respectivement 9,2 et 11,2, en parfait accord 
avec l’expérience. 
Lettre reçue le 20 octobre 1958. 

Note ajoutée à la correction. — Nous v s der ir 
des tables de M. E. Rose (‘* Internal hat 
cients ”, éditées par North-Holland Publishing Company, 
Amsterdam, 1958), où les coefficients de conversion ont 


donnés pour Z — 83 et calculés directement pour les trois 
niveaux L. 


On peut en déduire les rapports : Lr/Lxr — 10,5 et 
Lax/Larr = 10,4, plus proches des valeurs expérimentales 
que les résultats antérieurs de Ross. 
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REMARQUE A PROPOS D’UNE MÉTHODE 
DE CALCUL DES CHAMPS MAGNÉTIQUES 


Par Pierre GAUTIER, 
Laboratoire de Physique, Faculté des Sciences, Toulouse. 


M. Laudet a proposé [1] une méthode intéressante 
pour calculer l’induction et des dérivées en tout point 
de l’axe d’une lentille électronique magnétique de révo- 
lution, lorsqu'on connaît la valeur du flux ® à travers 
un cercle de même axe que la lentille. Cette méthode 
se prête particulièrement bien au calcul avec des 
machines électroniques arithmétiques. Elle est utilisée 
au Laboratoire d’Uptique Électronique du Centre 
National de la Recherche Scientifique et donne des 
résultats très satisfaisants. Le flux ® est mesuré au 
moyen d’un montage électrodynamique réalisé par 
M. Lurandeau [2] et les calculs sont effectués sur la 
machine élecrtonique arithmétique de l’Institut de 
Calcul Numérique de la Facul.é. d 

D’après le titre. donné à l’article [1] on pourrait 
penser que ce mode de calcul, très pratique pour les 
champs de révolution, n’est valable que dans ce cas. 
Nous voudrions attirer l’attention des utilisateurs 
sur le détail suivant : la méthode esi tout à jait générale. 
Elle peut être employée sans aucune modification 
chaque fois que l’on veut connaître la répartition, le 
long d’une droite donnée, de la composante sur cette 
droite d’un champ magnétique quelconque. C’est le cas 
en parviculier dans les lentilles cylindriques, dans les 
groupes quadrupolaires. 11 sufli de placer la bobine 
. de revolution qui sert à la mesure du flux, de façon que 
son axe soit confondu avec la droite en question. Il 
a y a pas a modifier le montage expérimental. Quant 
au programme de la machine à calculer, il est établi 


LETTRES À LA RÉDACTION 61 


une fois pour toutes et peut servir pour calculer 
n'importe quel champ magnétique. 
En effet, la méthode proposée repose uniquement 
sur l’équation : 
DEnCDe) 20h 


| 120. 
d=* dr? 


Aion ) 
à laquelle obéit la fonction flux. Or, cette équation 
n'est pas, comme on le croit souvent, une caracté- 
ristique des systèmes de révolution. Quelles que soient 
les symétries du champ, le flux 


rT 27% 
ds for [Be o,dde  () 


à travers un cercle d’axe Oz et de rayon r, satisfait 
à (1). Cela est un corollaire d’un résultat que nous 
avons établi ailleurs [3]. La vérification directe en est 
simple. Le calcul des dérivées successives de ® et leur 
substitution dans (1) conduit à : 


4 HS | 
f pde f 3: Pie, 9,2) de 


On tient compte ensuite du fait que B satisfait à 
l’équation de Laplace en remplaçant, dans (3), 


_ d? B,/dz? par son expression tirée de 


ROYAN 1#02B;. 202B 
ABS ( : ns FL 5 — 
(AB): z de (Q de Î pi Do ‘ 27 


et on obtient après intégrations 
fipes 

J pIde 

ce qui est bien une identité, puisque dB./d@ reprend 


la même valeur pour @ —0 et op — 27, quel que 
soit 7. 


0, 


2T 


= 0, 


0 


Lettre reçue le 23 octobre 1958. 
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ÉTUDE DU Ho! 


Par L. GRENACS et A. MEESSEN, 


Centre de Physique Nucléaire 
de l’Université de Louvain, Belgique. 


Nous avons étudié les rayons gamma du $6Dy1$5 se 
transformant en &Hc185. Le Ly fut irradié à la pile 
de Mol sous forme u’oxyue et évuaié sous torme solide 
ainsi que sous 1orme liquide. 

Dans le Ly1%5 de 2,3 heures, nous avons trouvé, en 
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accord avec Jordan [1] l’existence des deux cascades 
suivantes : 10 630-360 keV et 20 270-710 keV. 

La première cascade révèle une nature métastable 
pour le niveau de 360 keV. En variant le temps de 
résolution de l’appareïl, il est possible de conclure que 
letemps de vie de ce niveau est supérieur à 2,5.10-7 s 
Kane [2] a confirmé cette mesure par coïncidences 
retardées trouvant un temps de vide de 6,65.10-6s. 

La deuxième cascade permet la mesure de corré- 
lations angulaires, n’indiquant pas de délai interne. 
Les résultats de ces mesures sont reportés à la figure 1. 
La méthode des moindres carrés fournit les coefficients, 
corrigés pour les angles solides, suivants : 


Ao=1, A3 —=—0,040 + 0,007, A, — 0,011 + 0,006. 


100 a) 7 


ts S) 


(0) S 
è 
à 


09, 


09 
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Dans la figure 1 nous rapportons également les 
courbes théoriques correspondant à des transitions 
pures : 


EMA 2 0) PR D)? 
(bem 2e 2 
(DEEP ES 2e 2 
(c): 3/2 (2) 5/2 (1) 9/2 
Notons encore que la mesure de l’anisotropie pour 
une source liquide : À = — 0,034 + 0,060 corres- 
pond bien à celle pour la source solide : 
A2 0;:05322:070728 


dans les Hinites des erreurs expérimentales. 

Le même corps sera encore étudié dans notre labo- 
ratoire par d’autres techniques, avant de donner une 
interprétation générale des résultats et de proposer les 
spins et parités des niveaux. 


Lettre reçue le 7 novembre 1958. 
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LA RÉPARTITION DES IMAGES LATENTES 
DONNÉES PAR DES RAYONS P 
DE DIFFÉRENTES ÉNERGIES SUR LA SURFACE 
ET A L'INTÉRIEUR DES GRAINS 
D’'HALOGÉNURE D'ARGENT 
DANS LES PLAQUES NUCLÉAIRES 


Par Shinichi KiKkucur, 
Laboratoire de Physique et Chimie Nucléaires, 
Collège de France. 


Hoerlin et Hamm [1] ont étudié la distribution au 


sein des grains d’émulsion des images latentes données 
par les rayons «, les rayons X et la lumière et montré 


que les images données par les rayons « sont beaucoup 


plus dispersées que celles données par la lumière. 
D’autres auteurs ont étudié le même problème, mais 
leurs résultats n’étaient pas toujours en accord. 
Winand [2] a dit que les images latentes produites par 
les rayons « se trouvaient dans l’intérieur des grains 
d’halogénure d’argent, tandis que Haenny et Gail- 
loud [3] prétendaient qu’elles se situaient plutôt à la 
surface. Cüer [4] a fait une série d’études pour clarifier | 
la question en exposant les plaques nucléaires aux 
rayons « du Po, aux rayons $ du RaE et y du’Ra, au 
moyen de développements superficiels et internes et de 
post-exposition à la lumière infra-rouge. D’après lui, 
les images latentes données par les rayons $ et y à 
haute énergie sont plus dispersées au sein des grains 
d’émulsion que celles données par les rayons & ou la 
lumière. 


Le présent auteur et son collaborateur [5] ont fait 


une expérience de discrimination des images latentes 
données par les rayons f émis par S-35, P-32 et 1-131 
et ont montré que le rapport des images latentes 
internes aux images latentes superficielles est plus 
grand que dans le cas d’exposition à la lumière. Dans 
ces expériences ils ont choisi des radio-isotopes comme 
source de rayons £ de différentes énergies, mais ils 
n'ont pas défini leur énergie. L’auteur a donc entreprisle 
premier, indépendamment de Kartuzhanskij et Soltit- 
Ski] [6] qui ont étudié la sensibilité des émulsions 
photographiques aux rayons f de diverses énergies, la 
présente recherche pour préciser l’effet des rayons £ 
d'énergies définies, sur la répartition des images 
latentes au sein des grains en employant un spectro- 
mèêtre-f. 


PROCÉDÉS EXPÉRIMENTAUX., L'auteur se servant 


d’un spectromètre $ magnétique à 1800 [7], aexposé les 
plaques nucléaires d’Ilford G-5 dont l'épaisseur d’émul- 
sion est de 50 {2 aux rayons 6 de différentes énergies 
dans le domaine de 50-500 keV. En tant que sources 
de rayons $ il s’est servi de S-35 pour les énergies 
de 50-120 keV et de P-32 pour celles de 150-500 keV. 
L'intensité des sources était de l’ordre de 4 mC. Le 
porte-plaque dont le schéma se trouve dans la figure 1 
est émis à la place du détecteur du spectromètre. La 
plaque de T em X 1 em reçoit le faisceau de rayons 8 


ayant la forme du diaphragme, X-Y, qui est représenté 
à gauche sur la même figure. Devant ce diaphragme du 
porte-plaque sont appliquées deux couches super- 
posées de formvar dont l’épaisseur totale est de l’ordre 
de 90 1g/cm?. Les sources sont déposées sur un support 
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de LC 600 de même épaisseur que la fenêtre du porte- 
plaque. Le porte-plaque et l’intérieur du spectromètre 
sont mis sous vide pendant l'exposition des plaques 
photographiques. Le temps de pose est choisi de telle 
sorte que l’intensité totale de rayons 8 que les’ plaques 
reçoivent, soit identique dans tous les domaines 
d'énergie. Le temps de pose se répartit entre 1 heure et 
quelques heures. 

Les plaques ont été développées par les révélateurs 
superficiel et interne dont les formules et les temps de 
traitement sont les suivants : 


RÉVÉLATEUR SUPERFICIEL 


RÉVÉLATEUR INTERNE 

Glycin 2,5 g 

K:C03, 20 g 

Eau q.s. p. f., 80 ml 

10 ml de solution alcoolique à 

0,5 % de 6-nitrobenzimidazole 

Ajouté comme antivoile 

Développement, 10 mn à 20 °C 

Bain d’arrêt de 2% d’acide acé- 
tique 10 mn 


Génol, 3 g 

Na:SO;, 50 FX 
Hydroquinone, 6 g 
Na,CO:, 25 g 

KBr,3g 

Eau q.s. p. f., 1 000 ml 


10 mi de solution d’anti- 
voile ajouté à 80 ml du 
révélateur. 

Après l’exposition. 

Oxydation des images la- 
tentes superficielles par 
la solution aqueuse de 
3 % de K,Cr,0, pen- 
dant 10 mn. 

Lavage à l’eau courante 
10 mn. 

Développement interne 
pendant 6 mn à 20 °C. 

Bain d’arrêt et fixage. 


Fixage par fixateur acide. 


Vide 


F6. 4. — Le porte-plaque. 


La densité du noireissement des plaques est mesurée à 
l'aide d’un microphotomètre (1) le long de l'axe X-Y de 
la figure 1. Le microphotogramme est schématisé 
. dans la figure 2 où l’abscisse donne la direction de 
déplacement du microphotomètreet l’ordonnée, l’'inten- 
sité de lumière transmise par la plaque développée. 
y est la lumière transmise par le voile seulement et 7 
est celle transmise par l’image et le voile. Pour obtenir 
l'intensité totale de rayons B reçue par la plaque 
l’auteur à mesuré les surfaces Sy et S et il a calculé 
log { (S + Sy)/Sv } qui est proportionnel à l'intensité 
totale cherchée (fig. 3). La courbe 1 présente le noircis- 
sement obtenu par le développement superficiel tandis 
que la courbe 2 montre celui obtenu par le dévelop- 


(1) La mesure à été faite grâce à l'amabilité du 


Pr E,. Vassy. 


LETTRES À LA RÉDACTION 63 


pement interne. L’absorption et la diffusion dans les 
membranes du diaphragme du spectromètre ou du 
porte-plaque ont été négligeables. 


Noircissement d'une plaque 


Marche du microphotomètre 


FiG. 2, — Microphotogramme d’une plaque : 
Sy surface de densité du voile. 
S surface de densité au-dessus du voile. 


ages 


Jatentes 


4 


Intensite relative des im 


0 50 8010020 10 .200. 300 400  keÿ 500 
Énergie des rayons P 
Fic. 3. — Répartition des images latentes suivant 


l'énergie : 1) Développement superficiel. 2) Dévelop- 
pement interne. ; 


INTERPRÉTATION DES RÉSULTATS. -- l’expérience 
montre que les rayons $ donnent les images latentes 
superficielles et internes dans les grains d’halogénure 
d'argent. Le rapport des intensités de ces deux images 
n’est pas très différent dans cet intervalle d'énergie de 
celui obtenu pour la lumière bien que les résultats pré- 
cédents [5] aient indiqué que les images internes par 
les rayons B étaient plus intenses que dans le cas de 
l'exposition à la lumière. Les rayons f dont l’énergie 
est inférieure à 450 keV donnent beaucoup d'effets sur 
les émulsions. Ce phénomène indiquerait que ces 
rayons 8 se trouveraient près de la fin de leurs parcours 
dans la couche d’émulsion de 50 {1, en accord avec les 
résultats précédents [8]. Dans tous les domaines 
d'énergie et dans les conditions de développement de ce 
présent travail, les images latentes superficielles sont 
plus denses que les images latentes internes. La seule 
anomalie que l’on note pour les rayons £ du 50 keV 
serait due à la fluctuation de l'intensité totale des 
B reçus sur les plaques. Dans cette région d'énergie, 
une faible variation du courant dans la bobine électro- 
magnétique à une grande influence sur l'intensité des 
rayons @ du S-35. La diminution sensible de densité 
avec l'accroissement de l'énergie semble provenir de la 
diminution du pouvoir d’ionisation des rayons £ à 
haute énergie. 

Qu'il me soit permis de saluer respectueusement la 
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mémoire du Professeur F. Joliot qui a bien voulu me per- 
mettre de travailler dans son laboratoire de Physique et 
Chimie Nucléaires du Collège de France et n’a pas 
cessé de m’encourager. Je suis reconnaissant à 
Mlle Yuasa de son amabilité. Elle m’a suggéré l’usage 
du spectromètre 6 et m’a guidé dans ce travail. Je ne 
veux pas oublier de remercier les membres du Labora- 
toire de Physique et de Chimie Nucléaires et M. Gauvin 
de l’Institut du Radium de leur aimable collaboration 
et de leurs suggestions utiles. 


Lettre reçue le 14 novembre 1958. 
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CALCULS APPROCHÉS RELATIFS 
A L'ÉTAT FONDAMENTAL DE Hel 


Par Mme M. DORNSTETTER, 
G. Munscuy et P. PLUVINAGE, 


Institut de Physique 
de la Faculté des Sciences de Strasbourg. 


Dans une « Lettre à la Rédaction » précédente [1], 
était donné le résultat inattendu d’un calcul relatif à la 
valeur théorique de l’énergie de l’état fondamental 
de HeZ. Schwartz [2] a publié ensuite sur cette même 
question un résultat contradictoire qui est apparu 
immédiatement comme plus vraisemblable [3]. Entre- 
temps nous avions pu établir que les calculs numé- 
riques effectués sur les données utilisées étaient exacts 
mais qu’une erreur au moins s'était glissée dans les 
données elles-mêmes. La méthode de résolution de 
l'équation de Schrüdinger n’est pas en cause ; 1l s’agit 
seulement d’une erreur matérielle de transcription à 
laquelle d’ailleurs s’en est ajoutée une autre, moins 
grave. Dans une communication privée, C. L. Pekeris, 
Directeur du Lépartement de Mathématiques Appli- 
quées de l’Institut Weizmann à Rehovot (Israë:), nous 

-a fait part des résultats d’un calcul basé sur la méthode 
de variation simplifiée par la condition € = |/— #, 
c’est-à-dire rigoureusement équivalent au nôtre. 
Pekeris est en accord avec Schwartz et pour obtenir 
une conclusion définitive nous avons repris de notre 
côté tous les calculs en utilisant le système des poly- 
nômes orthogonaux qui, dans notre méthode, remplace 


les puissances successives des variables. La racine de 


f 


l'équation en e est 
£13 —= 1,704 003 3. 


D'où la valeur de l'énergie Æ,4 — — 2,903 627 1 


unités atomiques. 


Cette valeur est exactement celle de Pekeris. Les n 
coefficients du développement de la fonction d’onde » 
suivant les fonctions Umnp SOnt, avec les notations w 


de {1} 
%ooo — 1 &101 — — 0,001 3327856 
100 = — 0,043 953582 por = -— 0,000 525 545 
LR == 0,019 1436 528 Lx30 — 0,001 061 156 
&oo1 = — 0,053 585 319  asor — — 0,000 218 189 
&ooy — — 0,005 428 807 . &300 — — 0,000 173 329 
ag — 0,007 748 497 agi == = 0,000 024 7872 
&00 — — 0,001 182 105 x 


Pour avoir une idée de la différence entre lPapproxi- 


mation réalisée ainsi et celle qui serait obtenue par la | 


méthode de variation complète, où & est regardé comme 
un paramètre non lié à E, nous avons calculé l’énergie 


cinétique moyenne J'»t l'énergie potentielle V. Dans la 


seconde méthode, le théorème du viriel est rigou- 


reusement vérifié. 


4 dr = 
= V=—T 
Ici nous trouvons 
TT = 2,902 996 9 
LP — — 2,903 319 0. 


Le défaut de vérification porte sur le cinquième 


chiffre, et il semble que 1' et V tendent vers leurs 
limites par valeurs inférieures en module. Rappelons 
que la valeur Ex de Kinoshita, la meilleure des pré- 
visions théoriques connues pour l’énergie (*), est 


Er — — 2,903 722 5. 


._ La différence avec la valeur Æ;,, porte aussi sur le. 
cmquième chiffre mais elle est environ trois fois plus : 


faible que le défaut de vérification du théorème du 
viriel. Il paraît donc probable que l’étude de ce défaut 


fournit un moyen assez sensible de déterminer dans - 


quelle mesure on s’approche de la limite théorique. 
Le développement de la fonction d’onde réarrangé 
suivant les puissances des variables nous donne 


Ÿ — 1— 1,907 308 s + 0,444 977 re 


Les coefficients théoriques sont 2 et 0,5. La fonction 


à 38 termes de Kinoshita est meilleure, comme il 
allait s’y attendre après la découverte de l'erreur. 


Lettre reçue le 15 novembre 1958. 
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(#) Note ajoutée au moment de la correction des épreuves. 
Une prépublication de PEKeris indique la valeur 


Ep = — 2,903 724 3. 
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